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1,  INTRODUCTION 


The  ban  on  open-air  nuclear  testing  forces  the  U.  S.  Army  to  look  for  alternative  testing 
techniques  in  its  effort  to  qualify  military  equipment  as  nuclear  survivable.  There  are  currently 
two  techniques  used  to  simulate  the  blast  effects  of  nuclear  explosions  involving  either  high 
explosives  in  large  quantities,  or  special  shock  tubes.  The  simulation  of  nuclear  blasts  with  high 
explosives  (HE)  is  very  costly  and  limited  to  small  yields  ( <20  kt).  HE  tests  require  much  space 
and  their  set-up  is  very  time  consuming.  For  these  reasons,  the  use  of  specialized  shock  tubes  is 
becoming  increasingly  attractive.  Such  facilities,  called  large  blast-wave  simulators  (LBS),  are 
large  enough  to  accommodate  full-sized  tactical  equipment  such  as  trucks,  tanks  and  helicopters. 
A  few  LBS  facilities  exist  abroad,  the  largest  at  the  Centre  d’Etudes  de  Gramat  (CEG), 
France/*  The  U.  S.  Army,  in  concert  with  other  government  agencies,  is  presently  developing 
a  concept  of  such  a  facility  suitable  to  simulate  both  thermal  and  blast  effects  of  nuclear 
explosions  for  the  survivability  testing  of  military  equipment  and  for  research  studies. 

During  the  concept-phase,  questions  are  raised  about  the  necessary  size  and  the  expected 
performance  of  such  a  facility,  and  to  answer  these  questions,  a  parametric  study  was  initiated. 
A  numerical  study  seemed  preferable  to  an  experimental  one,  because  the  latter  would  take  up 
too  much  preparation  time,  and  yield  too  few  data  points.  Among  the  many  numerical  fluid- 
dynamic  schemes,  the  implicit  Beam-Warming  method"  stood  out  for  its  apparent  flexibility  and 
numerical  stability.  Using  this  scheme,  we  developed  a  quasi-one-dimensional  (Q1D)  Eulerian 
code  for  simulating  the  flow  in  shock  tubes  with  arbitrary  cross-sectional  areas.  In  this  report, 
the  BRL-Q1D  code  and  its  underlying  theory  are  described.  Computational  results  are 
compared  with  experiments  to  establish  confidence  in  the  performance  and  the  applicability  of 
the  code. 


2.  THEORETICAL  CONSIDERATIONS 


The  BRL-Q1D  code  incorporates  two  computational  techniques.  One  is  an  implicit 
finite-difference  technique  developed  by  Beam  and  Warming,5'^  the  other  is  an  explicit  finite- 
difference  technique  according  to  MacCormack.5  These  techniques  are  applied  to  the  quasi-one- 
dimensional  Euler  equations  in  their  weak  conservation  form.  The  weak  conservation  form  is 
retained  as  the  Euler  equations  are  transformed  to  a  uniform,  computational  grid.  Central 
spatial  differencing  casts  the  difference  equations  into  a  block-tridiagonal  structure  which  is 
solved  for  the  increments  in  the  dependent  variables  at  each  successive  time  step. 

2.1  GOVERNING  EQUATIONS. 

The  differential  Euler  equations  which  describe  the  one-dimensional  flow  with  variable 
area  may  be  written  in  the  following  form. 

-M.+  ^L  +  n=  o,  (2.i) 

dtdz 


*  References  are  listed  at  the  end  of  the  report 
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where  the  vectors  %  %  and  /?  are 


p 

pu 

0 

?  = 

pu 

e 

,  B  = 

(pu2  +  p) 

.  «  (e  +  p). 

,  8  = 

dA 
-  P~^— 
dx 

0 

This  set  of  three  scalar  equations  represents  the  conservation  of  mass,  momentum,  and  energy, 
per  unit  volume,  with  the  usual  notation  of  p  as  density,  u  as  velocity,  e  as  total  energy,  and  p 
as  pressure.  The  cross-sectional  area,  A,  may  vary  with  x  and/or  t,  where  x  is  a  linear 
dimension  and  t  is  time.  As  written,  the  equations  are  in  weak  conservation  form  because  of  the 
vector  $.  If  A  ^  f  (x),  M  vanishes  and  the  equations  revert  to  the  strong  conservation  form. 
Either  form  is  shown  by  Peyret  and  Viviand^  to  capture  shocks  accurately  in  the  grid. 


The  physical,  independent  variables,  x  and  t,  are  transformed  into  a  uniformly-spaced 
computational  grid  by  a  general  transformation  of  the  form 


,  and  T  =  t  -  t0  , 

(2.3) 

ion  (2.1)  is  then 

dB  r  _ 

- 1-  h  =  0  , 

d* 

(2.4) 

dr 

where  it  is  noted  that  the  weak  conservation  form  is  retained.  In  Equation  (2.4)  we  define  the 
transformed  vectors  as  follows. 


Q  =  = 


pA 

puA 

eA 


pAit  +  puAZt 

0 

2-  IMt  +  = 

puAZt  +  ( pu 2  +  p)Mt 

,  h  = 

dA 

~P~d£ 

eAZt  +  u  (e  +  p)AZx_ 

0 

(2.5) 


where  subscripts  x  and  t  imply  partial  differentiation  and  A  —  A/£r  The  system  of  Equation 
(2.4),  together  with  the  ideal-gas  equation  of  state 


P=(7- l)  (e--|u2), 


(2.6) 


where  7  is  the  ratio  of  specific  heats,  constitutes  the  governing  set  of  one-dimensional  Euler 
equations  with  arbitrary  geometry. 


2.2  IMPLICIT  NUMERICAL  SCHEME 


These  equations  are  numerically  applied  to  the  variable-area  shock-tube  problem  using 
the  implicit  ” Delta”  formation  of  Warming  and  Beam."  Implicit  time  differencing  and  central 
spacial  differencing  are  employed  to  evaluate  the  derivative  terms  in  Equation  (2.4)  while  the 
source  term,  li,  is  evaluated  explicitly.  Equation  (2.4)  then  takes  the  following  difference  form. 

<?;+1  -  Qj  +  (4+t1  -  %T)  +  At  (%  +  °(A^A^2)  = 0  (2-7) 

where  the  subscript  ”j”  refers  to  a  specific  grid  point  and  the  superscript  ”n”  to  a  time  step. 
The  terms  containing  the  vector  $  are  non-linear  functions  of  the  conserved  flow  variables  Q, 
and  the  system  of  Equation  (2.7)  can  be  solved  by  iteration.  However,  as  suggested  by  Lomax,7 
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Equation  (2.7)  can  be  solved  directly  with  only  one  inversion  when  these  terms  are  locally 
linearized.  This  can  be  shown  to  be  equivalent  to  one  iterative  step^  while  the  formal  accuracy 
of  Equation  (2.7)  is  retained.  Performing  the  local  linearization  and  defining  the  increment  in 
the  variable  Q  by 

&Q  =  3;+1  -  Qj ,  (2.8) 

one  obtains  the  ”  Delta1’  form  of  the  algorithm  in  matrix  notation: 

[l+  Ar5eA]"-(A^);  =  -Ar5f(4»-Ar(%,  (2.9) 

where  the  bold  character  indicates  a  k  x  k  matrix  and  8^  indicates  central  spacial  differencing. 
The  above  notation  is  favored  by  Beam  and  Warming0*^  and  clearly  points  to  the  solution  in 
terms  of  the  flow  variable  increment  A(J.  The  solution  of  the  dependent  variables  at  the  next 
time  step,  therefore,  is 


q;+1  -  q;  +  a 


(2.10) 


In  Equation  (2.9),  I  is  the  identity  matrix  and  A  is  the  Jacobian  of  the  convective  terms 
with  respect  to  the  flow  variables,  dS/d$.  The  notation  in  Equation  (2.9)  requires  the  dot 
product  to  be  carried  out  prior  to  the  spacial  differencing.  Although  the  ”  Delta”  formation 
leads  to  numerical  efficiency  and  analytical  simplicity,  it  should  be  pointed  out  that  Briley  and 
McDonald^  were  first  in  extensively  applying  the  local-linearization  concept  and  implementing 
the  algorithm  in  its  ” Non-Delta”  form.  In  the  actual  implementation  of  Equation  (2.9),  the 
left-hand  term  in  brackets  is  a  tridiagonal  system  which  is  easily  solved. 

2.3  EXPLICIT  NUMERICAL  SCHEME 

The  Beam-and-Warming  Algorithm  is  generally  preferred  as  the  robust  solution  technique 
for  stiff  partial-differential  equations.  For  the  unsteady  problems  treated  in  the  BRL-Q1D  code, 
however,  it  was  found  that  the  MacCormack  Explicit  scheme  captures  the  shock  within  the 
least  number  of  grid  points  and  consumes  the  lesser  amount  of  machine  time.  For  this  reason, 
it  is  presented  here  as  the  preferred  method  where  resolution  of  surface  discontinuities  is 
required. 

The  MacCormack  explicit  scheme  is  a  second-order,  non-centered,  predictor-corrector 
scheme5  that  alternatively  uses  forward  and  backward  differences  for  the  two  steps  as  follows. 
In  the  predictor  step,  Q  is  defined  by  . 

Qf*  =  Q;  -  Ar  (Af  fi;  )  .  (2.11) 

In  the  corrector  step,  the  new  Q  is  defined  by 

Q"+1  =  Q,*1  +  Qj  -  Ar  (vf  S;+1  )  +  S;  .  (2.12) 


In  both  equations,  the  barred  superscripts  refer  to  predicted  values;  in  particular,  Z?;n+1 
implies  that  the  flux  vector  is  evaluated  using  elements  of  the  predicted  Q-vector.  The  symbols 
A  and  v  are  used  as  the  standard  forward  and  backward  difference  operators,  respectively. 
The  quantity  D  represents  a  fourth-order  dissipation  term,  the  effect  of  which  is  governed  by  an 
empirical  constant  e.  The  dissipation  term  is  described  in  the  following  section. 
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2.4  STABILITY  CRITERION  AND  DISSIPATION 


The  typical  stability  limitation  for  the  system  of  Equation  (2.9)  is  the  Courant, 
Friedrichs  and  Lewy  (CFL)  stability  criterion.5  Application  of  this  criterion  prevents  any  small 
disturbance  (i.e.  sound  wave)  from  traveling  farther  than  one  cell  width  in  the  grid  in  one  time 
increment. 


CN  =  (|«|  +  a)  4“  -  L 


(2.13) 


The  CFL-criterion  is  used  in  the  BRL-Q1D  code  to  determine  the  time  step  At  from  the 
smallest  cell  width  in  the  grid  Ax  using  an  empirically  determined  value  for  the  Courant 
number  CN . 


To  control  phase  errors  associated  with  the  highest  frequencies,  a  fourth-order  dissipation 
term  is  explicitly  added  to  Equation  (2.1).  This  term  is  of  the  following  form. 


Q1+2  ~  ^Qj+i  +  dQj  -  4 4-  <2,-2 

(AO4 


=  (AeVe)2  Qj , 


(2.14) 


where  A  and  V  represent  forward-  and  backward-difference  operators  respectively.  The  final 
form  of  the  computational  algorithm  which  has  been  programmed  is 

[  I  +  Ar5f  A]  "  •  (A$)j  =  -  A r$(  (Stf  -  A r  (%  +  e  (AfVf)2  ■  (2.15) 

Since  the  order  of  the  dissipation  term  is  higher  than  the  order  of  the  truncation  error  of  the 
model  difference  equation,  formal  accuracy  is  maintained.  At  boundary  grid  points  the 
dissipation  is  of  second  order. 


2.5  INITIAL  AND  BOUNDARY  CONDITIONS 


The  variables  in  the  governing  equations  were  non-dimensionalized  by  the  following 
relationships. 


f  =  ral/L  ,  p  —  p/pA ,  p~p/pAa f,  £  —  £/L  ,  u  =  u/ ax  ,  e  =  ej pAa\  .  (2.16) 


The  subscripts  1  and  4  refer  to  the  initial  atmospheric  and  driver  conditions,  respectively.  L  is 
the  reference  length  of  the  shock  tube,  and  a  is  the  sound  speed.  The  computation  is  then 
initialized  with  the  following  conditions, 


.  r4i 

64  7  (7-1) 


(2.17) 


Pi  =  Ti 


7  Pu 


L  41 


7 (7-1) 


(2.18) 


where  Tu  and  P41  are  the  initial  temperature  and  pressure  ratio  across  the  diaphragm  of  the 
shock  tube. 


At  the  endpoints  of  the  grid,  boundary  conditions  were  defined.  The  reflective  boundary 
at  the  left-hand  side  of  the  grid  (i.  e.  the  closed  end  of  the  driver)  was  computationally  modeled 
by  means  of  image  points,  such  that  pj  =  r/jo3,  «j  =  -u3,  u2  =  0,  and  e1  —  e3.  At  the  right-hand 
side  of  the  grid  (at  the  open  end  of  the  expansion  tube)  boundary  conditions  for  outflow  as  well 
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as  for  inflow  were  considered.  For  outflow,  the  static  pressure  is  specified;  for  inflow,  the  static 
pressure  and  the  density  are  specified.  The  remaining  flow  variables  are  then  computed  from 
one-sided  differences  at  the  exit  plane  using  backward  differencing  in  space  and  forward 
differencing  in  time. 


3.  BRL-Q1D  CODE  DESCRIPTION 


The  BRL-Q1D  code  was  developed  to  study  the  flow  in  complex  shock  tubes  with 
arbitrary  area  changes.  This  development  was  done  in  view  of  simulating  in  a  simple  manner 
the  multiple-driver  CEG-LBS^  in  France  and  subsequently  to  study  other  LBS  designs.  This 
chapter  discusses  the  input  options  and  the  basic  structure  of  the  code.  The  user’s  input 
instructions  are  given  in  the  Appendix. 

3.1  SHOCK-TUBE  CONFIGURATIONS 

A  variety  of  shock-tube  geometries  can  be  accommodated  by  the  code  by  combining  and 
altering  the  basic  components  of  the  shock  tube:  Driver,  diaphragm  section,  and  expansion 
tube.  Figure  1  illustrates  some  of  the  possible  configurations. 

3.1.1  Drivers 

The  driver  may  be  shaped  as  a  cylinder,  a  series  of  up  to  four  cylinders,  each  with 
different  diameter,  or  as  a  frustrum  of  a  cone.  The  driver  diameter  is  independent  of  the 
diameter  of  the  expansion  tube  and  may  be  larger,  equal,  or  smaller.  The  cylindrical  driver  may 
be  equipped  with  a  series  of  baffles  each  of  which  is  modeled  as  a  parabolic  constriction.  The 
driver  length  may  be  chosen  at  the  discretion  of  the  user. 

3.1.2  Diaphragm  Section 

The  diaphragm  section  may  be  just  the  locus  of  the  diaphragm,  or  a 
convergent/divergent  nozzle  with  the  diaphragm  located  in  the  throat  area.  A  throttled 
diaphragm  opening  may  be  modeled  with  a  parabolic  area  constriction  about  the  diaphragm 
location.  Instead  of  an  instantaneously  opening  diaphragm,  or  in  combination  with  it,  the  user 
may  define  a  fast-opening  and/or  closing  blast  valve. 

3.1.3  Expansion  Tube 

The  expansion  tube  (or,  driven  tube)  may  be  of  any  length,  but  is  assumed  to  have  a 
constant  diameter  which  is  used  as  a  reference  value.  A  blockage  in  the  test-section  area  may 
be  modeled  as  a  parabolic  area  constriction.  At  the  open  end  of  the  expansion  tube,  a 
rarefaction-wave  eliminator  (RWE)  may  be  modeled  as  a  linear  area  constriction  which  may  be 
varied  linearly  with  time. 

3.2  GRIDDING  TECHNIQUES 

Three  gridding  methods  are  presently  incorporated  in  the  BRL-Q1D  code;  they  are  linear, 
multi-linear  and  non-linear  gridding  as  illustrated  in  Figure  1. 
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(a)  Conventional  Shock  Tube  with  25%  Blockage  in  Test  Section 


(b)  Shock  Tube  with  Valve  (modeled  as  Parabolic  Constriction)  -  Linear  Grid 


(c)  Shock  Tube  with  Baffled  Driver  and  Converging/Diverging  Nozzle 


(d)  Baffled  Driver  Detail  -  Multi-Linear  Grid 


[ 


\ 

y 


(e)  Shock  Tube  with  Stepped  Driver,  Nozzle  and  RWE 


(f)  Shock  Tube  with  Conical  Driver  and  Nozzle  -  Non-Linear  Grid 


Figure  1.  Quasi-One-Dimensional  Geometries  and  Numerical  Grids 
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3.2.1  Linear  Grid 

The  simplest  method  is  a  grid  of  constant  cell  width  (Figure  lb).  It  is  established  by 
dividing  the  number  of  grid  points  into  the  total  length  of  the  shock  tube.  This  method  is 
sufficient  for  geometries  with  moderate  and  gradual  changes  of  the  cross-sectional  area.  A  small 
number  of  grid  points  may  be  chosen  for  a  "quick  and  dirty”  calculation.  But  for  an  accurate, 
converging  solution,  a  high  number  of  grid  points  (>600)  is  required. 

3.2.2  Multi-Linear  Grid 


The  second  method  allows  clustering  the  grid  points  in  critical  areas  by  superimposing  a 
fine  grid  over  the  basic,  coarse  grid  (Figure  Id).  This  method  is  to  be  used  when  there  are  more 
than  one  critical  areas  in  the  shock  tube,  e.  g.  when  there  are  baffles  in  the  driver,  or  when  the 
driver  length  is  short  compared  to  the  total  length,  and  an  RWE  is  defined  at  the  other  end  of 
the  shock  tube.  This  method  saves  grid  points  but  requires  a  greater  number  of  time  steps  than 
the  first  method  because  the  time  step  depends  on  the  smallest  cell  size  in  the  grid  which  has 
been  generated  by  the  clustering. 

3.2.3  Non-Linear  Grid 

The  third  method  allows  clustering  about  a  single  critical  area  as  e.g.  the 
convergent/divergent  nozzle  (Figure  If).  Typically,  the  diaphragm  location  would  be  chosen  as 
the  cluster  point.  This  method  is  particularly  convenient  when  the  initial  pressure  ratio  across 
the  diaphragm  is  high  and  the  shock  has  to  be  stretched  over  a  few  cells  to  avoid  instabilities 
(i.e.  negative  density)  in  the  computation.  But  the  user  must  make  sure  that  there  are  no  other 
critical  areas  far  away  from  the  clustering  location  which  could  be  denuded  of  a  sufficient 
number  of  grid  points.  To  help  the  user  with  this  judgement,  a  graph  of  maximum  and 
minimum  cell  size  versus  the  clustering  parameters  was  developed. 


3.2.4  Choosing  the  Clustering  Parameters. 

The  clustering  function  used  for  the  j-th  grid  point  is 

+  Xc!ust  (3.1) 


X,  —  —  sinh 

J  /? 


In  G>  +  J-  2V  {In  Cx  -  In  C2) 

J  max  —  * 


where  Cx  and  C2  are  the  following  functions, 


C\  =  P  (1  -  Ac/usf)  +  \J ft2  (1  -  Xc/USf) 2  +  1 


and 


C2 - 0Xclust  +  y/W^dust  +  I 


(3.2) 


(3.3) 


Since  the  clustering  function  is  already  programmed  into  the  code,  it  was  easy  to  execute 
a  series  of  calculations  with  a  single  time  step  to  find  the  value  of  the  maximum  cell  size  at  the 
location  farthest  away  from  the  cluster  point,  as  well  as  the  value  of  the  minimum  cell  size  at 
the  cluster  point  as  a  function  of  the  two  input  constants  and  Xdiist .  These  values  were  then 
normalized  by  the  average  cell  size, 


Az  =  Lref  /  (w  -  2) 


(3.4) 


and  plotted  versus  ft  for  constant  Xchst  (Figure  2). 
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Figure  2.  Grid-Cell  Size  as  a  Function  of  /?  and  Xcitlst  . 
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For  given  f)  and  Xciust ,  the  graph  will  yield  the  smallest  and  the  largest  cell  size  in  the 
grid.  Conversely,  for  a  known  cluster  location  and  a  predetermined  critical  cell  size,  the  graph 
will  yield  the  required  /3-factor.  The  critical  cell  size  may  be  chosen  at  the  clustering  point 
(minimum),  or  at  the  point  in  the  grid  farthest  away  from  it  (maximum). 

3.3  CODE  STRUCTURE 

The  BRL-Q1D  code  is  an  algorithm  for  implicitly  or  explicitly  solving  the  quasi-one- 
dimensional  (Q1D)  Euler  equations  with  arbitrary  geometry.  The  code  was  written  in 
FORTRAN  77  to  be  executed  on  the  CYBER  system  at  the  BRL.  Another  version  of  the  code 
was  adapted  for  execution  on  the  Vax  11/780. 

3.3.1  The  Main  Program 

The  main  program  has  two  functions:  It  instructs  the  user  on  how  to  use  the  code,  and 
it  coordinates  the  information  flow  by  calling  the  appropriate  subprograms.  The  user 
instructions  are  contained  in  a  comment  section  which  is  located  at  the  head  of  the  main 
program.  They  are  presented  in  the  Appendix  (page  55-61).  Secondly,  the  main  program  serves 
as  a  calling  program  controlling  the  information  flow  which  is  charted  in  Figure  3. 

At  the  start,  the  program  reads  the  necessary  input  data  which  are  contained  on  from  six 
to  twenty-one  card  images  depending  on  the  complexity  of  the  shock-tube  configuration. 
Several  cases  of  flow  computations  may  be  stacked  ia  sequence;  the  program  will  continue  to 
read  a  new  set  of  input  data  every  time  a  computation  has  been  finished.  After  the  input  data 
have  been  read  they  are  normalized  and  reflected  back  to  file  OUTPUT  (see  pages  64-66  in  the 
Appendix).  Next,  a  finite-difference  grid  is  layed-out  according  to  the  specified  grid  parameters. 
Then  the  cross-sectional  areas  and  the  initial  flow  conditions  are  defined  at  all  grid  point.  The 
user  then  has  the  option  of  writing  these  initial  conditions  to  file  OUTPUT  and/or  to  file 
TAPE15  for  later  plotting. 

The  time  loop  is  entered  and  a  new  time  step  (At)  is  calculated.  Now  the  flow  variables 
are  calculated  for  the  new  time  at  all  interior  points  of  the  grid  using  either  the  implicit  Beam- 
Warming,  or  the  explicit  MacCormack  scheme.  Then  the  flow  variables  are  defined  at  the 
boundaries  and  the  flow  variables  may  be  written  to  files  OUTPUT  and/or  TAPE15  as  chosen 
by  the  user.  Finally,  the  time-dependent  areas  in  the  valve,  or  the  RWE  are  redefined  and 
control  then  returns  to  the  beginning  of  the  time  loop. 

There  are  two  checks  incorporated  in  the  time  loop  in  order  to  prevent  involuntary 
abortion  of  the  computation  by  the  central  processor.  The  first  is  a  check  on  the  size  of  the 
time  step.  We  observed  in  some  computations  that  the  time  step  would  go  to  zero  because  of 
excessive  local  flow  velocities.  When  this  happened,  an  overflow  condition  occurred  which 
caused  the  job  to  abort.  Therefore,  a  conditional  loop  exit  based  on  the  size  of  the  time  step 
was  introduced  which  prevents  this  overflow  condition  from  occurring  and  allows  the 
computation  to  conclude  normally.  The  second  check  is  made  on  the  flow  density.  A  negative 
density  is  physically  unreal  and  may  be  caused  by  too  large  a  time  step.  Once  a  negative 
density  occurs,  the  time  loop  is  exited  and  the  computation  concludes  normally. 

Upon  conclusion  of  the  time  loop,  a  set  of  blast-wave  characteristics  for  the  pressure 
history  is  determined  which  characterize  the  simulated  blast  wave  relative  to  free-field  blast 
waves.  The  pressure,  density  and  velocity  histories  (versus  time)  are  written  to  files  OUTPUT 
and/or  TAPE16  together  with  these  blast-wave  characteristics.  Control  then  returns  to  the 
beginning  of  the  main  program  where  an  attempt  is  made  to  read  the  next  set  of  input  data.  If 
an  end-of-record  mark  is  encountered  instead,  the  program  stops  executing. 
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Figure  3.  Basic  Information-Flow  Structure  in  the  BRL-Q1D  Code. 
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3.3.2  The  Subprograms 


There  are  presently  25  user-defined  subprograms  in  the  BRL-Q1D  code,  each  of  which 
performs  one  particular  aspect  of  the  computation.  All  but  one  subprogram  are  subroutines. 
The  first  eight  subroutines  are  related  to  the  input;  they  set  up  the  computation.  The  next  nine 
subroutines  constitute  the  QlD-algorithm  proper;  they  carry  out  the  computation.  The  last 
eight  subprograms  prepare  and  store  the  output  data.  The  names  and  purposes  as  well  as  the 
calling  routines  of  all  these  subprograms  are  as  follows. 


READIN 

This  SUBROUTINE  reads  the  input  data  from  card  images 
stores  them  in  common  and  echoes  them  to  output. 

MAIN 

NORMAL 

This  SUBROUTINE  normalizes  the  physical  input  data, 
areas:  Ar  — ►  Ar/L%ef 

lengths:  X  -►  X/Lre/ 

diameters:  D  DjLrSf 

time:  t  -+  taanb/Lref 

MAIN 

MULTLIN 

This  SUBROUTINE  defines  a  multi-linear  grid-point 
distribution  in  the  computational  scheme,  and  transforms 
the  x-coordinate  from  the  physical  to  the  computational 
(£-)  plane,  when  MTRX  =  2. 

MAIN 

METRIX 

This  SUBROUTINE  defines  the  grid-point  distribution  in 
the  computational  scheme,  and  transforms  the  x-coordinate 
from  the  physical  to  the  computational  (£-)  plane 
when  MTRX  7^  2. 

MAIN 

AREA 

This  SUBROUTINE  determines  the  cross-sectional  areas 
at  all  grid  points. 

MAIN 

ARVA 

This  SUBROUTINE  defines  the  areas  in  the  blast  valve  as 
a  function  of  time. 

MAIN 

ARWE 

This  SUBROUTINE  defines  the  areas  in  the  rarefaction 
-wave  eliminator  (RWE). 

MAIN 

DRIVOL 

This  subroutine  computes  the  reference  area  and  the  driver 
volume. 

MAIN 

INITIA 

This  SUBROUTINE  sets  the  initial  conditions  at  each  grid 
point.  The  dependent  variables  of  Euler’s  conservation 
equations  are  stored  in  the  Q-array: 

Q1  =  p-Ar  (mass) 

Q2  =  pu'Ar  (momentum) 

Q3  =  e-Ar  (energy) 

MAIN 

EIGEN 

This  SUBROUTINE  finds  the  largest  velocity  gradient  in 
the  grid  and  determines  the  time  step,  At,  from  the  CFL 
stability  criterion  for  compressible  flow. 

MAIN 
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EULERI 

This  SUBROUTINE  solves  the  Euler  equations  using  the 
Beam-Warming  implicit  scheme  which  involves  the  inversion 
of  the  block  tridiagonal  matrix  of  the  flow  equations.  This 
subprogram  calls  subroutines  SMOOTH,  FILTRX,  BTRI 
and  BC. 

MAIN 

MACCOR 

This  SUBROUTINE  solves  the  Euler  equations  using  the 
explicit  MacCormack  technique.  This  subroutine  calls 
subprograms  BCS,  SMOOTH  AND  BC. 

MAIN 

SMOOTH 

This  SUBROUTINE  does  fourth  order  smoothing  for 
interior  grid  points,  2  <  j  <  w«-l.  and  second  order 
smoothing  at  the  grid  points  j  =  2  and  Jmm-l. 

EULERI 

MACCOR 

FILTRX 

This  SUBROUTINE  fills  the  block-tridiagonal  A-matrix 
for  the  implicit  solution. 

EULERI 

BTRI 

This  SUBROUTINE  inverts  the  block  tridiagonal  matrix 
(Beam-Warming  solution).  Each  block  is  a  3x3  element 
matrix. 

EULERI 

BCS 

This  SUBROUTINE  determines  the  flow-field  parameters 
at  the  right  (=open)  and  left  (closed)  boundaries  for 
MacCormack’s  predictor  step. 

MACCOR 

BC 

This  SUBROUTINE  defines  the  flow-field  parameters  at 
the  right,  open-end  (with  RWE)  and  left,  closed-end 
boundaries  of  the  shock  tube  at  the  n-th  time  step. 

EULERI 

MACCOR 

WRIT 

This  SUBROUTINE  writes  the  geometrical  shock-tube 
parameters  and  the  gas  parameters  to  file  OUTPUT. 

MAIN 

BC 

WRITXP 

This  SUBROUTINE  prints  arrays,  10  numbers  per  line, 
mixed  with  lines  of  output  denoting  constant-value  strings. 

WRIT 

PLOTX 

This  SUBROUTINE  writes  gas  parameters  versus  X  on 
file  TAPE  15  for  later  plotting. 

MAIN 

OUTSTA 

This  SUBROUTINE  converts  gas  parameters  at  test  station 
XSTA(J)  to  dimensional  values  and  stores  them  in  arrays 
for  later  plotting  and/or  printing. 

MAIN 

IMPULS 

This  SUBROUTINE  finds  the  time  of  shock  arrival,  the 
shock  overpressure  and  the  positive-phase  duration  at  test 
station  XSTA(J).  It  determines  the  static  and  dynamic 
pressure  impulses  using  Simpson’s  rule  of  integration  and 
calculates  the  equivalent  nuclear  weapon  yield  for  both 
impulses.  This  subroutine  calls  subprogram  DVDINT. 

MAIN 

BASER4 

This  BLOCK  DATA  subprogram  specifies  the  data  base 
for  a  selected  nuclear  blast-wave  reference. 

IMPULS 
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DVD  I  NT 

This  SUBROUTINE  does  divided-difference  interpolation. 

The  arguments  have  been  tabulated  in  descending  order. 

IMPULS 

OUTPOT 

This  SUBROUTINE  writes  the  gas  parameters  which  were 
determined  at  the  test  station  by  subprogram  OUTSTA, 
to  files  OUTPUT  and  TAPE16  for  later  plotting. 

MAIN 

3.3.3  Code  Execution 

The  BRL-Q1D  code  is  typically  executed  on  the  CYBER  7600  at  BRL  which  runs  under 
the  SCOPE  2  operating  system  of  CDC.  It  requires  67,303  words  octal  in  Small-Core  Memory 
(SCM)  and  575,550  words  octal  in  Large-Core  Memory  (LCM).  The  CPU-time  required  for 
execution  is  less  than  1/8  ms  per  grid  point  per  time  step  when  the  code  is  compiled  with 
optimization  level  2.  Two  sample  job  streams  are  shown  in  the  Appendix  (page  62-63). 

The  output  is  stored  on  two  tape  files  (TAPE15  and  TAPE16)  which  may  be  stored  on 
disc  for  later  reference  and  serve  as  input  to  two  plotting  codes  which  are  executed  separate 
from,  and  subsequent  to,  the  BRL-Q1D  code.  This  separation  was  chosen  to  limit  the  SCM  field 
length  required  for  execution.  Although  a  request  for  greater  field  length  can  be  included  on  the 
job  card,  it  is  undesireable  because  it  delays  the  job  execution  until  the  requested  memory  is 
made  available  by  a  human  operator. 

3.4  OUTPUT  OPTIONS 

The  results  of  the  QlD  computation  are  printed  in  tabulated  form  or  stored  on  disc  for 
subsequent  plotting.  Sample  tabulations  and  plots  are  presented  in  the  Appendix. 

3.4.1  Tabulated  Output 

There  are  two  types  of  printed  output.  For  data  analysis,  the  complete  arrays  of  density, 
velocity  and  pressure  are  printed  out  at  user-defined  time  steps  together  with  the  arrays  of 
grid-point  location  (X),  d^jdx  (XIX),  areas  (AR)  and  its  derivatives,  dA/dx  (ARX)  and  A  (ART). 
See  the  example  in  the  Appendix  (page  67-74).  For  plotting,  these  arrays  are  written  to  disk  via 
local  file  TAPE15  at  user-defined  time  steps. 

The  pressure  history  is  printed  in  tabulated  form  at  the  end  of  the  computation  for  user- 
defined  time  intervals.  It  contains  a  heading  with  useful  summary  data  like  shock-arrival  time, 
peak  static  overpressure,  positive-phase  duration  and  static-overpressure  impulse.  In  the  table, 
the  static  (side-on)  overpressure  and  impulse,  the  dynamic  pressure  and  impulse,  the  flow 
velocity  and  the  density  are  listed  versus  time.  A  sample  listing  is  shown  in  the  Appendix  (page 
75-76).  The  pressure  history  is  also  stored  on  disc  for  later  plotting  if  the  user  so  chooses.  The 
local  file  for  this  record  is  TAPE16. 

3.4.2  Plotted  Output 

Two  plotting  codes  are  available,  PX1DPLT  for  plotting  pressure,  density  or  velocity 
versus  distance,  and  PT1DPLT  for  plotting  pressure  and/or  density  versus  time.  Both  plot 

codes  are  written  in  FORTRAN  77  using  the  DISSPLA-Version  8.2  plotting  package. ^  Both 
codes  contain  a  comment  section  with  the  necessary  user  instructions.  They  may  be  obtained 
from  the  authors  upon  request. 
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The  density-versus-distance  plot  is  very  useful  for  studying  the  wave  propagation  in  the 
shock-tube  flow.  The  incident  shock,  the  contact  surface,  the  recompression  shock  and  various 
shock  reflections  can  be  recognized  easily  in  the  sample  plot  shown  in  the  Appendix  (page  77). 
The  data  given  in  the  heading  serve  the  identification  of  the  plot.  Below  the  plot  title  are  listed 
from  1  to  6  test  stations  (XSTAj),  the  driver  volume  and  pressure  and  the  reference  length  of 

the  shock  tube.  In  the  second,  right- adjusted  column  are  listed  the  case  identifier  and  the  plot 
number,  the  offset  value,  Ay,  for  multiple-curve  plots,  and  comments  on  shock  stretching  and 
grid  clustering  when  applicable.  A  sketch  of  the  shock-tube  configuration  is  drawn  below  these 
data  and  above  the  plot  with  or  without  grid  indication. 

A  sample  plot  of  a  pressure  history  is  shown  in  the  Appendix  (page  78).  The  heading 
data  specify  the  shock-tube  dimensions  in  column  1,  the  test  conditions  in  column  2  and  the 
blast-wave  characteristics  in  column  3.  The  shock-tube  dimensions  listed  under  the  heading 
include  the  case  identification,  the  total  (i.e.,  Reference)  length  of  the  facility  (L-ref),  the  driver 
length  (L-drv)  and  volume  (V-drv),  the  length  of  the  expansion  tube  (L-dvn)  and  the  length  of 
the  RWE  (L-rwe). 

The  test  station  and  initial  conditions  are  listed  in  column  2.  The  pressure  history  is 
recorded  at  the  test  section,  X-sta,  the  location  of  which  is  given  in  metres  from  the  exit  of  the 
divergent  nozzle.  The  initial  conditions  include  the  driver  pressure  (P-drv),  the  ambient, 
atmospheric  pressure  (P-amb),  the  ambient  Temperature  (T-amb)  and  the  temperature  ratio 
across  the  diaphragm  (T^/T^). 

The  third  column  lists  the  blast-wave  characteristics,  shock  overpressure  (P-so),  shock- 
arrival  time  (t-a),  positive-phase  duration  (PPD),  the  static-overpressure  impulse  (I-so),  the  peak 
dynamic  pressure  (Q-s)  and  the  dynamic-pressure  impulse  (I-dyn).  The  values  in  parentheses 
are  the  related  free-field  data  for  PPD  and  weapon  yield.  They  were  obtained  by  interpolation 
from  the  reference  data  base  as  functions  of  the  shock  overpressure. 


4.  RESULTS  AND  DISCUSSION 


This  chapter  presents  the  results  of  several  computations  of  the  flow  in  conventional 
shock  tubes  and  in  shock  tubes  with  varying  crossectional  area. 

4.1  SIMULATION  OF  FLOW  IN  CONVENTIONAL  SHOCK  TUBES 

Figure  4  shows  a  comparison  of  four  computations  with  experimental  data  from  the 

BRL-24”  shock  tube.^'^  The  test  conditions  are  listed  in  Table  I.  Also  listed  are  three  blast- 
wave  characteristics  which  were  gained  from  the  pressure  curves  in  Figure  4.  These  are  the 
shock  overpressure,  p$oy  the  positive-phase  duration,  t+y  and  the  static-overpressure  impulse,  Iso. 

From  the  comparison  of  the  computational  and  experimental  values  of  these  blast-wave 
characteristics  we  learn  that  the  computation  predicts  a  higher  shock  overpressure  than  the 
experimental  value.  At  the  same  time,  the  computational  driver  empties  faster  leading  to 
shorter  positive-phase  durations.  This  is  to  be  expected  from  the  computation  because  the 
inviscid  Euler  equations  do  not  account  for  flow  losses. 

The  positive  area  under  the  static-overpressure  curve  in  Figure  4  was  then  integrated  to 
find  the  static-overpressure  impulse.  In  the  code,  this  was  done  numerically,  using  the 
trapezoidal  rule.  For  the  experimental  records,  it  was  done  graphically.  The  computational  and 
experimental  values  are  close,  although  differences  should  be  expected  because  of  the  different, 
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(a)  Driver  Pressure  1.87  atm  (pgo  =  35.9  kPa) 


(b)  Driver  Pressure  3.15  atm  (pgQ  =  72.2  kPa) 


Figure  4.  Computational  and  Experimental  Pressure  Histories  for  the  BRL-24”  Shock  Tube. 
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(c)  Driver  Pressure  4.69  atm  (pgQ  =  101  kPa) 


(d)  Driver  Pressure  6.91  atm  (pgo  =  142  kPa) 
Figure  4.  (concluded) 
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TABLE  I. 

TEST  CONDITIONS  AND  RESULTS  FOR  BRL-24”  SHOCK-TUBE  EXPERIMENT 

Dimensions 

-v 

Driver  Length 

-* 

10.85  metre 

Inner  Diameter  of  Driver  Tube 

0.575  metre 

Center  of  Square  Test  Section 

41.33  metre 

Height  and  Width  of  Test  Section 

0.508  metre 

Length  of  Expansion  Tube 

76.20  metre 

Total  Length  of  Shock  Tube 

87.05  metre 

Test  Conditions 

Test  Nr. 

1 

2 

3 

4 

p 

a  mb 

kPa 

102  6 

102.6 

102.6 

102.7 

T 

a  mb 

K 

293.7 

293.7 

293.7 

293.8 

p4/pi 

— 

1.87 

3.15 

4.69 

6.91 

VT1 

— 

1 

1 

1 

1 

Results 

P 

so,  exp 

kPa 

35.9 

71.5 

99.0 

140. 

P 

so, comp 

kPa 

36.7 

75.0 

110. 

148. 

t 

+,  exp 

s 

.103 

.131 

.203 

.213 

t 

+,comp 

s 

.090 

.118 

.142 

.156 

I 

so,  exp 

kPa-s 

2.50 

5.57 

9.43 

13.91 

I 

so, comp 

kPa-s 

2.49 

5.64 

9.11 

13.12 

i.e.  inviscid  versus  real,  flows.  As  it  appears,  the  higher  shock  overpressures  compensate 
somewhat  for  the  shorter  positive-phase  durations  in  the  computation. 

Figure  5  shows  a  comparison  of  the  explicit  MacCormack  scheme  with  the  implicit 

Beam-Warming  method  for  a  flow  simulation  in  the  BRL-2.44m  shock  tube.^  The  smoothness 
of  the  ”  implicit”  pressure  curve  bears  witness  to  the  stability  of  the  Beam-Warming  scheme. 
The  MacCormack  scheme  proved  to  be  much  more  sensitive  to  area  changes  and  prone  to 
develop  instabilities  early  in  the  computation. 
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FIGURE  5:  Comparison  of  the  Explicit  and  Implicit  Methods 


4.2  SIMULATION  OF  THREE  CEG-LBS  TESTS 


This  section  shows  how  a  geometrically  complicated  blast-wave  simulator  can  be  modeled 
in  a  quasi-one-dimensional  sense,  and  how  the  computational  results  from  the  BRL-Q1D  code 
computations  compare  with  the  experimental  static-,  and  dynamic-pressure  measurements  from 
the  CEG-LBS  facility. 

4.2.1  Blast-Wave  Simulator  Model 

The  CEG-LBS^  consists  of  seven  cylindrical  driver  tubes  of  various  lengths  which  feed 
compressed  air  through  nozzles  into  a  reinforced-concrete  expansion  tube  (Figure  6a).  The 
number  of  driver  tubes  was  chosen  such  as  to  limit  the  diameter  of  the  pressure  vessels  to 
certain  material-strength  and  manufacturing  criteria.^  The  convergent-divergent  nozzles  were 
chosen  to  retard  the  emptying  of  the  drivers  and  to  keep  the  diaphragms  at  a  managable  size. 

At  the  open  end  of  the  expansion  tube  sits  a  rarefaction-wave  eliminator  (RWE)  which 
consists  basically  of  an  array  of  louvers  (shown  removed  from  the  exit  of  the  expansion  tube  in 
Figure  6a).  The  purpose  of  the  active  RWE  is  to  regulate  the  outflow  of  air  from  the  expansion 
tube  as  a  function  of  time  by  closing  and  opening  the  louvers.  This  is  done  in  such  a  manner 
that  the  compression  and  expansion  waves  from  the  open  end  combine  within  the  tube  so  as  to 
generate  the  effect  of  having  a  much  longer  expansion  lube  behind  the  test  section. 

The  combination  of  variable  length  driver  tubes  and  active  RWE  produces  a  decaying 
blast  wave  in  the  expansion  tube  of  the  CEG  facility.  The  diaphragms  for  all  the  drivers  are 
located  in  the  throat  section  of  the  nozzles.  When  they  burst  simultaneously,  the  pressure 
discontinuity  across  the  diaphragms  generates  left-travelling  rarefaction  waves  which  reflect 
from  the  closed  ends  of  the  driver  tubes  at  different  times.  After  their  reflection  they  travel 
forward  and  weaken  the  incident  shock  progressively  in  the  early  part  of  the  pressure  history. 
At  later  times,  the  action  of  the  RWE  governs  the  shape  of  the  blast  wave. 

The  quasi-one-dimensional  computational  model  for  this  facility  is  shown  in  Figure  6b. 
The  seven  drivers  of  the  CEG-LBS  are  combined  into  one  by  lumping  the  cross-sectional  areas 
at  any  given  locations  into  a  single  area  and  the  stair-stepped  driver  results.  The  four  steps 
indicate  that  the  French  facility  uses  four  different  driver  lengths.  The  RWE  is  modeled 
computationally  as  a  single  lumped  opening  which  may  vary  with  time. 

Because  the  lumping  of  the  areas  changes  the  cone  angles  of  the  convergent  and 
divergent  nozzle  sections,  the  suggestion  was  made  to  use  a  1/7  lumped- area  model  in  the 
computations.  In  this  model,  the  single  nozzle  remains  unchanged,  the  expansion  tube  has  1/7 
of  the  cross-sectional  area  of  the  CEG-LBS,  and  the  stepped  driver  has  1/7  of  the  7-driver 
volume  of  the  French  facility.  This  stepped  1/7  driver  was  shaped  into  a  frustrum  of  a  cone  by 
determining  the  smaller  end  diameter  from  the  driver  volume  for  the  given  driver  length  and 
the  fixed  entrance  diameter  to  the  nozzle. 

4.2.2  Comparison  with  CEG  Data. 

Several  pressure  histories  from  the  French  blast- wave  simulator  at  CEG^5  are  available 
for  comparison  with  the  computational  results.  From  these,  three  experimental  records  (for 
peak  static  overpressures  of  30,  47  and  67  kPa)  were  chosen  to  be  simulated  computationally. 
Initially,  a  considerable  effort  was  spent  to  simulate  the  active  RWE  of  the  CEG-LBS  and 
results  of  this  work  were  reported  previously.^5  However,  the  number  of  adjustments  to  the 
experimental  RWE  settings  which  was  necessary  to  effect  a  matching  computation  and  our 
inability  to  find  a  correlation  between  the  experimental  and  computational  RWE  settings  has 
convinced  us  that  the  active  RWE  cannot  be  simulated  realistically  by  our  code. 
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(a)  The  CEG-LBS  at  Gramat,  France 


Figure  6.  The  CEG-LBS  and  its  Computational  Equivalent 


Figure  6.  (concluded) 


The  computational  results  shown  here  were  made  for  an  LBS  without  RWE  and  with  an 
extended  expansion  tube  instead,  to  make  up  for  the  effect  of  the  RWE.  The  initial  test 
conditions  and  the  substitute  lengths  are  listed  in  Table  II  and  the  results  are  compared  in 
Figure  7a-c.  In  these  figures  the  noisy  trace  represents  the  experimental  data  whereas  the 
computational  results  are  the  smooth  line.  The  solid  lines  are  for  the  static  overpressure  while 
the  dashed  lines  indicate  the  computed  dynamic  pressure. 


TABLE  II. 

INITIAL  CONDITIONS  FOR  CEG  TESTS 

Test 

p4/p. 

VTi 

^amb 

Tamb 

A  LEf‘ 

Nr. 

— 

— 

kPa 

'C 

m 

5b 

10.1 

.982 

99.0 

12 

200 

9 

15.4 

.993 

98.8 

4 

210 

8 

26.8 

.994 

99.0 

7 

200 

a  substitute  for  active  RWE  in  computation 
target  in  test  section:  25%  blockage 

The  comparisons  show  that  the  general  features  (wave  reflections  and  expansions)  of  the 
pressure  histories  are  replicated  in  the  computational  simulations.  However,  distinct  differences 
in  the  wave  patterns  can  be  observed  also  because,  due  to  the  absence  of  the  RWE  in  the 
computation,  all  the  wave  shaping  is  done  by  the  expansion  waves  coming  out  of  the  driver.  In 
the  first  comparison  (Figure  7a)  the  agreement  between  the  computational  and  experimental 
data  is  good.  The  computational  peak  overpressure  is  only  10%  below  the  experimental  value 
and  the  computational  curve  follows  the  experimental  trace  as  well  as  can  be  expected. 

Further  study  of  the  pressure  histories  in  Figures  7b  and  7c  reveals  however  that  the 
computation  fails  to  replicate  the  peak  overpressure  with  increasing  driver  pressure.  A 
comparison  of  the  pressure  histories  shows  that  the  computational  history  falls  initially  below 
the  experimental  history  the  more  the  higher  the  driver  pressure.  In  the  second  comparison 
(Figure  7b)  the  computational  peak  overpressure  is  25%  below  the  experimental  value  at  a 
driver  pressure  of  15.4  atm  and  in  the  third  comparison  the  peak  overpressure  is  39%  below  the 
experimental  value  at  26.8  atm  driver  pressure. 

More  blast-wave  characteristics  are  listed  for  comparison  in  Table  III.  From  the 
comparison  of  the  computational  and  experimental  data  we  learn,  that  none  of  the  blast-wave 
characteristics  are  truly  matched  by  the  computation,  although  the  computations  come  close  to 
the  experiment.  Most  significantly,  the  computational  impulses  fall  below  the  experimental 
values  and  because  the  peak  overpressure  values  are  lower,  the  weapon  yields  turn  out  relatively 
high.  The  differences  can  be  explained,  however. 

It  is  indeed  surprising  that  a  crude  quasi-one-dimensional  model  can  at  all  approximate 
the  complicated,  three-dimensional  flow  process  in  the  CEG-LBS.  Consider,  for  example,  the 
seven  drivers  of  different  lengths.  As  they  are  emptied  on  bursting  of  the  diaphragms, 
rarefaction  waves  empty  the  tubes  at  different  rates.  Subsequent  compression  and  expansion 
waves  from  the  throats  and  RWE  influence  the  flow  differently  in  every  tube.  The  shock  waves 
emerging  from  each  nozzle  form  a  shock  bubble  and  then  coalesce  forming  a  complicated  array 
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(a)  Test  5:  Driver  Pressure  10.1  atm  (psQ  =  30  kPa) 

Figure  7.  Three  Computations  Compared  with  CEG-LBS  Experiments 
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(b)  Test  9:  Driver  Pressure  15,4  atm  (pgQ  =  47  kPa) 


Figure  7.  (continued) 
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(c)  Test  8:  Driver  Pressure  26  8  atm  (p^o  =  67  kPa) 


Figure  7.  (concluded) 


of  Mach  stems  and  spherical  shocks.  The  flow  at  the  exit  of  the  driver  nozzles  experiences  a 
sudden  change  in  area.  The  global  nature  of  the  flow  is  to  form  a  smooth  transition  between 
the  nozzle  exits  and  the  downstream  flow  field  with  a  ” dead-water”  region  in  the  immediate 
vicinity  of  the  nozzles. 


TABLE  IB. 

COMPUTATIONAL/EXPERIMENTAL  BLAST-WAVE  CHARACTERISTICS  FOR  CEG  TESTS 

Test 

p 

Fso 

t 

+ 

I 

so 

Y 

so 

I 

q 

Y 

q 

Nr. 

kPa 

s 

kPa-s 

kT 

kPa 

kPa-s 

kT 

r  a 

exp 

30.2 

1.33 

11.8 

28 

3.10 

0.79 

21 

o 

comp 

28.1 

1.36 

10.4 

22 

4.01 

0.90 

41 

ft 

exp 

46.8 

1.44 

17.3 

42 

7.38 

0.64' 

- 

y 

comp 

35.5 

1.45 

16.1 

50 

4.33 

1.02 

19 

8 

exp 

67.3 

1.40 

30.5 

131 

14.6 

3.02 

31 

comp 

48.3 

1.42 

25.8 

127 

7.87 

2.57 

75 

a  test  section  blockage  25%. 
questionable  value. 

The  computational  model  lumps  the  seven  drivers  and  their  nozzles  into  one.  The  seven 
6°  cone  half-angle,  divergent  nozzles  of  the  CEG-LBS  are  lumped  into  a  single  16°  cone  half¬ 
angle  divergent  nozzle  with  an  area  discontinuity  where  the  nozzle  opens  into  the  expansion 
tube.  The  computation  cannot  generate  a  ’’dead-water”  region,  as  does  the  real  flow  at  this 
location.  The  emptying  of  this  one  driver  occurs  with  a  different,  single  system  of  rarefaction 
and  recompression  waves.  The  three-dimensional  effects  of  the  interacting  flow  from  the  seven 
drivers  are  absent  from  the  model.  The  flow  in  the  vicinity  of  the  nozzle  exit  is  indeed  not  well 
represented  with  a  one-dimensional  model.  It  is  therefore  desirable  to  compare  the 
computational  data  with  the  experimental  results  from  a  quasi-one-dimensional  LBS  model. 

4.3  SIMULATION  OF  BRL-LBS  MODEL-A  EXPERIMENTS 

It  should  be  expected  that  the  flow  from  a  shock  tube  with  variable  area  resembling  the 
computational  model  more  closely  than  the  CEG-LBS  will  eliminate  some  of  the  ambiguity 
introduced  into  the  comparison  by  the  complexity  of  the  mixing  of  the  seven  driver  flows  and 
associated  phenomena.  Therefore,  an  experimental,  axisymmetric,  1/37-scale  model  of  the 
computational  blast-wave  simulator  was  designed  (Figure  8)  and  built  at  the  BR.L  to  validate 
the  predictions  from  the  BRL-Q1D  code.  The  shock  tube  was  used  to  produce  decaying  blast 
waves  from  23  to  230  kPa  (3.3  to  33.3  psi).  Experiments  were  conducted  with  heated  and 
unheated  driver  gas.  The  data  were  recorded  at  a  test  station  located  7  hydraulic  diameters 

IS  17 

downstream  from  the  exit  of  the  divergent  nozzle.  Selected  results  from  these  experiments  ' 
covering  the  whole  shock-overpressure  range  are  compared  with  the  computational  data  in 
Figure  8a-f.  The  initial  test  conditions  are  listed  in  Table  IV. 

The  experiments  are  listed  in  ascending  order  of  shock  overpressure.  From  this 
arrangement  it  is  obvious  that  the  heated  driver  gas  produces  stronger  shocks  at  lower  driver 
pressures  than  the  unheated  driver  gas  at  even  the  highest  driver  pressure  (e.g.  Test  12). 
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Figure  8.  Experimental,  Axisymmetric  BRL-LBS  Model-A 


Initially,  all  computations  were  executed  with  the  nominal  throat  area  implying  that  the 
rupturing  diaphragm  would  open  fully  without  causing  a  partial  obstruction  of  the  throat  area. 
The  actual  throat  areas  were  obtained  from  the  ruptured  diaphragms  after  the  tests^  and  the 
equivalent  diameters  determined  which  are  given  in  the  table.  It  was  found  that  the 
computational  fit  improved  when  the  calculations  were  executed  with  the  actual  rather  than  the 
nominal  throat  area.  The  comparison  of  the  computational  and  experimental  pressure  histories 
in  Figure  8  shows  satisfactory  agreement. 


TABLE  IV. 

INITIAL  CONDITIONS  FOR  LBS-MODEL-A  TESTS 12 . 

Test 

Nr. 

P 

amb 

kPa 

Tamb 

K 

typi 

VT1 

dthr 

m 

a 

p 

^so,exp 

kPa 

D 

*so,comp 

kPa 

Unheated  Driver  Gas 

20 

105.7 

278.5 

6.75 

1.0 

0.0472b 

22.9 

22.0 

3 

101.8 

313.0 

28.37 

1.0 

0.0447 

81.3 

71.6 

12 

101.5 

296.9 

134.7 

1.0 

0.0396 

175.7 

180.5 

Heated  Driver  Gas 

25 

102.0 

299.8 

92.95 

1.593 

0.0463 

190.3 

192.7 

26 

101.4 

299.8 

107.8 

1.593 

0.0435 

201.6 

197.6 

28 

102.7 

299.8 

131.4 

1.778 

0.0472b 

229.6 

251.5 

a  measured  at  Test  Station  7 

b  nominal  throat  diameter 

4.3.1  Experiments  with  Unheated  Driver  Gas 

Three  examples  of  pressure  histories  for  unheated  driver  gas  were  chosen  from  Reference 
17.  Test  20  (Figure  9a)  represents  the  shock-tube  test  with  the  lowest  pressure  ratio  across  the 
diaphragm,  PjP\  =i  6.75.  The  computational-experimental  match  is  very  good  including  the 
spike  at  the  shock  front  even  though  the  computation  was  executed  for  the  nominal  throat  area. 
This  agrees  with  our  previous  observation  in  Section  4.2  (Figure  7a)  where  we  found  good 
agreement  between  computation  and  experiment  at  the  low  driver  pressure. 

Test  3  presents  a  81  kPa  (12  psi)  pressure  history  which  was  generated  with  a  diaphragm 
pressure  ratio  of  P4/P1  =  28.4.  As  observed  previously  in  Section  4.2  (Figure  7c),  the  BRL-Q1D 
computation  in  Reference  17  underpredicts  the  shock  overpressure  at  the  beginning  of  the 
experimental  record.  It  seems  reasonable  to  assume  that  this  discrepancy  is  caused  by  the 
inability  of  the  one-dimensional  computational  scheme  to  simulate  two-dimensional  phenomena 
which  occur  in  the  experimental  flow,  like  a  ’’dead-water”  region  at  the  nozzle  exit.  Perhaps  the 
flow  would  be  modeled  computationally  closer  to  reality  with  the  divergent  nozzle  extended  to 
the  full  diameter  of  the  expansion  tube.  This  was  done  for  Test  3  and  the  result  in  Figure  9b 
shows  an  improved  agreement  between  computation  and  experiment. 
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(a)  Shock  Overpressure  23  kPa  (Test  20)i7 


Figure  9.  Computational  and  Experimental  Pressure  Histories  for  the  BRL-LBS  Model-A 
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(c)  Shock  Overpressure  175  kPa  (Test  12) 


(d)  Shock  Overpressure  190  kPa  (Test  25) 
Figure  9.  (continued) 
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(e)  Shock  Overpressure  202  kPa  (Test  26) 
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(f)  Shock  Overpressure  241  kPa  (Test  28) 
Figure  9.  (concluded) 
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An  altogether  different  phenomenon  is  observed  in  Test  12  which  represents  the  highest 
diaphragm  pressure  ratio,  P4/P1  =  135,  used  in  this  series  of  experiments.  The  computation 
presented  in  Reference  17  overpredicts  the  experimental  overpressure  during  the  first  ten 
milliseconds  of  the  pressure  history,  and  thereafter  the  computed  overpressure  decays  much 
more  rapidly  than  the  experimental  curve.  This  computation  was  repeated  with  the  throat  area 
reduced  to  its  actual  size  which  was  measured  on  the  burst  diaphragm  after  the  test  (Table  IV). 
The  resulting  pressure  history  (Figure  9c)  shows  much  better  agreement  with  the  experiment.  It 
may  be  concluded  therefore,  that  at  high  driver  pressures  (which  necessitate  thicker 
diaphragms),  the  blockage  caused  in  the  throat  by  the  folded  diaphragm  petals  has  to  be 
accounted  for  in  the  computational  model.  This  conclusion  is  confirmed  in  the  following 
comparisons  for  heated  driver  gas. 

4.3.2  Experiments  with  Heated  Driver  Gas 

Three  comparisons  of  computational-experimental^"  pressure  histories  for  heated  driver 
gas  are  shown.  Test  25  and  26  (Figures  9d  and  9e)  present  a  good  match  of  the  computational 
and  experimental  pressure  histories  generated  with  a  diaphragm  pressure  ratio  of  P4/P1  =  93, 
and  108  resp.,  and  a  temperature  ratio  across  the  diaphragm  of  7y7\  =  1.59.  Both 
computations  were  made  using  the  actual  throat  area  given  in  Table  IV.  The  computed 
overpressure  in  Figure  9e  indicates  the  passage  of  a  second  shock  through  the  test  station  at 
about  17  ms  which  compensates  for  the  preceding  rapid  pressure  decay.  The  plot  of  density 
versus  distance  for  Test  26  (page  77  in  the  Appendix)  reveals  that  this  shock  is  the  reflection  of 
the  recompression  shock  which  is  reflected  downstream  from  the  area  discontinuity  at  the  nozzle 
exit. 


Test  28  (Figure  9f)  represents  the  pressure  history  for  the  highest  shock  overpressure 
obtained  in  this  series  of  experiments.  It  was  generated  with  a  diaphragm  pressure  ratio  of 
P4/Px  ==  131,  and  a  diaphragm  temperature  ratio  of  TAjTx  =  1.78.  The  computation  was  done 
with  the  nominal  throat  area  because  the  actual  throat  area  was  not  known.  Consequently,  the 
computation  predicts  the  shock  overpressure  too  high  and  the  positive-phase  duration  too  short. 
These  are  clear  indications  that  the  throat  area  in  the  experiment  was  necked  down  by  the 
petals  of  a  thick  diaphragm.  Although  the  computation  with  the  actual  throat  area  would  give 
better  agreement,  the  present  comparison  is  satisfactory  in  itself. 

4.4  MATCHING  THE  CEG  PRESSURE  HISTORIES 

In  this  final  section  we  will  discuss  two  methods  of  improving  the  computational- 
experimental  mismatch.  First,  we  will  discuss  the  influence  of  the  area  discontinuity  on  the 
overpressure  spike,  then  we  will  study  the  influence  of  the  driver  volume  on  the  initial  pressure- 
decay  rate  behind  the  pressure  spike  together  with  adjustments  of  the  driver  pressure. 

4.4.1  Pressure-Spike  Study 

It  is  known  that  the  flow  in  pipes  forms  a  dead-water  region  in  the  vicinity  of  a  sudden 
diameter  change  to  facilitate  a  smooth  transition.  Extending  the  divergent  nozzle  to  the  full 
diameter  of  the  expansion  tube  replicates  such  a  dead-water  region  in  the  computational  LBS 
model.  The  effect  of  doing  this  was  studied  computationally  and  the  results  are  presented  in 
Figure  10.  Two  curves  are  shown.  The  bottom  curve  represents  the  pressure  spike  for  the 
standard  nozzle  with  the  discontinuity  in  the  cross-  sectional  area,  whereas  the  top  curve  is  for 
the  extended  nozzle.  For  the  same  initial  condition,  only  the  cone  half  angle  is  varied.  The 
result  is  a  dramatic  difference  in  the  amplitude  of  the  initial  pressure  spike.  The  extended 
nozzle  more  nearly  approximates  the  pressure  spikes  found  in  the  experimental  records  shown  in 
Figures  7  and  9. 
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HEIGHT  OF  PRESSURE  SPIKE  (k  Pa) 


This  result  clearly  shows  in  the  following  figure.  Figure  11  presents  a  comparison  of  the 
computational  shock  strength  with  the  experimental  results  as  a  function  of  the  diaphragm 
pressure  ratio.  A  computational  curve  for  the  standard  nozzle  as  well  as  for  the  extended  nozzle 
is  given.  The  extended  nozzle  clearly  improves  the  prediction  of  the  shock  overpressure  relative 
to  the  experimental  data  in  the  lower  range  of  the  driver  pressure  (10  to  100  atm).  At  driver 
pressures  above  100  atm,  the  extension  of  the  nozzle  has  no  longer  an  effect  on  the  pressure 
spike.  But  without  doubt,  real-gas  effects  become  responsible  for  the  trend  of  the  experimental 
data  below  the  computational  curve  in  this  pressure  region. 

The  computational  simulations  of  the  CEG  pressure  histories  reported  in  Subsection  4.2.2 
were  repeated  for  the  Q1D  model  with  extended  nozzle  and  the  results  are  tabulated  in  Table  V. 
Those  data  which  now  match  the  experimental  value  are  printed  in  bold  type.  It  is  evident  that 
the  extended  nozzle  causes  the  amplitude  of  the  pressure  spike  to  increase.  This  applies  to  the 
static  as  well  as  to  the  dynamic  pressure.  However,  this  increase  is  somewhat  capricious.  The 
pressure  spike  increase  is  too  much  at  the  lowest  driver  pressure  (Test  5),  and  not  enough  at  the 
next  highest  driver  pressure  (Test  9).  The  increase  is  right  only  at  the  highest  driver  pressure 
(Test  8).  The  same  capriciousness  is  observed  for  the  other  blast-wave  characteristics,  and  the 
overall  impression  is  somewhat  disappointing.  Therefore,  another  procedure  was  followed  in 
trying  to  match  the  experimental  pressure  histories  for  the  CEG  tests. 


TABLE  V. 

MATCHING  THE  EXPERIMENTAL  BLAST-WAVE  CHARACTERISTICS  FOR  CEG  TESTS 

Test 

PA 

V, 

P  ~ 

t 

I 

Y 

%o 

I 

Y 

4 

4 

‘so 

-+■ 

so 

so 

q 

q 

Nr. 

atm 

% 

kPa 

s 

kPa-s 

kT 

kPa 

kPa-s 

kT 

Computation  with  Extended  Nozzle 

5a 

10.1 

32.5 

1.35 

10.5 

17 

5.16 

.902 

21 

9 

15.4 

100 

43.6 

1.45 

16.4 

37 

6.46 

1.08 

9 

8 

26.8 

66.7 

1.41 

26.9 

88 

14.6 

2.98 

28 

Computation  with  Increased  P4  and  Reduced  V4 

5a 

11.4 

95 

30.2 

1.35 

11.3 

24 

4.54 

1.08 

50 

9 

22.9 

73 

43.5 

1.42 

17.4 

45 

6.43 

1.37 

18 

8 

37.5 

84 

62.0 

1.42 

30.5 

142 

12.7 

3.71 

78 

a  test  section  blockage  25%, 

4.4.2  Adjusting  the  Driver  Pressure  and  Volume 

The  pressure  histories  of  free- field  blast  waves  such  as  are  generated  by  HE  or  nuclear 
explosions,  can  be  characterized  ideally  by  a  three-,  or  four-parameter,  exponentially  decaying 
function5  as  e.g.  the  modified  Friedlander  equation,  or  Ethridge’s  equation.  These  equations 
allow  matching  of  three  or  four  blast-wave  characteristics,  which  are  the  peak  overpressure,  p  q, 
the  positive-phase  duration,  t  ,  the  static-overpressure  impulse,  I  and  the  initial  decay  rate, 
dp/dt|t  ,  'at  shock  arrival  time.  It  is  therefore  reasonable  to  use  experimentally  determined 
blast-wave  characteristics  for  comparing  the  computational  and  experimental  pressure  histories. 
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Figure  11.  Shock  Strength  Versus  Diaphragm  Pressure  Ratio 
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Consequently,  we  modified  the  input  to  our  Q1D  code  and  repeated  the  computations  until  the 
blast- w aye  characteristics  matched.  This  was  done  for  the  three  CEG  tests  under  discussion 
and  the  results  are  listed  in  Table  V,  also. 

The  shock  overpressure  is  adjusted  obviously  by  raising  and  lowering  the  driver  pressure, 
P4.  The  positive-phase  is  curtailed  by  the  arrival  at  the  test  station  of  the  first  expansion  wave 
from  the  open  end  of  the  shock  tube.  It  is  not  quite  so  obvious  that  the  static-overpressure 
impulse  and  the  initial  pressure-decay  rate  behind  the  shock  depend  on  the  volume  and  shape  of 
the  driver.  More  precisely,  both  impulses,  for  the  static  as  well  as  for  the  dynamic  pressure 
history,  are  proportional  to  the  mass  of  the  driver  gas  which,  for  constant  driver  pressure  and 
temperature,  is  determined  by  the  volume  of  the  driver.  The  initial  pressure-decay  rate  is 
influenced  by  a  number  of  parameters,  among  them  the  volume  and  shape  of  the  driver  and  the 
size  of  the  throat  area.  In  the  present  study,  however,  the  throat  area  and  the  driver  length  are 
considered  given  constants  so  that  the  driver  volume  could  be  adjusted  only  by  changing  the 
conical  shape  of  the  driver. 

In  our  efforts  to  match  the  experimental  blast-wave  characteristics  computationally,  we 
gave  preference  to  matching  the  static-overpressure  impulse  and  the  positive-phase  duration. 
The  initial  pressure-decay  rate  was  matched  visually  if  possible  even  when  it  became  necessary 
to  compromise  the  matching  of  the  peak  static  overpressure.  We  feel  justified  to  do  so  because 
the  overpressure  spike  varies  somewhat  with  the  number  and  location  of  the  grid  points  in  the 
computational  grid,  whereas  the  experimental  overpressure  spike  is  most  likely  affected  by 
instrumental  noise.  (The  French  scientists  at  the  CEG,  for  instance,  related  during  meetings  at 
the  BRL  that  they  ignore  the  spike  in  their  data  evaluation.)  The  resulting  computational 
pressure  histories  match  the  experimental  histories  very  well.  They  would  match  the 
experimental  records  even  better  if  the  actual  throat  areas  of  the  burst  diaphragms  could  have 
been  accounted  for  in  the  computations.  With  only  the  nominal  throat  area  knowm,  the 
computational  pressure  histories  tend  to  be  just  above  the  data  in  the  initial  part  of  the  positive 
phase  and  slightly  below  the  data  toward  the  latter  part  of  the  positive  phase. 


6.  CONCLUSIONS 


The  computational  studies  discussed  here  have  shown  that  the  BRL-Q1D  code  is  a  useful 
tool  for  predicting  shock-tube  flows  with  complicated  geometries.  The  present  quasi-one- 
dimensional  model  represents  ideal  (not  real!)  flow  and  one  may  expect  that  the  computational 
predictions  would  be  improved  by  the  inclusion  of  losses  due  to  friction  and  heat  transfer.  But 
even  without  these  losses  modeled,  the  scheme  performs  well. 

The  scheme  is  for  instance  sensitive  to  area  changes  like  the  necked-down  diaphragm 
opening  formed  by  the  folded  petals  of  the  real  diaphragm.  It  has  been  shown  here  that  a 
knowledge  of  the  necking  effect  in  the  diaphragm  (i.e.,  throat)  area  improves  the  prediction 
markedly. 

When  decaying  blast  waves  are  simulated,  relating  the  blast-wave  characteristics  of  the 
computational  pressure  histories  to  those  of  free-field  blast  waves  is  very  helpful  in  qualifying 
the  predictions. 

In  the  case  of  the  axisymmetric  LBS  model  it  was  found  that  the  flow  in  the  vicinity  of 
the  sudden  area  change  at  the  nozzle  exit  can  be  more  realistically  modeled  by  extending  the 
divergent  nozzle  to  the  full  diameter  of  the  expansion  tube.  However,  a  fully  three-dimensional 
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(a)  Test  5:  Driver  Pressure  11.4  atm,  Driver  Volume  95% 


Figure  12.  Matched  Computations  for  the  CEG  Experiments 


(b)  Test  9:  Driver  Pressure  22.9  atm,  Driver  Volume  73% 


Figure  12.  (continued) 
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(c)  Test  8:  Driver  Pressure  37.5  atm,  Driver  Volume  84% 


Figure  12.  (concluded) 


computation  is  needed  to  more  accurately  treat  the  flow  details  in  the  mixing  region  behind  the 
exits  of  the  multiple  divergent  nozzles. 

Overall,  the  BRL-Q1D  code  performs  well  in  providing  good  predictions  and  correct 
trends  with  efficiency  and  low  cost. 
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*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  * 

QUASI  1  - D  EULER  IMPLICIT/EXPLICIT  CODE 


*  *  *  *  **  **  *  *  *  *  *  *  *  *  **  *»  **  **  *  *  *  * 

THIS  CODE  COMPUTES  THE  FLOW  IN  SHOCK  TUBES  WITH  VARIABLE  CROSS- 
SECTIONAL  AREA  BY  SOLVING  THE  QUAS I  -  ONE- DIMENS IONAL  EULER  EQUATIONS 
EMPLOYING  ONE  OF  THE  FOLLOWING  FINITE- DIFFERENCE  SCHEMES. 

(A)  BEAM  &  WARMING’ S  IMPLICIT,  BLOCK- S IMULTANEOUS  SOLUTION, 

(B)  MACCORMACK’S  EXPLICIT,  PRED I CTOR/ CORRECTOR  TECHNIQUE. 

ALL  INPUT  VARIABLES  ARE  NOND IMENS IONAL I  ZED  IN  THE  CODE: 
DIMENSIONLESS  =  DIMENSIONAL  VARIABLE 


TAU 

= 

AAMB* TIME /REFL 

TIME 

X 

— 

X/REFL 

PHYS.  LENGTH 

XI 

= 

XI /REFL 

COMPUT .  LENGTH 

A 

= 

A/REFL  *  *  2 

AREA 

RHO 

= 

RHO /RHO 4 

DENSITY 

U 

= 

U/AAMB 

VELOCITY 

P 

= 

P/ ( RHO  4  »  AAMB  *  *  2 ) 

PRESSURE 

E 

= 

E / ( RHO  4  * AAMB  *  *  2 ) 

ENERGY 

***  *  *  *  *  *  * 


***  OUTPUT  «**  ***  ***  ***  *** 


THIS  CODE  WRITES  TO  FILE  OUTPUT  (=TAPE6 )  AND  TO  TAPE 1 5  AND  TAPE 1 6 . 

(1)  OUTPUT  (=TAPE6 )  IS  WRITTEN  IN  SEVERAL  SUBROUTINES: 

READ IN  AND  AREA  ECHO  THE  INPUT  DATA. 

WRIT  WRITES  THE  PHYSICAL  SHOCK-TUBE  AND  GAS  PARAMETERS  AT 
A  GIVEN  TIME  STEP. 

OUTPOT  WRITES  THE  PRESSURE  VS  TIME  RECORD,  AND  RECORDS 
THE  BLAST-WAVE  PARAMETERS  (PSO , TA, PPD , I  SO , IQ)  . 


(2)  TAPE 15  IS  WRITTEN  IN  SUBROUTINE  PLOTX.  IT  CONTAINS  PRESSURE  & 
DENSITY  VS  DISTANCE  RECORDS  FOR  PLOTTING  WITH  PROGRAM  PX1DPLT. 
DATA  ARE  WRITTEN  TO  TAPE 1 5  IN  BINARY  FORMAT  FOR  EACH  TIME  STEP 
FOR  WHICH  N/NPLOX  IS  A  WHOLE  NUMBER.  TAPE 15  HAS  TO  BE  READ 
EXACTLY  THE  WAY  IT  WAS  WRITTEN  IN  SUBROUTINE  PLOTX,  NAMELY 


READ (15) 

+ 


READ (15) 
READ (15) 
READ (15) 
READ (15) 
READ (15) 


JM1  , KCS , KTUBE , LS ,  MTRX , NPL , NSTA , PRAT , TIME , DREF , REFL , 
TEST , VOLD , XCL ST , XSTA 

(NOTE  THAT  XSTA  IS  AN  ARRAY  OF  DIMENSION^)  !  ) 
(X(  J)  ,  J=1  ,  JMl  ) 

(FP( J)  , J=1 ,  JMl ) 

(Ql(J), J=1  , JMl) 

( FU ( J )  , J=1  , JMl  ) 

(ARA(J)  ,  J=1  ,  JMl  ) 


(3)  TAPE16  IS  WRITTEN  IN  SUBROUTINE  OUTPOT.  IT  CONTAINS  PRESSURE 
&  DENSITY  VS  TIME  RECORDS  FOR  PLOTTING  WITH  PROGRAM  PT1DPLT. 

DATA  ARE  WRITTEN  TO  TAPE  16  IN  BINARY  FORMAT  FOR  EVERY  IOUT-TH 
TIME  STEP.  TAPE  16  HAS  TO  BE  READ  EXACTLY  THE  WAY  IT  WAS  WRITTEN 
IN  SUBROUTINE  OUTPOT,  NAMELY 
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c  READ  (  16)  DVNL  ,  DRVL  ,  PRAT ,  TRAT ,  TSTA ,  JMAX ,  NMAX ,  KTUBE  ,  REFL  ,  VOLD  , 

C  +  PAMB  ,  TAMB  ,  RWEL  ,  PPD  ,  PSO  ,  QIM,  QS  ,  S1M,TAR,WYS  ,WYQ  ,  PPF  ,  NXK ,  TEST ,  VALV 

C  DO  NN  M=1  ,  NMAX 

C  NN  READ  (16)  TIME  (M)  ,  VFLO  (M)  ,  OPS  (M)  ,  PDYN(M)  ,  ARJX(M) 

C 

C 

C  *  *  *  *  *  *  *  *  INPUT  REQUIREMENTS  *  *  *  *  *  *  *  *  * 

C 

C  ALL  DIMENSIONS  MUST  ACCORD  WITH  THE  INTERNATIONAL  SYSTEM  OF  UNITS, 

C  I.  E.  METRE  (M)  FOR  LENGTH,  KILOGRAM  (KG)  FOR  MASS  (WEIGHT), 

C  DEGREES  CELSIUS  (C)  FOR  TEMPERATURE,  KILOPASCAL  (KPA)  FOR  PRESSURE, 

C  MILLISECONDS  (MS)  FOR  TIME,  METRE/SECOND  (M/S)  FOR  VELOCITY, 

C  CUBIC  METRE  (M**3)  FOR  VOLUME,  AND  KG/M**3  FOR  DENSITY. 

C 


c 

ri 

TO  CONVERT  FROM 

TO 

MULTIPLY 

BY 

\j 

c 

FAHRENHEIT  (TF) 

CELSIUS  (TG) 

TO=( 5/9 ) * (TF- 32 ) 

c 

FOOT ( /SECOND) 

METRE ( / SECOND ) 

.  3048 

E+0  0 

c 

I NCH( /SECOND) 

METRE ( /SECOND) 

.  0  254 

E+00 

c 

KILOMETRE /HOUR 

METRE/ SECOND 

.27778 

E+0  0 

c 

KNOT  (NAUT.  MILE/HOUR) 

METRE/ SECOND 

.51444 

E+0  0 

c 

LBM  (POUND  MASS) 

KILOGRAM 

.  45359 

E+00 

c 

LBM/FT*  *  3  (DENSITY) 

K I LOGRAM/METRE *  *  3 

.16018 

E+0  2 

c 

LBM/ IN* *3  (DENSITY) 

KILOGRAM/METRE  *  *  3 

.27680 

E+0  5 

c 

MILE,  NAUTICAL 

METRE 

.  1852 

E+0  4 

c 

MILE,  STATUTE 

METRE 

.  1  609  3 

E+0  4 

c 

MILE/HOUR 

METRE/ SECOND 

.44704 

E+0  0 

c 

POUND  MASS  (LBM) 

KILOGRAM 

.  45359 

E+  0  0 

c 

PS  I  ( POUND/ INCH*  *  2 ) 

KILOPASCAL 

.68948 

E+0  1 

c 

RANKINE  (TR) 

CELSIUS  (TC) 

TO=(  5/9 ) *TR-  273 . 1  5 

c 

SLUG 

KILOGRAM 

.  1  4594 

E+0  2 

c 

SLUG/FOOT*  *3 

KI LOGRAM/METRE ** 3 

.51538 

E+0  3 

c 

YARD 

METRE 

.9  144 

E+0  0 

c 

c 

c 

*  *  *  *  *  *  *  *  *  INPU 

T  DATA  FOR 

MAT  *  *  * 

*  *  *  *  *  * 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ALL  INPUT  IS  LIST  DIRECTED.  DATA  ARE  STORED  ON  CARD  IMAGES  EACH 
CONSISTING  OF  A  STRING  OF  VALUES  WHICH  ARE  DESCRIBED  BELOW. 

THE  VALUES  ARE  SEPARATED  BY  ONE  OR  MORE  BLANKS,  OR  BY  A  CONMA  OR 
SLASH,  EITHER  OF  WHICH  CAN  BE  PRECEEDED  OR  FOLLOWED  BY  ANY  NUMBER 
OF  BLANKS.  ALSO,  A  LINE  BOUNDARY  SERVES  AS  A  SEPARATOR. 

EMBEDDED  BLANKS  ARE  NOT  ALLOWED  IN  INPUT  VALUES  EXCEPT  CHARACTER 
VALUES  AND  COMPLEX  NUMBERS .  WHEN  THE  VALUE  SEPARATOR  IS  A  SLASH, 
REMAINING  LIST  ELEMENTS  ARE  TREATED  AS  NULLS  AND  THE  REMAINDER 
OF  THE  CURRENT  RECORD  IS  DISREGARDED. 

INPUT  VALUES  MUST  CORRESPOND  IN  TYPE  TO  VARIABLES  IN  THE  I/O  LIST. 
HOWEVER,  REAL  NUMBERS  MAY  LOOK  LIKE  INTEGERS.  THE  FORMAT  OF  VALUES 
IN  THE  INPUT  RECORD  IS  AS  FOLLOWS. 


INTEGERS 
REAL  NUMBERS 


FORMAT  IS  SAME  AS  FOR  INTEGER  CONSTANTS. 

ANY  VALID  FORTRAN  FORMAT.  THE  DECIMAL  POINT 
MAY  BE  OMITTED  IT  IS  ASSUMED  TO  BE  TO  THE 
RIGHT  OF  THE  MANTISSA. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CHARACTER  VALUES 


LOGICAL  VALUES 


***  *  *  *  *  *  *  IN 


C**  CARD  1:  PROGRAM  I/O 


A  STRING  OF  CHARACTERS  ENCLOSED  BY 
APOSTROPHES  . . .  CAN  ONLY  BE  READ  INTO 
CHARACTER  ARRAYS,  VARIABLES  AND  SUBSTRINGS. 

T  OR  F  FOLLOWED  BY  OPTIONAL  CHARACTERS 
EXCLUSIVE  OF  SEPARATORS  (,  AND  /)  OPTIONALLY 
PRECEEDED  OR  ENCLOSED  BY  PERIODS. 


PUT  DATA  CARDS  *  *  *  *  *  *  *  *  * 


CONTROL  PARAMETERS  (7  VALUES) 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1.  TEST  -  PROBLEM  IDENTIFICATION  CODE  (USE  FORMAT  A10) 

***  NOTE:  THE  CHARACTER  STRING  HAS  TO  BE  ENCLOSED 

BY  DELIMITERS,  E.  G.  ’A  NEW  TEST’. 


2.  METH  -  CHOICE  OF  COMPUTATIONAL  SCHEME 

=*  1,  BEAM-WARMING  IMPLICIT  TECHNIQUE 
-  2,  MACCORMACK’S  EXPLICIT  TECHNIQUE 

3.  I  OUT  -  FLAG  FOR  PRINTED  OUTPUT  OF  FLOW  PARAMETERS  VS  TIME 

AT  A  PREDEFINED  X- STATION  (XSTA) . 

FOR  I OUT  >  0,  EVERY  IOUT-TH  TIME  STEP  WILL  BE  PRINTED 
FOR  I OUT  <=0,  NO  PRINTED  OUTPUT ( ! ) 

4,.  I  PLOT  -  LOGICAL  FLAG  FOR  STORING  THE  PRES  SURE  -  TIME  RECORD 
ON  TAPE16  FOR  LATER  PLOTTING:  T(RUE)  OR  F(ALSE). 

5.  NWRIT  -  FLAG  FOR  PRINTED  OUTPUT  OF  PHYSICAL  AND  FLOW  PARAMETERS 
AT  ALL  GRID  POINTS,  AFTER  EVERY  NWRIT- TH  TIME  STEP. 

USE  SPARINGLY,  E.G.  SET  NWRIT  =  NMAX,  NMAX/2,  OR 
NMAX/3 ,  ETC.  SET  NWRIT  <=  0 ,  IF  NO  OUTPUT  WANTED. 


6.  NSTA  -  NUMBER  OF  FIXED  STATIONS  AT  WHICH  PARTIAL  OUTPUT 
IS  WANTED,  FROM  1  TO  6 . 


7.  NPLOX  -  FLAG  FOR  PLOTTING  P-X  DATA  EVERY  NPLOX-TH  ITERATION 
<=0 ,  FOR  NO  PLOTS . 


C*  * 

C 

C 

C 

C 

c 

C 


CARD  2:  OUTPUT  STATIONS  ("NSTA”  VALUES,  MAX.  6) 


XSTA(J)  -  DISTANCE  OF  THE  J - TH  OUTPUT  STATION  FROM  THE 
UPSTREAM  END  OF  THE  DRIVEN  TUBE,  METRE. 


C*  * 

C 

C 

C 

C 

C 

C 

c 

c 


CARD  3  :  COMPUTAT I ONAL  -GRID  PARAMETER  S  (  4  VALUE S  ) 


1.  JMAX  -  TOTAL  NUMBER  OF  GRID  POINTS  (MAXIMUM  1,0  02  ) 

(1)  FOR  A  SIMPLE,  LINEAR  GRID  WITH  CONSTANT  CELL  SIZE 

JMAX  =  REF . -LENGTH/CELL- SIZE  +  2 

(2)  FOR  A  LINEAR  GRID  WITH  FINE- GRID  OVERLAY  ABOUT 
REGIONS  OF  AREA  CONSTRICTIONS  (BAFFLES  OR  NOZZLE) 

JMAX  =  REF  . -LENGTH /CELL- SIZE  +  BET*NXK 
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c 

c 

c 

c 

c 

C  2 .  MTRX 

C 

C 

C 

C 

c 

c 

c 

c 

C  3 .  BET 

C 

C 

C 

C 

C  4.  XCLST 

C 

C 

C 

C 

C  *  *  CARD  4  : 

0  . 

C 

C  1 .  NMAX 

C 
C 
C 

C  2 .  NUST 

C 

c 

C  3  .  CN 

C 
C 
C 

C  4  .  SMU 

C 

c 

c 

C  5  .  LS 

C 

c 

c 

c 

C**  CARD  5: 


C 

C 


c 

1 . 

PAMB 

c 

c 

2  . 

TAMB 

c 

c 

3  . 

PRAT 

c 

c 

4  . 

TRAT 

c 

c 

5  . 

P51 

(3)  FOR  NON-LINEAR  CLUSTERING  OF  GRID  POINTS  ABOUT  A 
SINGLE  POINT  USING  SINK  FUNCTION,  CONSULT 
CLUSTERING  GRAPH  FOR  CHOICE  OF  CELL  SIZE  AND  BET. 
FOR  FURTHER  INFO,  SEE  MTRX. 

-  FLAG  FOR  CHOICE  OF  GRIDDING  FUNCTION. 

=  1,  LINEAR  GRID-DISTRIBUTION  FUNCTION  USED, 

CELL  WIDTH  EQUALS  REFL /  (  JMAX- 2  )  . 

=  2,  LINEAR  MULTIPLE-GRID  FUNCTION  USED  FOR  DEFINING 
FINE  GRID  ABOUT  PARABOLIC  AREA  CONSTRICTIONS, 
SHORT  DRIVER,  NOZZLE,  AND  RWE . 

=  3,  ARCS  I NH  FUNCTION  USED  FOR  CLUSTERING  ABOUT  XCLST 
NOTE:  ONLY  ONE  CLUSTER  REGION  IN  GRID  POSSIBLE. 

-  FACTOR  USED  FOR  GRID  CLUSTERING  (E.G.  BET=2  0 )  . 

FOR  MTRX=2 :  EQUALS  NO.  OF  POINTS  IN  PARABOLIC 

AREA  CONSTRICTIONS. 

FOR  MTRX=3 :  AUGMENTS  CLUSTER  EFFECT  ABOUT  XCLST. 

-  CLUSTER  ABOUT  PHYSICAL  X- LOCATION,  IN  METRE, 

MEASURED  FROM  DRIVER  END  OF  SHOCK  TUBE. 

USED  ONLY  WHEN  MTRX  =  3  . 


COMPUTATIONAL- SCHEME  PARAMETERS  (5  VALUES) 


-  NUMBER  OF  TIME  STEPS  IN  COMPUTATION 
NOTE:  ONLY  6000  TIME  STEPS  WILL  BE  STORED 

FOR  OUTPUT  AND  LATER  PLOTTING. 

-  NEW  TIME  STEP  COMPUTED  AFTER  EVERY  NUST- TH  ITERATION. 
<=0,  FOR  CONSTANT  TIME  STEP. 

-  COURANT  NUMBER 

=  1.0,  FOR  IMPLICIT  SCHEME  (METH=1  ) 

=  0.0,  FOR  EXPLICIT  SCHEME  (METH=2) 

-  SMOOTHING  COEFFICIENT  FOR  DISSIPATIVE  TERM 
=  0.04,  FOR  IMPLICIT  SCHEME  (METH=1  ) 

=  0.02,  FOR  EXPLICIT  SCHEME  (METH=2 ) 

-  NUMBER  OF  CELLS  OVER  WHICH  THE  SHOCK  IS  STRETCHED 
ON  EITHER  SIDE  OF  THE  DISCONTINUITY, 

=  0  FOR  STEP  SHOCK,  >0  FOR  STRETCHING. 


INITIAL  (TEST)  CONDITIONS  (8  VALUES) 


-  AMBIENT  PRESSURE  [KPA] 

-  AMBIENT  TEMPERATURE  [KELVIN] 

-  RATIO  OF  DRIVER/DRIVEN- TUBE  PRESSURES 

-  RATIO  OF  DRIVER /DRIVEN -TUBE  TEMPERATURE 

-  PRESSURE  RATIO  IN  THERM\L  REGION  5 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c** 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


6  .  T5 1  -  TEMPERATURE  RATIO  IN  THERMAL  REGION  5 

7.  XTH1  -  START  OF  THERMAL  REGION  (MIN  X-  COORD  I  NATE )  ,  METRE, 
MEASURED  FROM  THE  EXIT  OF  THE  DIVERGENT  NOZZLE. 


8.  XTH2  -  END  OF  THERMAL  REGION  (MAX  X- COORDINATE)  ,  METRE, 
MEASURED  FROM  THE  EXIT  OF  THE  DIVERGENT  NOZZLE. 


CARD  6:  SHOCK-TUBE  CONFIGURATION  (10  VALUES) 


} .  KTUBE  -  TYPE  OF  SHOCK  TUBE 

=  1,  CYLINDRICAL  SHOCK  TUBE  WITH,  OR  WITHOUT  PARABOLIC 
AREA  CONSTRICTIONS  (BAFFLES)  IN  DRIVER,  DIAPHRAGM 
LOCATION,  OR  DRIVEN  TUBE  (I.E.  TEST  TUBE). 

NOTE  THAT  A  DRIVER  DIAMETER  OTHER  THAN  THAT  OF  THE 
DRIVEN  TUBE  IS  PERMITTED  WITH  THIS  CONFIGURATION. 

=  2,  SHOCK  TUBE  AS  DESCRIBED  ABOVE  WITH  A  CONVERGENT 
-DIVERGENT  NOZZLE. 

=  3,  SHOCK  TUBE  WITH  STEPPED  DRIVER  WITH/WITHOUT 

CONVERGENT- DIVERGENT  NOZZLE  AND/OR  PARABOLIC  AREA 
CONSTRICTIONS  (BAFFLES). 

«  4,  SHOCK  TUBE  WITH  CONICAL  DRIVER  AND  WITH /WITHOUT 

CONVERGENT- DIVERGENT  NOZZLE  AND/OR  PARABOLIC  AREA 
CONSTRICTIONS  (BAFFLES) . 

ALL  OTHER  VALUES  DEFAULT  TO  1 

2.  VALV  -  LOGICAL  FLAG  FOR  BLAST  VALVE  IN  DIAPHRAGM  LOCATION 

(TRUE  OR  FALSE) 

3.  NXK  -  NUMBER  OF  PARABOLIC  AREA  CONSTRICTIONS  (FROM  0  TO  11). 

NOTE:  UP  TO  1 1  CONSTRICTIONS  MAY  BE  PLACED 

ANYWHERE  IN  THE  SHOCK  TUBE  EXCEPT  IN  THE  DIVERGING 
OR  CONVERGING  SECTIONS  OF  THE  NOZZLE. 

4.  IRWE  -  FLAG  FOR  RAREFACTION -WAVE  ELIMINATOR  (RWE)  . 

<  0,  NO  RWE  (ANY  NEGATIVE  VALUE) 

=  0,  STATIC  RWE  (ZERO  ONLY) 

>  0,  MOVING  RWE  (ANY  NON- ZERO,  POSITIVE  VALUE) 

5.  REFL  -  OVERALL  REFERENCE  LENGTH  OF  SHOCK  TUBE  [METRE] 

6.  ERRL  -  PERMISSIBLE  ERROR  IN  CHECK  OF  TOTAL  LENGTH  [METRE] 

7.  DREF  -  REFERENCE  DIAMETER  OF  SHOCK  TUBE,  I.E. 

DIAMETER  OF  DRIVEN  TUBE,  METRE. 

8.  DVNL  -  LENGTH  OF  DRIVEN  TUBE  [METRE] 

9.  DRVL  -  TOTAL  LENGTH  OF  DRIVER  (DIAPHRAGM  LOCATION),  METRE 

10.  DTHR  -  DRIVER  DIAMETER  AT  THE  DIAPHRAGM  LOCATION,  METRE, 

USED  WHEN  KTUBE  =  1  . 
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c  *  * 

c 

c 

c 

0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c*  * 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c*  * 

c 


CARD  7 : 
NOTE 

1  .  XBEG 

2  .  XEND 

3 .  DBEG 

4.  DMID 

5 .  DEND 

6  .  TV  1 

7  .  TV2 

8  .  TV  3 
G  .  TV  4 


BLAST  VALVE  (9  VALUES) 

THIS  CARD  ONLY  IF  VALV=  .TRUE. 


ON  CARD  6 ! 


X- COORD  I  NATE  AT  WHICH  THE  BLAST  VALVE  (PARABOLIC 
CONTRACTION)  STARTS,  MEASURED  FROM  THE  UPSTREAM  END 
OF  THE  SHOCK  TUBE,  METRE. 

X- COORDINATE  AT  WHICH  THE  BLAST  VALVE  (PARABOLIC 
CONTRACTION)  ENDS,  MEASURED  FROM  THE  UPSTREAM  END 
OF  THE  SHOCK  TUBE,  METRE. 

DIAMETER  OF  INITIAL  BLAST  VALVE  OPENING,  METRE. 

DIAMETER  OF  FULLY  OPEN  BLAST  VALVE,  METRE. 

DIAMETER  OF  FINAL  BLAST  VALVE  OPENING,  METRE. 

TIME  AT  WHICH  BLAST  VALVE  BEGINS  TO  OPEN,  SECONDS, 
MEASURED  FROM  START  OF  THE  COMPUTATION  (ZERO  TIME). 

OPENING  TIME  OF  BLAST  VALVE,  STARTING  AT  TV1 ,  SECONDS. 

TIME  PERIOD  DURING  WHICH  BLAST  VALVE  STAYS  FULLY  OPEN, 

CLOSING  TIME  OF  BLAST  VALVE  TO  DEND,  FOLLOWING  TV 3  , 
SECONDS . 


CARD  8:  AREA  CONSTRICTIONS  (”NXK”  CARDS,  4  VALUES  EACH),  METRE. 
NOTE: 


THIS  CARD  ONLY  IF  NXK  >  0  ON  CARD  6 
MAXIMUM  OF  11  CARDS  PERMITTED. 


1.  DK( J )  -  THROAT  DIAMETER  OF  THE  J - TH  AREA  CONSTRICTION. 

2.  XI  ( J )  -  X-COORDINATE  (MIN)  AT  WHICH  THE  CONSTRICTION  STARTS 

MEASURED  FROM  THE  UPSTREAM  END  OF  THE  SHOCK  TUBE. 

3.  X0(J)  -  CENTER  POINT  (THROAT)  X-COORDINATE  OF  CONSTRICTION 

MEASURED  FROM  THE  UPSTREAM  END  OF  THE  SHOCK  TUBE. 

4.  X2(J)  -  X-COORDINATE  (MAX)  AT  WHICH  THE  CONSTRICTION  ENDS 

MEASURED  FROM  THE  UPSTREAM  END  OF  THE  SHOCK  TUBE 


CARD  9:  COMPONENT  DRIVER  DIMENSIONS  (8  VALUES),  METRE. 

THIS  CARD  ONLY  IF  KTUBE  >  1  ON  CARD  6! 

\ 

>  DRIVER  SECTION  1 


c 

c 

NOTE: 

c 

c 

1 . 

XDR 1 

-  LENGTH 

c 

c 

2  . 

DDR1 

-  DIAMETER 

c 

c 

3  . 

XDR  2 

-  LENGTH 

c 

c 

4  . 

DDR2 

-  DIAMETER 

\ 


>  DRIVER  SECTION  2 
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c 

5  . 

XDR3  - 

LENGTH 

\ 

c 

c 

c 

c 

6  . 

DDR3  - 

DIAMETER 

> 

/ 

DRIVER  SECTION  3 

7  . 

XDR4  - 

LENGTH 

\ 

c 

c 

n 

8  . 

DDR  4  - 

DIAMETER 

> 

/ 

DRIVER  SECTION  4 

Kj 

c 

c 

NOTE :  FOR 

AND 

LESS  THAN  4  SECTIONS  USE  ANY  2  OR  3  SECTIONS 
SET  REMAINDER  EQUAL  TO  ZERO. 

c 

c 

c*  * 

c 

c 

c 


c 

c 

c 

c 

c*  * 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c  > 


CARD  10;  COMPONENT  NOZZLE  DIMENSIONS  (6  VALUES),  METRE. 

NOTE-  THIS  CARD  ONLY  IF  KTUBE  >  1  ON  CARD  6! 

LENGTH  OF  THE  CONVERG I NG- NOZZLE  SECTION. 

LENGTH  OF  THE  CYLINDRICAL  THROAT  SECTION. 

DIAMETER  OF  THE  THROAT  SECTION. 

LENGTH  OF  THE  DIVERGENT- NOZZLE  SECTION. 

DIAMETER  OF  THE  D I VERGENT- NOZZLE  EXIT. 

LOCATION  OF  DIAPHRAGM  MEASURED  FROM  END  OF  DRIVER, 
(I.E.  BEGINNING  OF  NOZZLE)  WHEN  CARD  9  IS  SPECIFIED 
NOTE:  I F  ( KTUBE . LT . 3 )  THEN  XDPH  =  0  ! ! ! 

CARD  11:  RAREFACTION  WAVE  ELIMINATOR  -  RWE  (7  VALUES). 

NOTE:  THIS  CARD  ONLY  IF  IRWE  >=  0  ON  CARD  6! 

LENGTH  OF  THE  RAREFACTION  WAVE  ELIMINATOR  (RWE)  . 


c 

1 . 

XCON 

c 

c 

2  . 

XTHR 

c 

c 

3  . 

DTHR 

c 

c 

4  . 

XDIV 

c 

c 

5  . 

DDI V 

c 

c 

6  . 

XDPH 

1  .  RWEL 


2.  AFIR  -  RATIO  OF  INITIALLY  OPEN  RWE  -  AREA  RELATIVE  TO 

THE  TOTAL  CROSSECTI ONAL  AREA  OF  THE  DRIVEN  TUBE. 

3.  TAU 1  -  TIME  PERIOD  DIR  I NG  WHICH  AFIR  IS  CONSTANT,  SECONDS, 

MEASURED  FROM  START  OF  COMPUTAT I  ON  (ZERO- TIME)  . 

4.  TAU 2  -  CLOSING  TIME  OF  RWE,  STARTING  AT  TAU1 ,  SECONDS. 

5.  AMIN  -  RATIO  OF  THE  MINIMUM  OPEN  RWE  -  AREA  AT  TAU  2  RELATIVE  TO 

THE  TOTAL  CROSSECTI ONAL  AREA  OF  THE  DRIVEN  TUBE. 

6.  TAU 3  -  TIME  PERIOD  DURING  WHICH  RWE  STAYS  CLOSED,  SECONDS. 

7.  TAU 4  -  OPENING  TIME  OF  RWE  TO  AMAX  FOLLOWING  TAU 3  ,  SECONDS. 

8.  AMAX  -  RATIO  OF  THE  FINAL,  FULLY-OPEN  RWE  -  AREA  RELATIVE  TO 

THE  TOTAL  CROSSECTI ONAL  AREA  OF  THE  DRIVEN  TUBE. 


6  1 


INPUT  JOB  STREAM  USING  THE  CALCOMP  UTILITY  FOR  PLOTTING: 


KOPALKA, STMFZ ,T400 ,P3 . 
ACCOUNT ,**»**. 

D I SPOSE( OUTPUT, ST=MFA, *PR 


PF=HATS  8 


COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 
COMMENT . 


******************************************** 

*  HATS  TEST  AT  CEG 

*  FRENCH  LBS  -  1/7TH  MODEL: 

*  SHOT  #8  (27  NOV  81)  -  PSO=  67  KPA 


* 

* 

(1) 

STEPPED  DRIVER- 

V=  43.22 

M*  *  3 

* 

’ HATS 8 /CEG ’  L= 

4  4  M, 

P4  1= 

=  26.8 

* 

(2) 

CONICAL  DRIVER- 

V=  36.49 

M*  *  3 

* 

* 

’ HATS  8/Q1D ’  L= 

44  M, 

P4  1= 

=  37,5 

DRIVEN-  L=  310.0  METRE,  NO  RWE ! 
OUTPUT  STATION  #7:  XSTA=  121.5  METRE 


******************************************** 

ATTACH ( Q 1 DUB , I  D=OPALKA ) 

MAP (OFF) 

Q 1 DUB . 

RETURN, Q1 DUB. 

REWI ND ( TAPE 1 5 , TAPE 1 6 ) 

ATTACH ( CAL CMP , I D=D I S  SPLA ) 

ATTACH ( D I SSPLA, ID=DISSPLA) 

L I BRARY ( D I S  SPLA , CALCMP ) 

REWIND ( TAPE 1 5 ,TAPE1 6 ) 

ATTACH ( PXPLT , I D=OPALKA) 

PXPLT . 

RETURN ( TAPE  1 5 , PXPLT ) 

ATTACH ( PTPLT , I D=OPALKA) 

PTPLT . 

BEG IN, PLOT, CALCOMP, TAPE 7 . 

*EOR 


’HATS 8 /CEG’ 

1 

30 

T 

-  99  9  9  1 

60 

66.5 

50  2  0 

27  00 

1  1 

0 

3 

.  04 

5  45 

.  65 

0 

0 

99.0 

280  .  1  5 

26.8 

1 

1  1 

0  0 

3  F 

0  -  1 

3  65.0 

0  1 

3.5886 

3  10.0 

9  .711  9 

.87  1 

7  1 

.  124 

19  1.33 

0.89 

2.52 

0  .  665 

7  .  59 

2.286 

1.65 

’ HATS 8 /Q1D ’ 

1 

30 

T 

■  -  9 

9  99  1 

60 

66.5 

502  0 

27  0  0 

1  1 

0 

3 

.  04 

5  45 

.  65 

0 

0 

99.0 

2  80  .  1  5 

37.5 

1 

1  1 

0  0 

4  F 

0  -  1 

3  65.0 

.  01 

3 .5886 

310.0 

0  0 

0  0 

0  0  . 

66 

44 

1.33 

0.89 

2.52 

0.665 

7.59 

2.286 

1  .  65 

*  EOR 

8.5  11 

4 

47  0.01 

5  ’ 

CALCMP’ 

0  .  1  . 

.  1 

5 

.00  .5 

.  05 

5 

*  EOR 

11.  8  . 

5  2.1 

1  .  1 

7 

4  .  5 

45  .  65 


45  .  65 


0  1.75  .25  5 

-20  100  20  5 

0  6.4  2.2  .19 


’ CALCMP ’ 


0.665 


0.665 
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INPUT  JOB  STREAM  USING  THE  DISSPOP  UTILITY  FOR  PLOTTING: 


OPALKA, STMFZ ,T400 ,P3 . 

ACCOUNT ,**»**. 

DI SPO SE( OUTPUT, ST=MFA, *PR) 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

COMMENT . 

ATTACH , Q 1 DUB , I D=OPALKA . 

MAP (OFF) 

Q1DUB . 

RETURN, Q1 DUB. 

REQUEST, PLF I LE , *PF . 

ATTACH, DISSPLA, ID=DISSPLA. 

ATTACH  ,  COMPRS  ,  I  D=D  I  S  SPLA . 

L I BRARY , D I S  SPLA , COMPRS . 

REWIND ( TAPE  1 5 , TAPE 1 6 ) 

ATTACH , PTPLT , I D=OPALKA . 

PTPLT . 

SAVEPF  ( PLF  I  LE ,  PTA2  6  ,  UN=OPALKA ,  UP=*  *  *  *  *  ,  ST=MFA ,  CH=*  *  *  *  *  ,  PN=*  *  ,  TY=SAVE ) 
RETURN ( TAPE  1 6 , PTPLT ) 

REWIND, PLF I LE. 

ATTACH , PXPLT , I D=OPALKA . 

PXPLT . 

SAVEPF ( PLF I LE , RXA2 6 , UN=OPALKA , UP=*  *  *  *  « , ST=MFA , CH=* ♦ *  *  * , PN=* * , TY=SAVE ) 
*EOR 


************************************** 

*  MATCHING  BERTRAND’S  EXPERIMENTS 

*  LBS-MODEL  A2  -  SCALE  1:37.4 

*  TEST  #  24-84-26 


* 

* 

STEPPED  DRIVER- 

P  4  1= 

107.8 

* 

* 

L=  0.884  M 

T4  1= 

1.593 

* 

EXPANSION  TUNNEL - 

L=  5  . 

.461  M; 

* 

* 

TEST-STA  #  7  (X= 

1.778 

METRE) 

* 

JMAX=  502  WITH 

MTRX= 

=  3 

************************************** 


’MOD . A2 -26’  I 
1.778 


502 

0  4320 

5 

.  928 

0 

10  1.4 

299 

.  8 

3 

F  0 

-  1 

.241 

.  0503 

.024 

.068 

EOR 

1  1 

8  .  5 

2  .  1 

0 

035 

005 

-  50 

250 

50 

0 

6  .  4 

2  .  2 

1  EOR 

8  .  5 

1  1 

4 

0  . 

1  . 

.  1 

0  . 

.55 

05 

45  T  2100 

1  1.0  0  3 

0 

1  07.8  1  .  593  1 

6.64  .001  .25 

241  .0617  .187 

435  .202  .1618 

1.1  7.0  4.5 

5 

5 

.19  2  ’COMPRS’ 

49  .01  5  ’COMPRS’ 

5 
5 


90 
.  04 
0  0 

5.461  .928 

0795  .215  .0942 

.044 


0  435 
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input  parameters  for  case  no. 

SCARD1 


TEST 

*  ’  MOD  • AS— 26  *  f 

METH 

*  If 

I  OUT 

-  4Sf 

IPLOT 

*  T, 

NVTRIT 

*  21 00 i 

NSTA 

*  If 

NPLOX 

=  90? 

SEND 


SCARD2: 

tPJTPUT  (TE2T)  STATIONS 

XSTAm  *  .17780E+01 

SCARES 

05 

»u 


J  M  AX 

*  5029 

J  ADD 

*  Of 

NM  AX 

*  4320? 

MUST 

*  1? 

CN 

*  • lE+Olf 

IS 

*  Of 

MTP.X 

«  3f 

5  MU 

*  • 4E—01f 

SEND 


SCAP.D4 

BET 

*  .5E+01? 

XCLST 

=  .928E+00f 

US 

**  Om  Of 

SLPE 

*  0.  Of 

SEND 


N  o  t  e  :  STRETCHING  FUNCTION  USED  IS 

XI  =  ALPHA*  AP.CSINHt  BET*  l  X-XO  J ) 

SCARDr 


PAMB 

- 

•1014E+03, 

IAMB 

* 

• 2998E+03, 

PRAT 

•M 

•  10T8E+0  3# 

TRAT 

■= 

•1593E+0I, 

P51 

a 

.1E+01, 

T51 

a 

•  1E+01, 

TTH1 

a 

•* 

o 

• 

o 

XTH2 

SEND 

a 

0.0, 

*** 

NO  THERMAL  Li 

SCARD6 

ktube 

a 

3, 

VALV 

a 

F# 

NXK 

a 

0, 

TRVfE 

- 

-1, 

REEL 

- 

.664E+01, 

ERRL 

a 

.IE-92, 

OREF 

a 

•254E+00, 

DVNL 

a 

■5461E+01, 

DRVL 

- 

•928E+00, 

OTHR 

SEND 

SCARD7 

a 

.435E-01, 

**  NO  BLAST  VALVE  SPECIFIED!  ** 

S CARDS 

**  NO  AREA  CONSTRICTIONS  SPECIFIED,  NXK.LE.O 


SCARD9 


XDR1 


241E+00, 


00R1 

S 

.503E-01, 

XDP.2 

- 

• 241E+00* 

DDP.2 

* 

.617E-01f 

XDR3 

* 

.1 87E+00* 

DDR3 

at 

.795E-01* 

XDR4 

m 

•215E+00* 

DDR4 

at 

.942E-01# 

SEND 

SCARD10 

XCOM 

m 

.24E-01, 

XT  HR 

- 

.68E-01# 

DTHR 

* 

*  435E—01 j 

XDTV 

- 

• 202E+00? 

ODIV 

■m 

.1618E+00* 

XDPH 

- 

.44E-01, 

SEND 

SCARD11 

! 

***  MO  RVJH  SPECIFIED!  *** 
NORMALIZED  DIMENSIONS  * 


XDR1*  XDR2*  XDR3=  XDR4* 

•03630  • 07259  .10075  ,13313 

NORMALIZED  DIAMETERS  * 

DDR1*  tD DR 2*  )  { DDR3* ) 

•  7575E-02  • 9292  E— 02  .1197E-01 

OUTPUT  GAS  PARAMETERS  AT  STATIOMS): 


XCOM* 

.13675 

XDPH* 

.13976 

XTHR*  XDI V3 

.1^699  .17741 

R4* 

1419E-01 

DTHR* 

.6551E-02 

DDI V* 

• 2437E-01 

.44518 


**  GRID  CLUSTERING  ABOUT  XCLST*  .13976  ** 

THE  REFERENCE  AREA  (COMPUTED  FROM  DREF)  IS  . 5067E-01  SQUAREMETRES . 


THE  DRIVER  VOLUME  WAS  COMPUTED  FROM  AREA(J)  TO  . 3739E-02  CUBICMETRES. 


XRWE* 

.99985 


x: 


END  OF  INPUT  DATA  -  RESULTS  FOLLOW 


SSC3XSX *  r a  3 aaa  =  *** s *■ * 


**  *1  ~  0  AND  TIME t SCC)  *  .OOOOOE+OO 


"X" 

1 

10 

11 

20 

21 

30 

31 

40 

41 

50 

51 

60 

61 

^0 

71 

00 

81 

90 

100 

101 

no 

111 

120 

121 

130 

131 

140 

141 

150 

151 

160 

161 

1^0 

171 

180 

181 

190 

191 

200 

201 

210 

211 

220 

221 

230 

232 

240 

241 

250 

2  51 

260 

261 

2^0 

271 

280 

281 

290 

291 

300 

301 

310 

311 

320 

321 

330 

331 

340 

341 

350 

351 

360 

361 

370 

371 

380 

381 

390 

391 

400 

401 

410 

411 

420 

421 

430 

431 

440 

441 

450 

451 

460 

461 

470 

4“T1 

480 

481 

490 

491 

500 

501 

502 

"XIX 

1 

10 

11 

20 

21 

30 

31 

40 

41 

50 

51 

60 

-.13724E-02 
* 12197E— 01 
. 25367E— 01 
•3E173E— 01 
« 50656 E— 01 
.62S57E-01 
.74813E-01 
• 86563E— 01 
.98145E-01 
.10959E+00 
• 120950+00 

•  132240+00 
.14351E+00 
.15479E+00 
• 166120+00 
•177540+00 
. 1790^0+ 00 
•290760+00 
• 21264E+00 
• 22476E+00 
• 237140+00 
» 249840+00 
• 262830+00 
• 276310+00 
• 290180+00 
• 304530+00 
. 319390+00 
• 33483E+00 
.350890+00 
.367620+00 
•385070+00 
«  40330E+00 
.422378+00 
.442330+00 
. 46326E+00 
. 485210+00 

•  5  08260+00 
. 532488+00 
.55795E+00 
•58474E+00 
.612950+00 
• 642668+00 

•  6~*396E+00 
•70697E+00 
•74177E+00 
•778 48 E+00 
.817230+90 
•858120+00 
•90129E+00 
.94688E+00 
•995048+00 


•72980E+03 
•748240+03 
.7^03^E+03 
•791258+03 
.810660+03 
• 82835E+03 


• OOOOOE+OO 
•13532E— 01 

•  26663E— 01 

•  394350—01 

•  51888E— 01 
.64063E-01 
• 7599^8— 01 
•87729E— 01 
•99295E— 01 
•11073E+00 
• 122088+00 
•13337E+00 
.144  64E+00 
•155928+00 
•16726E+00 
.  178688+00 
•190238+00 
.201948+00 

•  21384E+00 
*225988+00 

•  238408+00 

•  25112E+0Q 

•  264208+00 
.277638+00 

•  291590  +  00 
.305998+00 
. 32091E+00 
. 33641E+00 
«  352538+00 

•  36933E+00 

•  386868+00 

•  405178+00 
• 42433E+00 
•44438E+00 
•465418+00 
. 48747E+00 
•51063E+00 
.5349^E+00 
*560578+00 
.  58750E+00 

•  61535E+00 
.  645^28+00 
• 67719E+00 

•  71 036E+00 
•74535E+00 
.78226E+00 
.8  2122E+00 
• 362338+00 
•90574E+00 
• 951 58E+00 
•100008+01 


•72863E+03 
•75050E+03 
•^7251E+03 
•79327E+03 
• 312518+03 
•  83001E  +  03 


• 13724E— 02 
.148628-01 
« 27956E— 01 
« 40694E— 01 
.531180-01 
.652660-01 
.771790-01 

•  88892E— 01 
•10044E+00 

•  11187E+00 
.123210+00 
.134500+00 
.145770+00 
.157058+00 
• 16840E+00 
. 17983E+00 
.191390+00 
.203120+00 
.215050+00 
.227210+00 
.239650+00 
. 25241E+00 
.265530+00 
.279050+00 
.293010+00 
.307460+00 
•32243E+00 

•  33**  990  +  00 
.354180+00 
.371050+00 
. 38865E+00 
.407050+00 
.4262^8+00 
.446440+00 
•46757E+00 
.489730+00 
.513010+00 
.  5374TE  +  00 
.563200+00 
.590270+00 
.618770+00 
•64879E+00 
.680430+00 
.713780+00 
•74896E+00 
,786070+00 
. 82523E+00 
.866570+00 
•91021E+00 
• 95630E+00 


.729800+03 
.752750+03 
.774650+03 
.795260+03 
• 81434E+03 
.831660+03 


.27405  E— 02 
. 16189 E-01 
.29245 E-01 
.419500-01 
.543440-01 

•  66468  E— 01 

•  78358  E— 01 
.90054 E-01 
.10X590+00 
.113010+00 
.124340+00 
.135620+00 
.146890+00 
• 15319E+00 
.169540+00 
.180980+00 
.192560+00 
.204300+00 
• 21625E+00 
.228440+00 
.240920+00 
. 25371 E+00 
.266870+00 
• 28043E+00 
. 29443  E+OC 
. 30893  E+00 
.323960+00 
.33959 E+00 
« 35  584  E+00 
, 37278  E+00 
.390460+00 
.408930+00 

•  42826  E+00 
.448510+00 
.469740+00 
. 49201 E+00 
• 51540E+00 
• 53999E+00 
.565840+00 
.593050+00 
• 62170E+00 
.651880+00 
.683680+00 
.717210+00 
.752580+00 
• 78989E+C0 
.829260+00 
.87082 E+00 
.91471 E+OC 
.96105 E+00 


•  73  214E+03 

•  75499  E+03 
•77677E+03 
•79724E+03 

•  81616  E+03 
♦ 83328E+03 


•41042E— 02 

•  175  HE— 01 
• 30530E— 01 

•  432  03E— 01 
.555  68E-01 
.676670-01 
.795360-01 
.91214E— 01 
•10274E+00 
.114140+00 

•  12547E  +  00 

•  13675E+00 
.  14C02E+00 

•  159  32E  +  00 
.170680+00 
.182130+00 
•19372E+00 
.205490+00 
.217460+00 
.229680+00 
.242180+00 
.255010+00 
. 268  ?0E  +00 
.28181E+00 
.295860+00 
• 31041E+00 
.325500+00 
.341130+00 
.357500+00 
.374510+00 

•  392  27E+00 
.410830+00 
.430240+00 
.450590+00 
.471910+00 
.494300+00 
.517810+00 

•  542  51E  +  00 
.  56850E  +  00 
.595850+00 
•62465E+00 
.654980+00 
.686960+00 
.  72066E  +  00 
.756  22E+C0 

•  793  73E+00 
.833320+00 
.875110+00 
.919230+00 
•96583E+00 


.734460+03 

.757220+03 

.778830+03 

.799210+03 

•  31796E+03 

•  83489E  +  03 


.546350-02 

.188300-01 

.318120-01 

•  44452  E— 01 
.567890-01 
.638630-01 
•30712E— 01 

•  92  373E-01 
.10388E+00 

•  11528E  +00 
.12660E+00 

•  13788E  +00 
.149150+00 
. 16045E  +00 
.171820+00 
.183280+00 

•  19489E  +00 

•  20667E  +00 
.218670+00 
.230910+00 

•  24  345E+00 

•  25631 E  +00 
.260540+00 
.283190+00 
.297290+00 

•  311 89E  +  00 
.327040+00 
.342780+00 
.359170+00 
.376250+00 
.394090+00 
.412730+00 
•43224E+00 
.452670+00 

•  47410E+00 
.496600+00 
.520220+00 
.54505E+00 

•  57117E  +00 

•  59866E  +  00 
.6276X0+00 
.658110+00 

•  69025E  +00 
•72414E+00 
•75988E+CC 
.797590+00 

•  83739E  +  00 
.879410+00 
.923780+00 
.970630+00 


•  73  678E  +03 

•  75944E  +  03 
.78  098E  +03 
.801160+03 

•  81974E  +03 

•  83647E  +03 


.681870-02 
•20145E— 01 
.330910-01 

•  4  5  69 9E— 0 1 
.58008E-01 
•70058E— 01 
•81886E— 01 
« 9  35  300-01 
.105030+00 

•  1 1642E  +  00 
.  1 2773E  +  00 

•  1 3900E  +  00 
•15028E+00 
•1 6158E+00 
.172960+00 
.184440+00 
• 19606E+00 
.207860+00 
.219880+00 
.232150+00 
.244720+00 
.257620+00 
.270890+00 
.284580+00 
.298730+00 
.313380+00 
.328590+00 
.344390+00 
.360840+00 
.378000+00 
.395910+00 
.414640+00 
.434240+00 
.454770+00 
•47630E+00 
•498910+00 
.522650+00 
.  547610  +  00 
.573860+00 
.601490+00 
•63059E+00 
.661240+00 
.693560+00 
.727620+00 
.763560+00 
.801480+00 
.841490+00 
.683740+00 
.928350+00 
.975460+00 


.739090+03 
.761650+03 
.783060+03 
.803090+03 
•  82150E  +  03 
.838040+03 


.816960-02 
.21456E-01 
. 34367E-01 
• 46943E-0 1 
.59224E— 01 

•  712  50E— 01 
«  8 30 58E— 0 1 

•  9468  6E-01 
« 1061 7E+00 
.  1175  5E  +  00 
.  128  86E  +  00 

•  1 40 1 3E  +  00 
.151400+00 
• 16272E+00 
. 1741 OE+OO 
.185590+00 
• 19723E+00 
• 20905E+00 
. 221100+00 
.233400+00 
.245990+00 
.258930+00 
•27224E+00 
.  285 97E  +  00 
.300170+00 
.314880+00 
. 330140+00 
.346010+00 
.362530+00 
.379760+00 
.397750+00 
.416560+00 
.436250+00 

•  456  880  +  00 
.478520+00 

•  501 2  3E  +  00 
.  525090  +  00 
.550170+00 
.576560+00 

•  604  3  3E  +  00 
,633580+00 

•  664  40E  +  00 
. 69688E+00 
.731130+00 
» 76726E  +  00 
. 805380+00 
.845620+00 
. 88809E+00 
.932940+00 
.980310+0 0 


.741390+03 
.763850+03 
.785130+03 
.305010+03 
. 82324E+03 
.839580+03 


.951630-02 
.227636-01 
.356390-01 
.481836-01 
•60437E— 01 

•  72  4  40E-01 
.842280-01 

•  95  8  40E-01 

•  10731E  +  00 
.118680+00 
.129990+00 
.141260+00 

•  15253E+00 
.  163  85E  +00 
.175240+00 
•18675E+00 
.198410+00 
.210250+00 
.22231E+00 
.234646+00 
.24727E+00 
.260240+00 
.273590+00 
•28737E+00 
.301620+00 
.316380+00 
.331700+00 
.347630+00 
.364220+00 

•  38 1 52E  +  00 

•  399  59E  +  00 
.418490+00 
.  438270  +  00 
.459000+00 

•  48074E  +  00 
.503560+00 
.527540+00 
.552750+00 

•  57  927E  +00 
.607190+00 
.636590+00 
.667576+00 
.700220+00 
•73466E+00 
•77098E+00 
.809310+00 
•84976E+00 
.892470+00 
.937560+00 
. 985190+00 


.743680+03 

•  766  03E  +  03 

•  78719E  +03 
.806910+03 
.824960+03 

•  84 1 1 0E  +  03 


•  i 08  59  E-0 1 
• 24067E-01 
.36907 E-01 
•49421 E-01 

•  616  4  8  E-0 1 

•  7362  7  E-0 1 

•  853  97E-0 1 

•  969  9  3  E— 0 1 

•  108  45  E  +  03 
•11982E+00 
•13111E+00 
.  14238E  +  00 
•15366E+00 

•  1 649  9  E  +  00 
.  176  39E  +  00 

•  187  9 1 E  +  00 

•  199  5  8  E  +  00 
•21144E+00 
.2 2353E+00 
.  235890  +  00 
.243556+00 
.261560+00 

•  27495  E  +  00 
.288770+00 
•30307E+00 
•31783E+00 
.33326E+00 

•  349  26E  +  00 
. 36592E+00 
, 38329E+00 

•  401 44  E+00 

•  420  4  2  E  +  00 

•  440  30E  +  00 

•  461 12  E  +  00 

•  482  976  +  00 
.50591 E+00 

•  5  3001 E  +  00 
.555340+00 

•  5  82  00  E  +  00 
.610060+00 
. 639620+00 
•67076E+00 
•70359E+00 
.  738200  +  00 
.774720+00 

•  8 13  26  E  +  00 
. 85393E+00 

•  896  87E  +  00 
. 94221 E+00 
. 990 1 0  E  +  00 


•  745  97  E  +  03 
• 76821E+03 
• 78923E+03 
.  80  8  79E  +  03 
•92666E+03 
.84261 E+03 


61 

70 

• 84409E+03 

« 84554  E+03 

• 84698E+03 

•84839E+03 

71 

90 

• 85766E+03 

«  85389E+03 

• 86009E+03 

•86I27E+03 

81 

90 

•86887E+03 

•86985E+03 

• 87081E+03 

•S71^4E+03 

oi 

100 

•87755 E+03 

•87S27E+03 

•  87897E+03 

• 87964E+03 

101 

110 

•88357E+03 

•88403E+03 

• 88445E+03 

•88485E+03 

111 

120 

•88685E+03 

•88702E+03 

.88717E+03 

•80728 E+03 

121 

130 

• 88732E+03 

• 88721E+03 

• 88707E+03 

«  88691 E+03 

131 

140 

• 88498E+03 

•88459E+03 

• 88417E+03 

•8B373E+03 

141 

150 

• 87986E+03 

• 87920E+03 

• 87851E+03 

•  87780  E+03 

151 

160 

• 87205E+03 

• 87113F+03 

•87018E+03 

•  86920  E+03 

161 

1^0 

•86167E+03 

•86050E+03 

• 85930E+03 

•  85808  E+03 

171 

180 

• 84887E+03 

•84746E+03 

• 84603E+03 

•  84458 E+03 

151 

190 

•  83383E+03 

•93222E+03 

• 83058E+03 

• 82892E+03 

191 

200 

.81677E+03 

•81496E+03 

• 81314E+03 

•81129E+02 

201 

210 

•79791E+03 

•  “r9594E+02 

• 79395E+02 

•79194E+03 

711 

220 

•77749E+03 

•77538E+03 

•77324E+03 

• 77110E+03 

221 

230 

•75575E+03 

•  ~>5352E+03 

•75127E+03 

• 74901E+03 

231 

240 

•73293E+03 

•73060E+03 

• 72826E+03 

•  72591E+-03 

241 

250 

•70926E+03 

«  70636E+03 

•70445E+03 

•  70203  E+03 

251 

260 

.6  8497E+03 

•68251E+03 

*  68005E  +  03 

•67759E+03 

261 

2^0 

• 66025E+03 

•65777E+03 

•  65528E+03 

•  65  279  E+03 

2“T1 

280 

•63531E+03 

•  63  281E+03 

• 63031E+03 

•  62781 E+03 

281 

290 

•61032E+03 

•60783E+03 

• 60533E+03 

• 60284E+03 

291 

300 

•58544E+03 

•  58296E  +  03 

• 58048E+03 

•57801E+03 

301 

310 

•56079E+03 

•55C34E+03 

• 55590E+03 

• 55346E+03 

311 

320 

•53650E+03 

•53409E+03 

•53169E+03 

• 52930E+03 

321 

330 

•51266E+03 

•51031E+03 

•  50~796E+03 

• 50562E+03 

331 

340 

•45937E+03 

•  48707E+03 

•  484“f8E+03 

•  48  250E+03 

341 

350 

•46668E+03 

•  46445E  +  03 

• 46222E+03 

• 46000E+03 

351 

360 

•44466E+03 

•44249E+03 

•44C34E+03 

.43 818 E+03 

361 

370 

• 42333E+03 

•42124E+03 

• 41916E+03 

•41708 E+03 

3"*1 

380 

•  402’74E+03 

•40C72E+03 

•  398*TlE  +  03 

• 39671E+03 

381 

39  0 

<•  38290E+03 

•  38096E+03 

• 37902E+03 

•3^710E+03 

391 

400 

•36382E+03 

•  36196E+03 

• 36010E+03 

•  35  825E+03 

401 

410 

•34551E+03 

« 34372E+03 

• 34194E+03 

•  34016  E+03 

411 

420 

• 32796E+03 

•  32625E+03 

•  32454E+03 

•32284E+03 

421 

430 

•31117E+03 

•  30953E+03 

• 30790E+03 

•  30628  E+03 

431 

440 

• 29512E+03 

•  29355E  +  03 

•29200E+03 

• 29045E+03 

441 

450 

•  2_?980E+03 

•27831E+03 

• 27682E+03 

•27534E+03 

451 

460 

• 26519E+03 

•  26377E  +  03 

•  262  36E  +  03 

•  26C95  E+03 

461 

420 

•25128E+03 

•  24992E  +  03 

«  24857E+03 

• 24723E+03 

471 

480 

• 23903E+03 

•  23674E+03 

«  23546E+03 

.  23  418  E+03 

481 

490 

• 22543E+03 

•22421E+03 

•  22299E  +  03 

•22177E+03 

491 

500 

« 21345E+03 

•  21229F+03 

•21113E+03 

• 20998E+03 

501 

502 

•  20208  E+03 

•2009^E+03 

"TUBE  AREA” 

VALUES  FROM  J  *  1  THROUGH  20  ARE  ,4507033E-04 


21 

30 

•  450“*0F— 04 

•45070E— 04 

•45070E— 04 

•45070E— 04 

VALUES  FROM 

J  - 

31 

THROUGH  50 

ARE  •6781489E— 04 

51 

60 

•67815E— 04 

.67815E-04 

•67815E-04 

•67815 E-04 

VALUE 

S  FROM 

J  * 

61 

THROUGH  80 

ARE  .112587 

2E-03 

81 

90 

.11259E-03 

•11259E-03 

•11259E— 03 

• 15807E— 03 

VALUE 

S  FROM 

J  * 

91 

THROUGH  110 

ARE  •1580726E— 03 

111 

120 

•15807E— 03 

• 14715E— 03 

•10040E-03 

•  62568  E-04 

121 

130 

•33708E— 04 

.33708E-04 

• 33708E— 04 

.33708 E-04 

131 

140 

•9T155E-04 

• 10906E— 03 

•12167E— 03 

. 13496 E-03 

141 

150 

•24765E— 03 

• 2665^E— 03 

• 28620E— 03 

•  30655  E— 03 

VALUES  FROM  J  =*  151  THROUGH  502  ARE  .1149270E-02 

”APr  =  DA/DX” 

VALUES  FROM  J  «  1  THROUGH  20  ARE  .OOOOOOOE+OO 

21  30  ,000000+00  .OOOOOE+OO  .OOOOOE+OC  .OOOOOE+OO 

VALUES  FROM  J  *  31  THROUGH  50  ARE  .OOOOOOOE+OO 


84978E+03 
862  43  E+03 
872  65E+03 
88028E+03 
895  22E+03 
88737E+03 
886  72E+03 
883  26E+03 
877  06E+03 
86820E+03 
85683E+03 
84311E+03 
827  24E+03 
80943E+03 
73992E+03 
76894E+03 
74674E+03 
72355E+03 
69961E+03 
67512E+03 
650  29E+03 
6  25  31E+03 
600  35E+03 
575  54E  +  03 
551 02E+03 
526  90E+03 
50328E+03 
480  22E+03 
45779E+03 
436  04E+03 
415  01E+03 
39471F+03 
375 18E+03 
35640E+03 
3  38  40E  +  03 
32115E+03 
304  66E+03 
28890E+03 
27387E+03 
259  54E+03 
24590E+03 
232  91E  +  03 
22C57E+03 
20883E+03 


45070E-04 

678 15E— 04 

158C7E-03 

33708E-04 
40217E— 04 
14896E-03 
327  61E-03 


000  OOE+OO 


«  85115 E  +03 

•  86357E  +03 

•  87353E  +03 

•  88090E  +  03 
« 885 56E  +03 
« 88743E  +03 
•88650E+03 

•  88276E  +03 

•  87629E  +03 

•  86718E  +  03 

•  85  556E  +03 
•84162E+03 
• 82554E+03 

•  80755  E  +03 
• 78788E+03 

•  76677E  +03 
•74446E+03 
*721 19E  +03 

•  69718E  +03 

•  67265E  +03 

•  64780E  +03 
•62282E+03 

•  59786E  +  03 
«  57307E  +03 
• 54859E+03 
«  52452E  +03 
• 50094F+03 

•  47795E  +03 

•  45558E  +03 
•43391E+03 

•  41294E  +  03 

•  39273E  +03 

•  37327E  +03 
•35457E+03 
•33664E+03 
« 31947E  +  03 

•  30305E  +03 

•  28737E  +03 
•27241E+03 
•25815E+03 

•  24457E  +03 
•23165E+03 

•  21937E  +03 
•20769E+C3 


•45070E— 04 

•  678 15  E  — 04 
•15807E-03 

« 33708E— 04 

•  47985  E-04 
•16365E-03 

•  34940E-03 


•  OOOOOE+OO 


.35250E+03 

•  8646  8E  +03 

•  87439E  +  0? 
•88149E+03 

•  88587E  +  03 

•  8  3747E  +  03 
•88625E+03 
•88224E+03 

•  87549E  +  03 

•  8  661 2  E  +  03 
•85427E+03 
•84010E+03 
.82383E+03 

•  80566E  +  03 
•78583E+03 
•76459E+03 
•74217E+03 
•71882E+03 
•69475E+03 

•  67018F  +  03 
•64531E+03 

•  62032E  +  03 

•  5  9537E  +  03 
•57061E+03 

•  546 16E  +  03 
•52214E+03 

•  4  9862E  +  03 
.A7568E+03 
•45339E+03 
♦43178E+03 
.  4 1089E  +  03 
•39075E+03 

•  371 36E  +  03 
•35274E+03 
•33489E+03 
« 3 1779E  +  03 
•30145E+03 
•28584E+03 
•27095E+03 

•  25676E  +  03 

•  24325E  +  03 

•  23039E  +  03 
•21817E+03 

•  20656E  +  03 


. 45070E— 04 
.  6781 5E-04 

•  1  5807E-03 

•33708E-04 

•  5  6440E— 04 
•17904E-03 

•  3719  OF— 03 


.OOOOOE+OO 


•85382F+03 

•  865  76E  +  03 
• 87522E+03 
• 88205E+03 
.  88616E  +  03 

•  88747E  +  03 
.  38597E  +  03 
.  98168E  +  03 
• 87467E+03 
• 36505E+03 
•85295E+03 

•  838  57F  +  03 
.  822  09E  +  03 
. 90375E+03 
• 78377E+03 
•76240F+03 

•  7  398  3E  +  03 
.7164^6+03 
•6923IE+03 
•66770E+C3 

•  642  81F  +  03 

•  61782E  +  03 
•59288E+03 

•  568 1 5E  +  03 
• 54374E+03 
. 51976E+03 
. 496  30F  +  C3 

•  473  42E  +  03 
•45119 E+ 03 
. 42965E+03 
•40884E+03 
. 38877E+03 
• 36946E+03 
•35092E+03 

•  333 14F  +  03 
•3161 3E+03 
•29986E+03 
. 28432E+03 
.  269  5  0E  +  03 
•25538E+03 
• 24194E+03 
•22914E+03 
• 21698E+03 
. 20543E+03 


•  4507 OE- 04 
• 67815E-04 

•  1 58  07E  —  03 

•33708E-04 
• 65583E-04 

•  195148-03 

•  3951 3E-03 


• 00000c+00 


•85512E+03 

•86682E+03 

•87602E+03 

•88259E+03 

•88642E+03 

•  88745E  +  03 
. 88567E+03 
.88110E+03 
•87383E+03 

•  86  3  9  5E  +0  3 

•  8  5  1 6  IE  +0  3 

•  837  018 +03 
•82034E+03 
•80182E+03 
•78169E+03 
• 76020E+03 

•  73  7  57E  +  03 
•71405E+03 
•68986E+03 
•66522E+03 

•  64  0  31c  +03 

•  61  532E  +  03 

•  59040E  +  03 
•5656 9E +03 
•54132E+03 
•51739E+03 
•49398E+03 

•  471 17E  +  03 
•44901E+03 

•  42754E  +  03 
. 40680E+03 
.  38681E  +  03 
•36758E+03 

•  34  9  1  IE  +  03 
•33141E+03 
•31447E+03 
• 29827E+03 
•28281E+03 
•26806E+C3 

•  2  5401E  +  03 
. 24063E+03 

•  22790E  +03 
. 21580E+03 
•20431 E +03 


•  450708—04 

•  67  8 1 5E— 0  4 
•15807E-03 

•33708t-04 

•  75  4178—04 
•21194E-03 

•  4 1 9  08E  — 0  3 


.113726-04 


.  856406  +  03 
. 367d6F+03 
.876308*  >3 
•883096+33 
. 886656+03 
•887408+03 
.885348+03 
•38050E+33 
•87295E+03 
•86282E+03 
. 85025E+03 
• 83543E+03 
.818566+03 
•79987E+03 
.779606+03 
. 75798E+03 
•735266+03 

•  71166E  +  03 
•68742E+03 
•66274E+03 

•  6373 1 E  +  03 

•  61 2  82  E  +  0  3 
•58792E+03 
• 56324E+03 

•  538  9 1 E  +  03 
. 51502E+03 
•49167E+03 
.468926+03 
•44683E+03 
.425436+03 
•404778+03 
•38435 E+03 
•36569E+03 

•  347  3 1 E  +  03 

•  32968  E  +  03 
. 312  3 1 E  +  03 
• 29669E+03 
•28130E+03 

•  26662  E  +  0  3 
• 25264E+03 

•  239  3  3  E  +  03 
.  22666  E  +  03 
•21462E+03 

•  203  19E  +  03 


•  678 1 5  E-04 
.  112  59  E  — 03 
•15807E-D3 

.33703 E-04 

•  8  59  4  0  E-04 
•22944 E- 03 

•  44376E— 03 


•11372 E-04 


51 

60 

•OOOOOE+OO 

•OOOOOE+OO 

•OOOOOE+OO 

•OOOOOE+OO 

• OOOOOE+OC 

•OOOOOE+OO 

.OOOOOE+OO 

.OOOOOE+OO 

VALUES 

^RUM 

J  =  61  - 

THROUGH  80 

ARE  • OOOOOOOE+OO 

81 

90 

•OOOOOE+OO 

•OOOOOE+OO 

•22743E-04 

•22743E-04 

•OOOOOE+OO 

•OOOOOE+OO 

.OOOOOE+OO 

. OOOOOE+OO 

VALUES 

:  FROM 

J  »  91 

THROUGH  110 

ARE  .OOOOOOOE+OO 

111 

120  - 

• 54619E— 05 

— . 20835E— 04 

— .42291E— 04 

— •  33347F— 04 

—  •  14430E— 04 

.OOOOOE+OO 

.OOOOOE+OO 

• OOOOOE+OO 

121 

130 

•OOOOOE+OO 

•OOOOOE+OO 

•OOOOOE+OO 

.32547E-05 

•71384E— 05 

« 811 13E— 05 

•  8  7994E  — 05 

•  943  8 3E— 0 5 

131 

140 

•11562 E-04 

. 12255E— 04 

•12950E-04 

.13646 E-04 

•14344E-04 

•  15043  E— 04 

•15744E— 04 

. 16446F-04 

141 

150 

•18566E-04 

.19277E— 0^ 

•19990E— 04 

•20705 E-04 

«  21422E— 04 

•  221 42E-04 

.  22865E-04 

•  235906-04 

151 

160 

•352T5E— 03 

•OOOOOE+OO 

.000006+00 

. OOOOOE+OO 

•OOOOOE+OO 

•OOOOOE+OO 

.OOOOOE+OO 

. OOOOOE+OO 

VALUES  FROM  J  *  161  THROUGH  500  ARF  •OOOOOCOE+OO 

501  502  •OOOOOE+OO  .69389E-17 


.223866-04 

•  00  0  OCE  +  0  0 

•  OOOOOE+OO 
.10178E-04 

•  17 1 5  IE— 04 

•  243 19E  — 04 
•OOOOOE+OO 


AREA  Til 

1 

10 

11 

20 

21 

30 

31 

40 

41 

50 

51 

60 

61 

*■>0 

71 

80 

81 

90 

91 

100 

101 

110 

111 

120 

121 

130 

131 

140 

141 

15  0 

151 

160 

161 

1^0 

171 

180 

181 

190 

191 

200 

201 

210 

211 

220 

221 

230 

231 

240 

241 

250 

251 

260 

261 

270 

271 

280 

281 

290 

291 

300 

301 

310 

311 

320 

321 

330 

331 

340 

341 

350 

351 

360 

361 

3T0 

371 

380 

381 

39  0 

391 

400 

401 

410 

411 

420 

421 

430 

431 

440 

441 

450 

451 

460 

461 

470 

471 

480 

481 

490 

491 

500 

501 

502 

E  *  TO  BE  Af 
•61757E-07 
•60235E— 07 
•58505E— 07 
• 85706E— 07 
• 83654E— 07 
• 81868E— 07 
• 13338E— 06 
• 13127 E— 06 
•12958E— 06 
•T8013E-06 
•17890E-06 
« 17824E-06 
•37989E— 07 
« 10978 E— 06 
• 28146E— 06 
•13179E— 05 
•13338E— 05 
.13539 E-05 
• 13783E— 05 
•14071 E— 05 
•14403E— 05 
•1478 2 E-05 
•15207E-05 
•15680E-05 
•16204E— 05 

•  167"79E— 05 
•1?407E— 05 
•18090E— 05 
•18831 F— 05 
•19631E— 05 
•20494E— 05 
• 21422E— 05 
• 22418E— 05 
• 23485E— 05 
• 24626E— 05 
•25 846 E-05 

•  271 48E— 05 
• 28536E— 05 
•30015E-05 
•315 89 E-05 
.33263 E-05 
•35043 E-05 
•36934E-05 

•  389 43 E— 05 

•  410*75E— 05 
•43337E— 05 
• 45737E— 05 
•48282E-05 

•  50981 E— 05 
.  53841E— 05 
•56873E— 05 


EA/XIX 
•61856E— 07 
•60054E— 07 
• 58342E— 07 
• 85488E— 07 
•83463E-07 

•  81703E— O'7 
•13315E-06 

•  13109E— 06 
•12943E-06 
•17998E— 06 
• 17881E— 06 
•16589E— 06 

•  37993E-07 

•  12329E— 06 

•  30320E— 06 
•13193E-C5 
•13356E— 05 
.13561 c— 05 
•13810E— 05 
•14102E— 05 
•14439E— 05 
•14822E— 05 
•15252E— 05 
•15731E— 05 
•16259E-05 
•16839E— 05 
•17472E— 05 
•18161E-05 
•18908E— 05 
•19714E— 05 
•20584E-05 
• 21518E— 05 
•22521E— 05 
• 23595E— 05 
• 24745E— 05 
. 25973E— 05 
•27283E-05 
• 28680E— 05 
• 30168E— 05 
• 31752E— 05 
•33436E-05 
• 35227E— 05 
«  37130E— 05 
. 39150E— 05 
•41295E-05 
•43571E— 05 
•45985E— 05 
•43545E-05 
• 51260E— 05 
•54137E-05 
•57185E-05 


• 61757E— 07 
•59874E— 07 
• 58181E— 07 
• 85274E— 07 
• 83275E— 07 
•81542E-07 
•13293E— 06 
• 13090E— 06 
• 12929E— 06 
•17984E-06 
•17872E— 06 
•11317E— 06 

•  37999E— 07 
• 13760E— 06 
•32578E-06 
• 13207E— 05 
• 13375E— 05 
•13584E— 05 
•1383TE— 05 
• 14134E— 05 
•14475E— 05 
•14863E— 05 
•15298E-05 
•15781E— 05 
•16315E— 05 
•16900E— 05 

•  17539E— 05 
• 18233E— 05 
•18986E-05 
•19799E-05 
• 20674E— 05 
• 21615E— 05 
•22625E-05 
.  2370"rE— 05 
• 24864E— 05 
. 26100E— 05 
•27419E-05 
•28825E-05 
• 30322E— 05 
• 31916E— 05 
• 33611E— 05 
• 35412E— 05 
•37326E-05 
• 39359E— 05 
• 41517E— 05 
• 43806E— 05 
•46234E— 05 
• 48810E— 05 
• 51540E— 05 
•54434E— 05 


• 61560F— 07 

•  59696 E-07 

•  58022  E— 07 
• 85062E— 07 
.83090 E-07 

•  81383  E— 07 
. 13271 E-06 
• 13072 E-06 
•18133E— 06 
• 17970E— 06 
• 17864E— 06 
•70 516 E-07 
. 38006c— 07 
•15272E-06 
• 34923E— 06 

•  13  222  E— 05 
•13394E— 05 
• 13608E— 05 
• 13865E— 05 
• 14166E— 05 
« 14512  E— 05 
•14904E— 05 
, 15344 E-05 
•15832E-05 
•16371E— 05 

•  16961 E— 05 
•17606E-05 
• 18 3 06 E— 05 
•19064 E-05 
• 19883E— 05 

•  20765  E— 05 
• 21713E— 05 
•22730E-05 
•23819E-05 

•  24984  E— 05 
• 26228E— 0* 
•27555 E-05 
•28970E-05 
• 30477E— 05 
• 32080E— 05 
•33786E— 05 
.35598 E-05 
• 37524E— 05 
• 39569E— 05 
• 41739F— 05 

•  44042  E— 05 
« 46485E— 05 

•  49076  E— 05 
•51 8 22 E-05 
•54733 E-05 


•61365E-07 

•59521E-07 

•  578  65E-07 
•84853E-07 
• 82908E— 07 
•81226E-07 
♦13249E-06 

•  13055E— 06 
•18114E-06 
.179  57E-06 
•17857E-06 
• 37986E— 07 
•45355E-07 

•  168  64E-C6 

•  373  54E-06 
•13237E-05 
•13413E-05 
.13631E-05 
•13893E-05 
•14198E-C5 
•14549F-05 
•14946E-C5 
.15390E-05 
•15884E— 05 
•16427E-05 

•  170  23E-05 
• l7t 73E— 05 
•18379E-05 
•19143E-05 
•19O69E-05 

•  208  5^8-05 
•21812E-C5 
•22836F-05 
•23932E-05 
•25105E-05 
«  26357E— 05 
•27693E-05 

•  29117E— 05 
•30633E-05 
• 32246E— 05 

•  33962E— 05 

•  35786E— 05 
•37723E-05 
• 39780E— 05 
«  419  64E— 05 
•44280E-05 
•46737E-05 
•49343E-05 
•52105E-C5 
•55033E-C5 


•61172E-07 
.  59347E— 07 
•57710E-07 

•  84646E— 07 
• 82728E— 07 
.81072E-07 

•  13  228  E— 06 
• 13037E— 06 
•18096E-06 
•17944E-06 
•17850E-06 

•  37984E— 07 

•  54 1 28E— 07 

•  18538E— 06 

•  39872E— 06 

•  13  2  53E -05 
• 13433E— 05 

•  13655E— 05 
.139 21 E-05 

•  14231E— 05 
•14587E-C5 

•  14988  E-05 

•  15  438E— 05 
• 15936E— 05 

•  16485E— 05 
• 17086E -05 
• 17741E— 05 

•  18  453E— 05 

•  19223E— 05 
« 20054E— 05 

•  20949E— 05 
•21911E-05 

•  229  42E— 05 

•  24046E— 05 
•25226E-05 

•  26  487E-05 

•  27  8  31E— 05 

•  29264E— 05 
• 30790E— 05 
•32413E-05 
• 34140E— 05 

•  35974  E— 05 
• 37923E— 05 
• 39993E— 05 
•42189E-05 
• 44519E— 05 
• 46991E— 05 
• 49612E-05 
•52390E-05 
•55335E-05 


•  60981E— 07 

•  5  9175E— 07 
•57556E— 07 
.R4442E-07 
• 82551E— 07 

•  8  0921 E-07 

•  13207E-06 
•130216-06 
•18078E— 06 

•  17932E-06 

•  178  44E-06 

•  37982E-07 
•63684E-07 

•  2  0294E-06 

•  4  2479E-06 

•  1 3  269E-05 
• 1 3453E— 05 
•13680E-05 

•  1 395  0E-05 

•  14  265E-05 
•14625E-05 

•  1 503  IE  — 05 

•  1 5485E— 05 

•  1 5988E— 05 
• 16542E— 05 
•17149E-05 
. 1 7  81 OE— 05 
•18527E-05 
.  19303E— 05 
•20141E-05 

•  2 1043E-05 
•22011E-05 
•23049E-05 

•  2  4 1 61 E  — 0  5 
•25349E-05 
•26617E-05 
•27970E-05 
.294126-05 
• 30948E— 05 
•32581E-05 

•  3  43 1 8E— 05 

•  3  6 164E-05 
. 38125E— 05 
• 40207E— 05 
•42416E-05 

•  4  4760E-05 
•47246E-05 

•  4  9883E-05 

•  5  2677E-05 

•  5  5  639E-05 


• 60791E— 07 
• 59004E-07 

•  574  05E-07 

•  842  41E-07 
• 82376E-07 
• 80772E— 07 
• 13186E-06 
• 13004E— 06 
. 18061E— 06 

•  179  21E-06 
• 17838E-06 
. 37982E-07 

•  74024E-0  7 

•  221 32E-06 

•  4  51 74E-0  6 
•13286E-05 
• 13474E— 05 
•13705E-05 
• 13980E-05 
• 14299E-05 

•  14663E-0  5 

•  15074E-05 

•  155  33E-0  5 
•16041E-05 

•  16601 E-05 
.172126-05 
•17879E-05 
• 18602E-05 
•19384E-05 
.202286-05 

•  2 11 36E-0  5 
.221126-05 
. 231 57E-05 
. 242766-05 
.254726-05 
•26749E— 05 

•  2  81 1 0E-05 

•  295  61 E— 0  5 

•  311066-05 
•327506-05 
. 344986-05 

•  363  55E-05 
.  383  27E-0  5 

•  404  22E-05 
•42644E-05 
• 45002E-05 

•  475  03E-05 

•  501 5  5E— 05 
. 52366E-05 

•  559456-05 


•60604E-07 
.  5 8 836E-07 
•572556-07 

•  84043E-07 
.822046-07 

•  80626E-07 

•  13166E-06 
•12989E-06 
• 18044E— 06 
• 17910E-06 

•  17  83  3E— 06 

•  37983E-07 
•85152E-07 

•  24  053E-06 
•479596-06 
•13303E-05 
•13495E-05 
•13731E-05 
•14010E-05 

•  14333E— 0  5 
• 14702E-05 
•15118E-05 
.155826-05 
.160956-05 

•  16659E-05 

•  172  77E-0  5 
.179496-05 

•  18678E-05 
.194666-05 

•  20316E— 05 

•  21231E— 05 
.222136-05 

•  23  265E-05 
.24392E-05 
•25596E-05 

•  26  8 8 IE  — 0 5 

•  28  2  51E— 05 
•29712E-05 
. 312666-05 
.  32920E-05 
.346786-05 
•36547E-05 
•38531E-05 
.406336-05 

•  42  8  74E  — 0  5 
.45246E-05 

•  47761E-05 
.50429E-05 
.53  2  56E-05 
•56252E-05 


•22386 E- 04 

.OOOOOE+OO 

.000006+00 
.  10869E-04 
•17857 E-04 
•  365 10E-03 
.000006+00 


•  60  4 1 9  E— 0  7 

•  5  8670E-07 
.85926  E-07 
.83847E-07 
•82034E— 07 
•133626-06 

•  13 1 47E-06 

•  1  2973E-36 

•  18028  E-06 

•  179  00E-06 

•  178  2  8  E-06 

•  37985 E-07 
, 97070E-07 

•  260  58  E-06 

•  508  3 5E-05 
. 13320E-05 
.  135 1 7E-05 

•  13757E-05 
•14040 E-05 
. 14363 E-05 
.147426-05 

•  1 5 1 6  2  E— 0  5 
. 15631E-05 
• 16149E-05 

•  1 671 9  E-0  5 
• 17341E-05 
. 18  019  E— 0  5 

•  18754E— 05 
.19548 E-0 5 

•  2040  5  E-0  5 

•  21326E-05 

•  223 1 5  E-0  5 
.  233  7  5  E— 0  5 

•  245  09E-05 
•25721E-05 

•  2  70 1 4E-0  5 
•28393E-05 

•  29863  E-05 
.  31427E-05 
•33091E-05 
.  34  8  60  E-0  5 
• 36740E-05 
.  38736  E-0  3 
. 40856 E-05 
•43105 E-05 
.45491 E-05 

•  48021 E-05 
. 50704E-05 
. 53543E-05 
. 56562E-05 


OEMS TTY 

VALUES  FROM  J  *  1  THROUGH  110  ARE  .1000000E+01 

HI  120  .10000E+01  •10000E+01  .10000E+01  .10000E+01 

VALUED  FROM  J  *  121  THROUGH  502  ARE  .1477737E-01 

VELOCITY 

VALUES  FROM  U  *  1  THROUGH  502  ARE  .OOOOOOOE+OO 

PRESSURE 

VALUES  FROM  J  *  1  THROUGH  110  ARE  .1137613E+01 

111  120  .11376E+01  •11376E+01  •11376E+01  .11376E+01 

VALUES  FROM  J  =  121  THROUGH  502  ARE  •1055300E-01 

SOTMD  S’PEED 

VALUES  FROM  J  ~  1  THROUGH  110  AP.E  •1262141E+01 

111  120  .12621E+01  •12621E+01  .12621E+01  ,12621E+01 

VALUES  FROM  J  =*  121  THROUGH  502  ARE  .1000000E+01 


o 


* 


t 


100C0E+01  .10000E+01  .10000E+01  •10000E+01  .14777E-01  .14777E-01 


11376E+01  •  11376E  +01  .11376E  +  01  .11376E+DI  ,  105536-01  .10553E-01 


.  12621E+01  .12G21E+01  .12621E  +  01  •12621E  +  01  .100006+01  ,10OOOEO1 


**  n  *  4200  AND  TIME (SEC ) 


36513E-01 


DENSITY 


1 

10 

11 

20 

21 

30 

31 

40 

41 

50 

51 

60 

61 

70 

71 

80 

81 

90 

91 

100 

101 

110 

111 

120 

121 

130 

131 

140 

141 

150 

151 

160 

161 

170 

171 

180 

181 

190 

191 

200 

201 

210 

211 

220 

221 

230 

231 

240 

241 

250 

251 

260 

261 

2^0 

2TL 

280 

281 

290 

291 

300 

301 

310 

311 

320 

321 

330 

331 

340 

341 

350 

351 

360 

361 

370 

3-n 

380 

381 

390 

391 

400 

401 

410 

411 

420 

421 

430 

431 

440 

441 

450 

451 

460 

461 

470 

471 

480 

481 

490 

491 

500 

501 

502 

• 57873E— 01 
•57929 E— 01 
• 5 8546 E— 01 
.60413E-01 
•61456E— 01 
*  629 77 E— 01 
•65984E— 01 
•66933E— 01 
•66884E— 01 
• 66854E— 01 
• 668 3 3 E— 01 
•66639E— 01 
• 37153E— 01 
.81071 E-02 
•16304E— 01 
.17945E-01 
• 17462E— 01 
•16756E— 01 
• 1*5690 E— 01 
•14680E— 01 
•14392E-01 
•14231E-01 
•13932E-01 
•13307E— 01 
•12991E— 01 
« 11953E— 01 
•92462E— 02 
• 83940E— 02 
* 86251E— 02 
.88314E-02 
•90641 E— 02 
•92933 E-02 
.94738 E-02 
.95671 E-02 
• 9 62 57 E— 02 
•96401E— 02 
.96541 E— 02 
•98234E— 02 
•10019E— 01 
•11302E-01 
•14352E— 01 
•15335E— 01 
« 159  5  2  E— 01 
•15917E— 01 
•15811E— 01 
• 15410E— 01 
•15600E-01 
•16202E-01 
• 16979E— 01 
•16725E— 01 
•14910E— 01 


•57867E— 01 
•57950 E-01 
•53675E— 01 
• 60357E— 01 
•61586E-01 
•63192E— 01 
•65972E-01 
•66955E— 01 
•66972E— 01 
•66852E— 01 
•66822E-01 
• 66885E— 01 
•38201E-01 
•71036E— 02 
•16479E— 01 
•18267E-01 
• 17341E— 01 
• 16712E— 01 
• 15563E— 01 
•14620E— 01 
•14382E— 01 
•14210E-01 
• 13883E— 01 
•13245E— 01 
•12974E— 01 
•11707E-01 
•90466E— 02 
•84065E— 02 
• 86467E— 02 
•38532E— 02 
•90882E— 02 
•93141E— 02 
•94875E-02 
•95726E— 02 
• 96315E— 02 
•96369E— 02 
•96654E— 02 
•98429E— 02 
•10038E-01 
•12109E-01 
•14503E— 01 
•15412E-01 
•15976E— 01 
•15906E-0I 
•15785E— 01 
•15379E— 01 
• 15663E— 01 
•16272E— 01 
• 17044E— 01 
•16554E— 01 
• 14715E— 01 


• 57873E— 01 
• 57976E— 01 
• 58816E— 01 
•60450E— 01 
•61719E-01 
•63429E— 01 
« 66135E— 01 
•66966E— 01 
• 66759E— 01 
•66851E-01 
• 66811E— 01 
• 66390E— 01 
• 38116E— 01 
• 63172E— 02 
• 16830E— 01 
•18053E— 01 
•17217E— 01 
•16650E-01 
• 15441E— 01 
•14568E— 01 
• 14370E— 01 
•14188E— 01 
• 13830E— 01 
•13190E— 01 
•12948E— 01 
• 11441E— 01 
• 88751E— 02 
• 84255E— 02 
•86676E— 02 
• 88755E— 02 
• 91121E— 02 
• 93344E— 02 
•95002E— 02 
• 95780E— 02 
• 96366E— 02 
•96338E— 02 
• 96786E— 02 
• 98624E— 02 
• 10054E— 01 
• 12851E— 01 
•14632E-01 
• 15488E— 01 
•15990E-01 
• 15896E— 01 
• 15755E— 01 
•15356E-01 
• 15726E— 01 
.  16345E— 01 
•17097E— 01 
• 16361E— 01 


• 57872E— 01 
.58010 E-01 
.58973 E-01 
•60600 E-01 

•  61 8 54 E— 01 
• 63695E— 01 
• 66328E— 01 
.66968 E-01 
.66 819 E-01 
•66850E-01 
• 66803E— 01 
• 64700E— 01 
.33350 E-01 
• 54124E— 02 
.17051 E-01 
• 17999E— 01 
•17100E— 01 

•  16568  E— 01 
•15324E— 01 
• 14526E— 01 
• 14 3 57 E— 01 
.14165 E-01 
•13772E— 01 

•  131 42  E— 01 
•12909E— 01 
• 11158E— 01 

•  87323  E— 02 
•84488E— 02 
. 86880E— 02 
• 88981E— 02 
.91359  E— 02 
,93 541 E-02 
. 95118  E— 02 

•  95834  E— 02 
.96408 E-02 
• 96312E— 02 

•  96935  E— 02 
• 98817E— 02 
.10064 E-01 
• 13369E— 01 
•14741 E— 01 
. 15562E— 01 
•15995E— 01 
.15888  E— 01 
• 15718E— 01 
• 15345E— 01 
• 15786E— 01 
.16421 E-01 
.17135 E-01 
.16150 E-01 


.57875 E-01 
• 58051E— 01 
•59143E-01 
. 60719E— 01 
•61991E— 01 

•  63996E— 01 
• 66467E— 01 
• 66963E— 01 
• 66977E— 01 
.  66849E— 01 
•66799E— 01 
. 571 07E— Cl 

•  256  35E-01 
• 53992E— 02 

•  172  22E— 01 
• 17967E— 01 
. 17001E-01 
. 16468  E— 01 
* 15212E-01 
• 14492E— 01 

•  143  42E— 01 
. 141 40E— 01 
• 13711E— 01 
.131 01E-01 

•  128  51E— 01 
.10866E-01 

•  861 75E— 02 

•  847  44E— 02 

•  870  81E— 02 

•  89212E-0  2 

•  915  94E— 02 
• 93733E-02 

•  952  24E— 0  2 

•  958  89E— 02 
.96440E-02 
• 96294E— 02 
• 97098E— 02 

•  990  09E— 02 
•10065E-01 
. 13642E— 01 
. 14836E— 01 
• 15634E-01 
•15993E— 01 
.15881E— 01 
• 15678E— 01 
• 15347E-01 
« 15845E— 01 

•  165  00E— 01 
• 171 55E— 01 
. 159  31E— 01 


•57879E— 01 

•  58 103E-01 

•  59291E— 01 
.608346-01 
•62133E— 01 
.64  290E-01 

•  665  89E— 01 
.669546-01 
• 66866E-01 
. 66849  E— 01 

•  66788  E-01 
•44108E— 01 
.19807E-01 

•  40  886E— 02 
• 17408E— 01 

•  17  898  E— 01 
• 16923E— 01 
•16352E-01 

•  15107E— 01 

•  14  4  64E— 01 

•  143  26E— 01 

•  14 1 14  E— 01 
. 13646E— 01 
.130686-01 
• 12770E-01 
.105696-01 
• 85293E— 02 
.850116-02 

•  872  80E— 02 

•  89445  E— 02 

•  918  27E— 02 
.93919E-02 
. 95320E-02 
• 95947E-02 

•  96461 E— 02 
. 962  89E— 02 

•  97  274E-02 

•  992  01E— 02 
.  10055E— 01 
.  13756E— 01 

•  14924E— 01 
.  157  03E— 01 
.159866-01 
• 15875E-01 
.156346-01 
.153626-01 

•  15902E— 01 
.165816-01 
.171526-01 
.157156-01 


.578846-01 

.581656-01 

•  59486E-01 
.609566-01 
•62280E— 01 

•  6  4630E— 01 
.666956-01 
.669426-01 
.66856E— 01 
.66849E— 01 

•  6  6747E-01 

•  3  8  46  5E-01 
• 1 5830E— 01 

•  32100E— 02 
.175726-01 
.17828E— 01 
• 16868E-01 
.1 6225E— 01 

•  1 5007E-01 
• 1 4443E-01 
• 1 4309E— 01 
. 14084E— 01 

•  1 3578E-01 

•  1 3044E— 01 

•  1 2  66  4E— 01 
. 1 0276E— 01 
. 84655E— 02 
.852776-02 
.874816-02 

•  8  9681 E  — 02 
.92056E-02 
•94098E-02 
. 95405E— 02 
. 96007E— 02 

•  96470E-02 
.962996-02 

•  97458E-02 

•  99394E— 02 
,  10045E-01 
•1 3818E— 01 
.15009E-01 
• 15768E-01 
.15974E— 01 
.158686-01 
.155876-01 
.153916-01 

•  1 5  95  8E— 01 
.166646-01 
• 17125E— 01 
•15513E— 01 


•  5789  IE-01 

•  5  82  40E-01 
.598496-01 
.610786-01 
.624  3  5E-01 
.651956-01 
.667816-01 
.669316-01 
.668646-01 
. 668  48E-01 
. 66660E-01 
• 36656E— 01 

•  1 3031E-01 

•  88672E-02 
.177156-01 
. 177556-01 
.168326-01 

•  1609  2E-01 
. 1491 3E-01 
.144276-01 
.142916-01 
• 14052E-01 

•  135  09E— 01 
.130266-01 
• 12530E-01 
.999176-02 

•  842  31E-Q2 
.855376-02 
,  876  84E-02 
• 89920E— 02 
.922826-02 
.942706-02 
.954826-02 
.96069E-02 
,964676-02 

•  963  29E— 02 
• 97648E-02 

•  995  89E-02 
.100786-01 
. 13902E-01 
.150926-01 
.  158  26E-01 
.159606-01 

•  158  59E-01 
.155396-01 
,154316-01 
.160156-01 
.167466-01 
. 17069E-01 
.153406-01 


.579016-01 

.583286-01 

•  59645E-01 
.612016-01 

•  62  601E-01 
.650716-01 
.  66848E— 01 
. 66922E— 01 

•  66  8  58E-01 
,  668466-01 

•  66572E-01 

•  361 40E-01 

•  10956E— 01 
.153056-01 
.180926-01 

•  17  6  7 IE  — 01 
.168086-01 
.159576-01 
.148276-01 
.144146-01 
•14271E-01 
.140166-01 

•  13  440E— 01 
.130136-01 
.123676-01 
.97221E-02 
.839926-02 
.857876-02 
.878906-02 
.901596-02 
.925046-02 

•  94435E  — 02 
.955516-02 
.961336-02 
.964536-02 
,963786-02 
.978426-02 
.997866-02 
.102416-01 
.140296-01 
.151756-01 

•  15  877E— 01 
. 15945E-01 
. 15847E-01 
.154936-01 
.154816-01 

•  16075E— 01 

•  16  8  28E-01 

•  16984E  — 01 
• 15196E-01 


.  57913  E-01 
,58430  E-0 1 

•  599  3  8  E-01 
.613  27  E-0 1 
.627816-01 

•  6  5400  E-0 1 
.66398  E-0 1 

•  6  68  98  E— 0 1 

•  668  5  5  E-0 l 
.668  4  2  E-0 1 

•  665  5  4E— 01 
. 362536-01 
.937506-02 
.16670  E-0 1 
.  1768  3  E-0 1 
,17574  E-01 
.167856-01 

•  1582  2  E-01 
.  147  4  9  E— 0 1 

•  14403  E-01 
.142516-01 
,139766-01 

•  133  72  E-01 
.130026-01 
.121746-01 

•  947  23  E— 02 
.  93906 E-02 
.860246-02 

38100E-02 
.904006-02 
, 92721 E-02 
.945916-02 

•  956 1 3  E-02 
.9619  6  E-02 
. 964  30  E-0 2 
.  964  4  9  E— 0  2 
.98038E-02 
• 99986E-02 
.10640 E-01 
.1418  8  E-0 1 

•  15256E-01 

•  15919 E-01 

•  15931 E-01 
.15831 E-01 

•  15  449  E— 01 

•  1 5  5  3  8  E-0 1 
• 161 37E-01 
.  169  06  E-01 
. 16869E-01 

•  1 5 06 5  E-0 1 


VELOCITY 
1  10 
11  20 
21  30 

31  40 

41  50 

51  60 


30512E— 02 
27037E-01 
51804 E— 01 
44410E-01 
64313E-01 
71953E-01 


OOOOOE+OO 
29917E— 01 
53654E— 01 
48232E— 01 
65798E-01 
71356E-01 


30512E-02 
32742E—01 
55394E— 01 
50445E— 01 
67178E-01 
70381E— 01 


.61341 E-02 
• 35496E— 01 
• 57100E— 01 
• 51984E— 01 
•  68437  E— 01 
.69066 E-01 


.91558E-02 
•381666-01 
• 58775E— 01 
.538786-01 
•  695  57E— 01 
•67460E— 01 


.121676-01 
• 40738E— 01 
.600356-01 
• 55779E— 01 
•  70516E— 01 
.65106E-01 


•  15178  E— 01 
.432006-01 
• 6 1787E-01 
.575896-01 
.712906-01 
.  62801E-01 


.  18176E-01 
.455406-01 
•  653  63E— 01 
.5936 1 E— 0 1 
.718  5 2E-0 1 
.626426-01 


.211 56E— 0 1 
•  477  54E-01 
.603586-01 
. 61 0  84E— 01 
.721726-01 
.531896-01 


•241 1 3E-01 

•  498  39  E-3 1 
,48248 E-01 

•  627  3  8E-0 1 
.72215  E— 01 
.365  9  3  E-0 1 


61 

70 

71 

80 

81 

90 

91 

100 

101 

110 

111 

120 

121 

130 

131 

140 

141 

150 

151 

160 

161 

170 

171 

180 

181 

190 

191 

200 

201 

210 

211 

220 

221 

230 

231 

240 

241 

250 

251 

260 

261 

270 

271 

28  0 

281 

290 

291 

300 

301 

310 

311 

320 

321 

330 

331 

340 

341 

3*>0 

351 

360 

361 

370 

3"?1 

380 

381 

390 

391 

400 

401 

410 

411 

420 

421 

430 

431 

440 

441 

450 

4*51 

460 

461 

470 

471 

480 

481 

490 

491 

500 

501 

502 

PRESSURE 


1 

10 

11 

20 

21 

30 

31 

40 

41 

50 

51 

60 

61 

^0 

71 

80 

81 

90 

91 

100 

101 

110 

111 

120 

121 

130 

131 

140 

141 

150 

151 

160 

161 

170 

•28545E-01 
•22554E-01 
•34692E-01 
•41629E— 01 
. 62172 E-01 
• 83852E— 01 
•77965E+00 
•12577E+01 
• 23906E+00 
•62019E— 01 

•  64653E— 01 
•79324E-01 
•94030E— 01 
•10398E+00 
•10973E+00 
•11844E+00 
•13014E+00 
•14315E+00 
• 15696E+00 
•17093E+00 
•18483E+00 
•19930E+00 
•214S8E+00 
•23184E+00 
•25038E+00 

•  27057E+00 
•29204E+00 

•  31386E+00 
•33530E+00 
.35674E+00 
• 37876E+00 
•40118E+00 
•42278E+00 
•34562 E+00 
. 15204E+00 
•11648E+00 
•92162E— 01 
•99211E— 01 
•15031E+00 
•20534E+00 
• 21318E+00 
• 20356E+00 
• 19079E+00 
•18231E+00 
•20076E+00 


•21098E— 01 
•21128E-01 
•21444E-01 
•22419E— 01 
• 22972E— 01 
•23768E-01 
•25409E-01 
•25910E-01 
•25896E-01 
• 25883E— 01 
• 25872E— 01 
•25699E-01 
•11152E-01 
•13173E-02 
•63020E-02 
•67822E— 02 
•68895E-02 


•  28  546E— 01 
• 23265E— 01 
•37381E-01 
•43614E-01 
•64315E-01 
•9 44 09 E—01 
•75875E+00 
•12835E+01 

•  22433E+00 
•45425E-01 
•66081E-01 
• 80839E— 01 
•95332E— 01 
•10461E+00 
•11041E+00 
•11951E+00 
•13140E+00 
•14450E+00 
• 15836E+00 
•17232E+00 
• 18624E+00 
• 20080E+00 
«  21651E+00 
• 23362E+00 
•25233E+00 
• 27267E+00 
•29422E+00 
•31603E+00 

•  33743E+  00 
•35891E+00 
•38100E+00 
•40341E+00 
•42493E+00 
• 28232E+00 
•14534E+00 
•11396E+00 
.90470E-01 
•10294E+00 
•15632E+00 
•20913E+00 
•21167E+00 

•  20252E+00 
•18936E+00 
•18358E+00 
• 20063E+00 


•21095E-01 
•21138E-01 
•21511E-01 
• 22390E— 01 
• 23039E— 01 
•23883E-01 
•25404E-01 
•25921E-01 
• 25947E— 01 
•25880E-01 
• 25868E— 01 
•25566E-01 
•11595E— 01 
•10986E— 02 
•63389E— 02 
•69562E-02 
•  68927E— 02 


•26804E-01 
•24143E-01 
• 35927E— 01 
•45625E-01 
•66481E— 01 
• 13745E+00 
•77204E+00 
•12985E+01 
•20329E+00 
•53813E-01 
• 67515E— 01 
•82356E-01 
•96576E-01 
•10521E+00 
•11113E+00 
•12062E+00 
•13266E+00 
•14586E+00 
•15977E+00 
•17370E+00 

•  18765E+00 
• 20231E+00 
•21816E+00 
• 23542E+00 
•25429E+00 
•27479E+00 
• 29641E+00 
•31818E+00 

•  33956E+00 
•36109E+00 
• 38324E+00 
«  40563E+00 
•42731E+00 
• 22841E+00 

•  14010E  +  00 
•11148E+00 
• 89196E— 01 
•10711E+00 
•16233E+00 
•21230E+00 
• 21029E+00 
•20139E+00 
•18792E+00 

•  18545E+00 


• 21098E— 01 

•  21152E— 01 
• 21583E— 01 
• 22439E— 01 
• 231 07 E—01 
•24010E-01 
•25494E— 01 

•  25927E— 01 
•25835E-01 
• 25878E— 01 
• 25862E— 01 
• 25111E— 01 
•11530E-01 
•91883E— 03 
•64718E— 02 
•68584E— 02 
• 68959E— 02 


•  24858  E—01 
.25160 E-01 
.303 12 E-01 
• 47656E— 01 
• 68663E— 01 
•25189E+00 
• 85018E+00 

•  13338  E+01 
•18892 E+00 

•  55517 E—01 
• 68957E— 01 
• 83874E— 01 
• 97754E— 01 
•10576E+00 
•11189E+00 
•12174E+00 
« 13394E+00 
• 14723E+00 
•16117E+00 
• 17509E+00 

•  18908  E+00 
• 20384E+00 
•21981E+00 
• 23723E+00 
•25627E+00 
•27691E+00 
• 29860E+00 
•32034E+00 

•  34169  E+00 

•  36327  E+00 
• 38548E+00 
• 40784E+00 
• 43020E+00 
• 19412E+00 
•13612E+00 
•10900E+00 

•  88413  E—01 
• 11166E+00 
• 16831 E+00 
•21476E+00 
• 20912E+Q0 
•20017E+00 
•18649E+00 
•18782E+00 


•  21098 E-01 

•  21169  E—01 

•  21664 E—01 
.22518 E-01 
•23 176 E-01 
• 24151E— 01 
• 25600E— 01 
• 25929E— 01 
•25871 E-01 
.25876 E-01 

•  25855  E—01 
.24381 E-01 

•  95435  E— 02 
• 77773E— 03 
« 65  383  E— 02 

•  68  565  E-02 

•  68991 E-02 


•23618E-01 

•26293E-01 

•  294  26E-01 
•49701E-01 
•70854E-01 
•46017E+00 

•  951 32E+00 
• 131 86E+01 
•17719E+00 

•  558  82 E—01 
•70407E-01 
•85386E-01 
•98862E-01 

•  106  31E+00 
•11271E+00 

•  122  89E+00 

•  135  23E+00 

•  148  60E+00 

•  162  57E+00 
• 17647E+00 

•  190  51E  +  00 

•  205  38E  +  00 
•22149E+00 
• 23906E+00 

•  258  27E  +  00 
•27905E+00 
• 30079E+00 
• 32249E+00 

•  343  83E+00 

•  365  47E  +  00 

•  387  72E+00 
• 41004E+00 
• 43396E+00 
. 178  60E+00 

•  132  93E  +  00 

•  106  50E+00 

•  881 79E-01 

•  116  54E+00 

•  174  23E+00 

•  216  46E+00 
•20815E+00 
•198918^00 

•  185 11E  +  00 
• 19055E+00 


• 21099E— 01 
• 21190E— 01 
• 21750E-01 
•22581E-01 
• 23247E— 01 
•24309E— 01 
• 25675E— 01 
•25929E-01 
•25959E— 01 

*  25875E— 01 
•25846E-01 
• 20640E— 01 
•66283E-02 
•66513E— 03 

•  65817E-02 
•68731E— 02 
•69025E-02 


•22689E-01 

•27523E-01 

•32245E-01 

•  51755E— 01 
•73055E— 01 

•  65069E  +00 

•  10393E  +  01 

•  13377E  +01 
. 16659E  +00 
.  57479E— 01 
•71868E-01 
•86889E-01 
•99897E— 01 
•10684E+00 

•  11356E  +00 
•124 06 E+00 

•  13652E  +00 

•  14999E+00 
• 16397E+00 

•  17786E  +00 
•19195E+00 

•  20693E  +00 
•22318E+00 

•  24090E  +00 

•  26028E  +00 

•  28119E  +00 
• 30297E+00 
• 32463E+00 

•  34597E+00 

•  36767E  +00 

•  38997E  +00 
•41223E+00 

•  43866E+00 

•  17457E  +00 

•  13008E  +00 

•  10397E  +00 
•88533E— 01 

•  12170E  +00 
•18003E+00 

•  21738E  +00 

•  20735E  +00 
•19761E+00 
« 18384E  +00 

•  19338E  +00 


•  21102  E—01 
• 21217E— 01 
•21824E-01 

•  22643E— 01 

•  23321E— 01 

•  24460E— 01 

•  25739E— 01 
• 25927E— 01 
.258998-01 
•25874E-01 
♦25837E— 01 

•  14370E— 01 
•46186E— 02 
• 69436E— 03 
«  66323E— 02 
•68739E— 02 

•  69059E— 02 


•22054E-01 
. 28  83  LE— 01 
• 34104E— 01 
•53814E-01 
•75245E— 01 
• 74917E+00 
•11063E+01 
•1 691 8E+01 
.1 561 IE  +  00 

•  5  899  3E  —  01 

•  7  334  0E— 01 

•  8  8376E— 01 
•10086E+00 

•  10737E  +  00 

•  1 1 446E  +  00 
•12525E+00 
•13783E+00 
•1 5137E+00 

•  16537E  +  00 
•17924E+00 
•1934 OE +00 
•20849E+00 
•22488E+00 

•  24277E  +  00 
•26231E+00 

•  28335E  +  00 
•30516E+00 

•  3  2677E  +  00 
•34811E+00 

•  36987E  +  00 
•39222E+00 

•  41 43  9E  +  00 
•44325E+00 

•  1 741 IE  +  00 
•12730E+00 

•  10142E  +  00 

•  8  9496E— 01 

•  12711E  +  00 

•  18567E  +  00 
•21755E+00 

•  20665E  +00 

•  1 9628  E  +  00 
•18277E+00 
• 1 9608E+00 


•21104E— 01 
•21249E-01 
•21924E-01 

•  2  2709E— 01 
•23398E— 01 

•  2  4640E-01 
• 25794E-01 
•25923E— 01 

•  25893E-01 
•25875E— 01 
•25826E-01 
•11751E— 01 
•33697E-02 

•  0  0  00  OE  +  00 

•  6666 IE  — 02 
• 68760E— 02 

•  6  9  09  5E—  02 


•  217  596  —  01 
. 30226E-01 
•35827E-01 

•  5  58  80E— 01 
• 77401E— 01 
• 78357E+00 
• 11574E+01 

•  685  2  7E  +  00 
•14619E+00 
•60372E-01 
•74822E-01 
.  898 41 E  —  01 
.  10174E  +  00 
• 1 0792E+00 

•  115  40E  +  00 
• 12645E+00 
•13915E+00 
. 15276E+00 
•16677E+00 
. 18063E+00 
• 19486E+00 
• 21007E+00 
•22660E+00 
• 24465E+00 
• 26435E+00 
. 28551E+00 
• 30734E+00 

•  3  2891 E  +  00 

•  3  50  26E  +  00 
• 37209E+00 
• 39446E+00 
. 41652E+00 
•44389E+00 

•  172  08  E  +  00 
. 12452E+00 

•  988  9  5E  — 01 
• 91067E-01 
• 13271E+00 

•  191 09E  +  00 

•  21706E  +00 
• 20599E+00 
• 19494E+00 
• 18198E+00 
•19839E+00 


•21108E-01 
. 21287E-01 
•22115E-01 
• 22774E— 01 

•  234  80E-01 
•24950E-01 
• 25837E-01 
• 25919E-01 
•25895E-01 

•  258  75E-01 

•  258 14E-01 
« 10947E— 01 
• 25635E-02 

•  2  79  0  OE  — 0  2 
•66868E-02 
• 68799E-02 
•69132E-02 


.21767E-01 
. 31721E-01 
.  37  7  49E  — 0 1 
•57957E-01 
.797076-01 
.794366+00 
•11 974E+01 
• 28378E+00 
•1476 4E +00 
.617956-01 
.763146-01 
.912776-01 
.102556+00 
.108496+00 
.116386+00 

•  12767E  +00 
.  140476+00 
.154166+00 

•  16  8 16E  +00 
•182036+00 

•  1 96  3  3E  +0  0 
•21166E+00 

•  22833E  +  00 
•246546+00 
•26641E+00 

•  287686  +  00 

•  309  52E  +  00 

•  33 1 04E  +  00 
• 35241E+00 

•  37  4  31E  +  00 
•39670E+00 
•41863E+00 
•432586+00 
« 16697E  +00 
•12177E+00 
•96460E-01 

•  93  2  31E-01 
•13848E+00 

•  196  23E  +  00 
• 21605E+00 
•20527E+00 
.193  58E  +00 
.181 57E  +00 

•  20007E  +00 


•2111 3E-01 
•21332E-01 

•  22  01 7E— 01 

•  22  8  39E— 01 
.23  5  68E— 01 
•24903E-01 

•  25  870E— 01 
.  25  9 1 5E-0 1 
.25  8  90E— 01 

•  25  875E-0 1 
.25  8  08E— 01 

•  1 07  25E-01 
• 20099E— 02 
•60037E-02 
.68506E-02 
.688  3  IE  — 0  2 
.69170E— 02 


.22043  E—01 
.33111 6-0 l 
•396826-31 

•  600  5  3  6-31 
.819176-01 

•  793  23  6  *  >J 
.  123036  +  01 
.228226+30 
. 11430E+00 
.63227  E-01 
.77815E-01 
.926766-01 
.  103306  +  00 
.  109096  +  00 
. 117  39  E  +  00 
.12890E+00 
.141816+00 

•  1 55  5  6  E  +  00 
.169556+00 
.18342 E+00 
•19781E+00 
.  21326E  +  00 

•  23008  E  +  00 
.248456+00 
.268  48  E  +  00 

•  28  9  8  6  E  +  00 
. 3 11 69E+00 
.  3  33 17E  +  00 
. 35  45  7E  +  00 
. 37653 E+30 
.  3 98 94E+00 

•  42071 E  +  00 
•40G11E+00 

•  1 59  70E  +  00 
.119  08  E  +  00 

•  94 1 8  7  E-0 1 
.959586-01 
.  144  35  E  +  00 
.20101 E+00 
.  21 469E  +  00 
.20447E+00 
.  192  20E  +  00 

•  18166E  +  00 

•  200  3  5  E  +  00 


.211 2  0E-Q1 

•  213  84  E-0 1 

•  2217  3  E-0 1 

•  229  0  5  E-01 

•  23663  E-01 

•  25093E-0  L 

•  25894  E— 0 1 

•  2  59  03  E-0 1 
.25  8  86  E-0 1 
• 25874E-01 
• 25790E-01 

•  1 07  69  E— 0 1 
.161556-02 

•  6  5  7  7  2  E-0  2 
. 6  64 1 4  E-0  2 
•68863 E-02 
.  6  92  09E-02 


171 

180 

181 

190 

191 

200 

201 

210 

211 

220 

221 

230 

231 

240 

241 

250 

251 

260 

261 

270 

271 

280 

281 

290 

291 

300 

301 

310 

311 

320 

321 

330 

331 

340 

341 

350 

351 

360 

361 

3^0 

371 

380 

381 

390 

391 

400 

401 

410 

411 

420 

421 

430 

431 

440 

441 

450 

451 

460 

461 

470 

471 

480 

481 

490 

491 

500 

501 

502 

SOUND  SPEE 

1 

10 

11 

20 

21 

30 

31 

40 

41 

50 

51 

60 

61 

70 

71 

80 

81 

90 

91 

100 

101 

110 

111 

120 

121 

130 

131 

140 

141 

150 

151 

160 

161 

170 

171 

130 

181 

190 

191 

200 

201 

210 

211 

220 

221 

230 

231 

240 

241 

250 

251 

260 

261 

270 

271 

280 

.692 49 E-02 
•69587E— 02 
.69321 E-02 
•68540E— 02 
.68068E— 02 
• 67885E— 02 
.67782E-02 
•67695 E— 02 
•67551 E— 02 
•67357E— 02 
• 671 86 E— 02 
• 67049 E— 02 

•  6  69  5  5  E— 02 
•66917E-02 
.66933 E-02 
•66961E-02 

•  66904E— 02 
•66681E-02 
•66316E— 02 
.65B54E-02 
• 65257E— 02 
.64481 E-02 
•73463E—02 
•99212E— 02 
• 10713 E— 01 
•11080E— 01 
•10826E-01 
•10122E-01 
•93481 E— 02 
•93276E— 02 
.96685E-02 
•10068E— 01 
•10425E— 01 
.105 05 E-01 

D 

•71449E+00 
•71464E+00 
•71617E+00 
•72088E+00 
•72348E+00 
•72696E+00 
•73432E+00 
•73625E+00 
•73633E+00 
•73630E+00 
•73626E+00 
•73486E+00 
•64833E+00 
•47700E+00 
•73570E+00 
•72750E+00 
•74328E+00 
•76073E+00 
•78806E+00 
•81317E+00 
• 81662E+00 
•81839E+00 
• 82603E+00 
•84457E+00 
•85423E+00 
•88959E+00 
•10100E+01 
•10587E+01 


•69289E-02 
•69598E— 02 
• 69251E— 02 
•68473E— 02 
•68041E-02 
•67873 E-02 
.67774 E-02 
•67684E— 02 
•67533E— 02 
•67339E— 02 
•67171E-02 
•67038E— 02 
•66949E— 02 
• 66916E— 02 
• 66937E— 02 
• 66961E— 02 
•66890E— 02 
•66649E— 02 
•66274E— 02 
•65802E— 02 
•65187E— 02 
•64388E— 02 
•80689E— 02 
•10042E— 01 
•10772E— 01 
•11083E— 01 
• 10772E— 01 
• 10039E— 01 
•92974E— 02 
•93664E— 02 
•97043E— 02 
•10111E— 01 
•10434E-01 
•10553E— 01 


•71447E+00 
.71470E+00 
•71649E+00 
•72074E+00 
•72376E+00 
•72748E+00 
•73432E+00 
•73628E+00 
•73656E+00 
•73627E+00 
•73626E+00 
•73161E+00 
•65194E+00 
• 46536E+00 
•73393E+00 
•73023E+00 
•74606E+00 
•76194E+00 
•79134E+00 

•  81443E+00 
•81652E+00 
•81884E+00 
•82740E+00 

•  84646E+00 

•  85470E+00 
•89875E+00 
•10209E+01 
•10578E+01 


• 69330E— 02 
«  69602E— 02 
• 69175E— 02 
«  68409E— 02 
• 68017E— 02 
• 67861E— 02 
•67765E— 02 
• 67673E— 02 
•67513E— 02 
• 67320E— 02 
• 67156E— 02 
• 67026E— 02 
•66943E— 02 
• 66916E— 02 
• 66940E— 02 
• 66960E— 02 
•66873E— 02 
• 66617E— 02 
• 66231E— 02 
•65749E-02 
• 65114E— 02 
• 64271E— 02 
•87372E-02 
•10144E-01 
•10828E— 01 
• 11078E— 01 
•10713E— 01 
•99542E-02 
•92566E-02 
.94047E-02 
• 97415E— 02 
• 10154E— 01 
•10436E— 01 


•7144QE+00 
•71476E+00 
• 71684E+00 
•72096E+00 
•72405E+00 
•72805E+00 
•73471E+00 
• 73630E+00 
•73614E+00 
•  .73624E+00 
•73624E+00 
• 72777E+00 
•65083E+00 
• 45130E+00 
•73382E+00 
•72938E+00 
•74891E+00 
•76359E+00 
•79449E+00 
• 81542E+00 
•81647E+00 
•81933E+00 
•82892E+00 
• 84818E+00 
«  85548E  +  00 
•90903E+00 
•10306E+01 
•10565E+01 


.69371 E-02 

•  69598  E— 02 
• 69096E— 02 
• 68351E— 02 
•67996 E-02 
•67850E-02 
.67757 E-02 

•  67  661E-02 
• 67494E— 02 
•67303 E-02 
.67141 E-02 
• 67016E— 02 
•66938 E— 02 
• 66917E— 02 
• 66944E-02 
• 66958E— 02 
• 66855E— 02 

•  66583  E— 02 
• 661C7E— 02 
• 65694E— 02 
• 65039E— 02 
•64102 E-02 
•91956 E-02 
• 10229E— 01 
• 10879E— 01 
. 11067E— 01 
• 10650E— 01 
• 98697E— 02 
• 92270E— 02 
«  94411 E— 02 
.97799 E-02 
•10197E— 01 
•10433E— 01 


• 71449E+00 
•71484E+00 
• 71722E+00 
.72133  E+00 
• 72435E+0C 
• 72866E+Q0 
• 73516E+0C 
• 73633E+00 
•73632E+00 
•73622E+00 
•73618E+00 
• 72642E+00 
. 63302 E+00 

•  44  857  E+00 
• 73277E+00 
• 73037E+00 
•75163E+00 
.76570E+00 
•79749E+00 
• 81614E+00 

•  81649  E+00 
• 81987E+00 

•  83059  E+00 

•  84969  E+00 
. 85672E+00 
.92034 E+00 

•  10389  E+01 
•10549E+01 


• 69410E— 02 

•  695  85E— 02 
• 69014E— 02 
• 68298E— 02 
• 67976E— 02 
•67839E-02 
•67749E-02 
•67648E-02 
• 67474E— 02 
• 67285E-02 
•67127E— 02 
•6T006E-02 
• 66933E— 02 
• 66918E— 02 

•  669  47E— 02 
•66955E-02 

•  668  36E-02 
•66548E— 02 
•66143E— 02 
• 65637E— 02 
• 64963E-02 
•63850E-02 
• 94216E— 02 
•10304E— 01 
.109  25E— 01 

•  110  50E— 01 
.105  83E— 01 
•97858E-02 
•92096E-02 
. 94755E— 02 
. 98192E— 02 
• 10240E-01 
.104  26E— 01 


•71450E+00 
•71495E+00 
• 71760E+00 
•72164E+00 
• 72465E+00 

•  729  31E+00 
. 735  46E+00 
•73635E+00 
.73670E+00 

•  736  21E+00 
•73608E+00 
.711 41E+00 
• 601 73E+00 

•  415  34E+00 
•731 53E+00 
•73190E+00 
•75402E+00 

•  768  25E+00 
• 80033E+00 
•81663E+00 

•  816  59E  +  00 
.820  46E  +  00 
« 83238E+00 
• 85096E+00 

•  858  57E  +  00 
• 93251E+00 
•10456E+01 

•  105  32E  +  01 


•69448E-02 

•  69563E— 02 
•68931E— 02 

•  68  249E-02 
•67958E— 02 

•  678  29  E— 02 
• 67740E— 02 
• 67634E— 02 
• 67454E— 02 
. 67268E— 02 
• 67113E— 02 
• 66996E— 02 
• 66929E— 02 

•  66920E— 02 

•  66951 E— 02 
• 66950E— 02 
•66814E-02 
•66512E-02 
•66097E— 02 
• 65579E-02 
.64  884  E-02 

•  63  5  06E-02 
• 94962E-02 
. 10374E-01 
.  10966E-01 

•  11 026E— 01 

•  105 13E-01 

•  97  034E-02 

•  92  046E-02 
• 95079E— 02 
•98592E-02 

•  10281 E-01 

•  10420E— 01 


.71 451 E  +00 
•71507E+00 

•  71793E  +00 
.72195E+00 

•  72497E  +00 
.72991E+00 

•  73  571E  +  00 

•  73  637E  +00 

•  73  646E  +00 

•  73620E  +00 

•  73601E  +00 

•  67543E  +00 

•  57142E  +00 
• 48766E+00 

•  73041E  +00 
•73334E+00 
•75593E+00 

•  77118  E  +  00 

•  80300  E  +00 

•  81690E  +  00 

•  81675E  +00 
.  82113E  +00 

•  83430E  +00 

•  85 197E  +00 

•  86118E  +00 

•  94  534E  +00 

•  10509  E  +01 
.105 14E  +01 


• 69484E-02 
• 69531E— 02 

•  6  8849E— 02 
♦68205E—02 
. 67941 E— 02 
•67819E— 02 
•67732E-02 
•67619E-02 
.  67435E— 02 
•67251E-02 
•67099E— 02 
.66987E-02 

•  6  692  5E— 02 
.66922E-02 

•  6695  4E-02 
•66944E— 02 
.6  679  IE-02 
• 66474E-02 

•  6  605  IE— 02 
.  65519E— 02 
.64805E-02 
.  6  3 172  E-02 
.  9523  4E-02 
• 1 0444E— 01 
• 1 1 002E-01 
• 1 0997E— 01 
•1 0440E-01 

•  962  32E-02 

•  92 117E-02 

•  95  391E— 02 
• 98998E— 02 
• 1 03 19E-01 
• 1 041 8E-01 


.71453E+00 
.71 523E+00 
•71840E+00 

•  72227E  +  00 
•7253 1E+00 

•  73  065E  +  00 
•73591E+00 

•  73  63  8E  +  00 
• 73642E+00 

•  73  621E  +  00 
•73608E+00 

•  65405E  +  00 
•54597E+00 
.  OOOOOE  +  OO 

•  7288  4E  +  00 
.73491E+00 
• 75735E+00 

•  7743  8E  +  00 
•80549E+00 
. 8170  OE  +  OO 

•  8 1 69  8E  +  00 

•  8  218  8  E+00 
. 83630E  +  00 
• 85272E+00 
. 86468E+00 
.95860E+00 
.105476+01 
• 10497E+01 


•  695 17E-02 

•  694  91 E— 02 

•  687  67E-0  2 

•  681 6  5E-02 

•  679  2  6E— 0  2 

•  678  09E-02 
• 67723E-02 
• 67603E-02 

•  6741 5E— 02 

•  672  34E— 02 

•  670  86E— 02 

•  6  6978  E— 02 

•  6692  2E-02 
• 66924E— 02 
.669  5  7E-02 

•  669  37  E— 02 
• 66765E— 02 
.  664  3  6 E— 02 
• 66003E— 02 

•  654  5  6E— 02 
• 64725E-02 

•  632 1 2E— 02 
. 95725E-02 
• 10514E-01 
.11031 E-01 

•  1 0962E-01 
• 10364E— 0 1 
. 95464E-02 
• 92297E-02 

•  957  01E-02 
.994  09  E— 02 

•  103  54E-01 

•  1 04  26E-Q1 


•  7145  4E  +  00 

•  715  41E  +  00 

•  719  3  3E  +  00 
• 72258E+00 
•72568E+00 
• 73204E+00 

•  7360  5E  +  00 
• 73639E+00 

•  736  42E  +  00 
• 73622E+00 

•  7363  8E  +  00 
. 64668E+00 

•  5  24  8  5E  +  00 
. 663  77E  +  00 
• 72702E+00 

•  736  62E  +  00 

•  758  3  6E  +  00 
• 77776E+00 
• 80777E+00 
. 81698E+00 
• 81727E+00 
. 82274E+00 
.  838  38E  +  00 
.853  2  5E+00 
. 869  20E  +  00 
. 972  01E  +  00 
.10572E+01 
« 10480E+01 


• 69545E-Q2 
.69442E-02 

•  686  8  8E  — 02 
•68129E-02 
.  679 11E-02 
.57  800E  — 02 

•  677 14E— 02 
.67587E-02 
.67395E-02 
.672186-02 

•  67073E-Q2 
.66970E-02 
.  669  20E-02 
.66927E-02 

•  669  59E  — 0  2 
.669286-02 

•  667  39E-02 
• 66397E-02 
.  6595  5E-02 

•  65  392E— 02 

•  64645E-02 

•  64  395E-02 
. 96651E— 02 

•  10582E-01 
•11054E-01 
.10922E-01 

•  102  85E— 0 1 
• 94740E-02 

•  92  566E-02 
. 96016E— 02 

•  998  26E-02 

•  103  84E-01 

•  10443E-01 


•  71457E  +00 

•  71 563E  +  00 

•  71 896E  +  00 

•  72289E  +00 
•72607E+00 

•  73205E+00 

•  73  615E  +00 
•73638E+00 

•  736  37E  +  00 
•73623E+00 

•  73  679E  +0  0 

•  64  4  64  E  +0  0 

•  50685E  +  00 
•74115E+00 
• 72817E+00 
.  73  8  54E  +  00 
•75912E+00 
.781226+00 
.  809826  +  00 
.816886+00 

•  81761E  +  00 
• 82370E+00 
. 84048E+00 
.853626+00 
.874826+00 
.985256+00 
.105866+01 
.  10463E+01 


,695696-02 

•  693  8  5  E— 02 
• 68612E— 02 
.680976-02 
.57898 E-02 

•  6779 1 E-02 

•  677  05  E-02 
,  67569  E— 0  2 
• 67376E-02 

•  672  02  E-02 
. 67061 E— 02 

•  669  6  2  E  — 02 

•  669 18  E-02 

•  6 69  30  E-0 2 
. 66960 E— 02 

•  66917  E-0  2 
.  66710E— 02 

•  663  57  E-02 

•  6590  5  E— 02 

•  653  26E-02 
.  64  564E-02 

•  67690E-02 
.9789  2  E-02 

•  106  49  E-0 1 
.11071 E-01 

•  10876E-01 
.10205  E-0 1 
.  9407  5  E-02 

•  929  02  E-02 

•  963  4 3  E— 0 2 
.  1002  5  E— 0 1 
.  10408  E-01 
.10471 E-01 


• 7146GE+00 
.  715  8  8  E  +  00 
,  71974  E  +  00 

•  72319  E  +  00 

•  726  50E  +  00 

•  73298  E  +  00 
.73621 E+00 

•  736  3  4E  +  00 
. 7363  3  E  +  00 

•  736  24  E  +  00 
.  736  63  E  +  00 
•64494E+00 

•  491 2  2E  +  00 
.743  31E  +  00 

•  7  25  22E  +  00 

•  74075E  +  00 
•75985E+00 
.78467E+00 
.81163  E  +  00 
.  816  7  5  E  +  0  0 

•  817  98  E  +  00 
.  824  80  E  +  00 

•  842  56  E  +  00 

•  853  91E  +  00 
• 88160E+00 

•  998  01 E  +  00 
.105  9  0  E  +  0 1 
.10448  E  +  01 


281 

290 

.10433E+01 

.104196+01 

•10406E+01 

•  10  393E+01 

291 

300 

• 103046+01 

.102906+01 

.10277E+01 

. 10264  E+01 

301 

310 

•10168E+01 

.101546+01 

.101416+01 

.10128  E+01 

311 

320 

•10043E+01 

• 10032E+01 

• 10021E+01 

• 10011 E+01 

321 

330 

•99486E+00 

. 99414E+  00 

.993476+00 

• 99284E+00 

331 

340 

•98957E+00 

.989186+00 

.988786+00 

• 98837E+00 

341 

350 

.984906+00 

•98438E+00 

• 98388E+00 

•  98341 E+00 

351 

360 

•98147E+00 

.98132E+00 

.981166+00 

. 98098E+00 

361 

370 

•97734E+00 

•97639E+00 

.97532E+00 

.97416 E+00 

371 

380 

•96448E+00 

•96300E+00 

• 96152E+00 

• 96003E+00 

361 

390 

•94935E+00 

.947766+00 

• 94613E+00 

• 94439E+00 

391 

400 

•95402E+00 

•96595E+00 

• 97572E+00 

.981426+00 

401 

410 

•98386E+00 

.984676+00 

.985286+00 

• 98575E+00 

411 

420 

• 98905E+00 

.989306+00 

.989416+00 

• 98937E+00 

421 

430 

• 98622E+00 

.985616+00 

.984986+00 

.984326+00 

431 

440 

.97593 E+00 

.973826+00 

.971456+00 

.968836+00 

441 

450 

•94682E+00 

•94366E+00 

.940616+00 

•  93769  E+00 

451 

460 

•92166E+00 

.920106+00 

•91874E+00 

•  91760  E+00 

461 

470 

•91502E+00 

.9150^E+00 

.91512E+00 

•  91513  E+00 

471 

480 

.914126+00 

.913856+00 

• 91356E+00 

• 91323E+00 

481 

490 

•91120E+00 

.911426+00 

• 91195E+00 

.91 286 E+00 

491 

500 

.934256+00 

•  9394~'6+00 

. 94511E+00 

.95108 E+00 

501 

502 

.993276+00 

•10021E+01 

•  103  80E+01 
•10250E+01 
•10115E+01 

•  100  01E+01 
.992  276+00 
.98794E+00 
.982  996+00 
• 980746+00 
.97292E+00 

•  958536+00 
• 94250E+00 
.983416+00 
•986176+00 

•  989196+00 

•  983  60E+00 
•966006+00 

•  934906+00 

•  916  686+00 

•  915 10E+00 

•  912  87E+00 
•91424E+00 
.957  29E+00 


•  10368E  +01 
.  10236E  +01 

•  10102  E  +01 

•  99911E  +  00 
.  99173  E  +  00 
•98748E+00 
•98261E+00 

•  98042E  +00 

•  97162E  +  00 
•95702E+00 

•  94043E+00 
*98  3  21  E+00 

•  98662  E+00 
• 98888E+00 
• 98278E+00 

•  96299E  +00 

•  93227E  +  00 
•91598E+00 
. 91502  E+00 
•91248E+00 
♦91613E+00 
•96359E  +  00 


•10355E+01 

*  1 022  2E  +  01 
• 10089E+01 
.99818E+00 
.99124E+00 
. 98  700E+  00 
•98229E+00 
•98003E+00 
•97025E+00 
.9555 1E+0Q 
.93843E+00 
.  98238E  +  00 
.  98712E  +  00 
•988466+00 
.981 82E+00 
•95984E+00 
•92980E+00 

•  91 5  4  9E  +  00 
.91490E+00 
•91209E+00 
•91859E+00 
.969756+00 


. 103426+01 

•  102  09E  +  01 
• 10077E+01 
. 997  29E  +  00 
.990796+00 
.  98649E  +  00 
. 98203E+00 

•  979  5  3E  +  00 

•  968  8  5E  +  00 
. 95398E+00 
• 93721E+00 
. 981 94E  +  00 
.  98766E  +  00 
. 98796E+00 
• 98069E+00 

•  956  61E  +  00 
• 92750E+00 
. 91513E+00 
. 91474E+00 
. 91172E+00 
.921646+00 
.975556+00 


.  1O33QE  +  0L 
.10195E+01 
.100656+01 
•99643E+00 
.990376+00 
.98597E+00 
.98181E+0Q 
.978916+00 
.967416+00 
.952456+00 
■9383 4E +00 
•9821 9E +00 
.98819E+00 
•9874 OE +00 

•  97  9  3  5E  +00 
.95333E+00 

•  92  5  3  7E  +  00 
•91503E+00 
.91456E+00 
.91142E+00 
.92528E+00 

•  98  098E  +  00 


.103176+01 
.  10181E  +  01 
.  10054E  +  01 
• 99562E+00 

•  98996  E  +  00 
. 98544E+00 
•98163E+00 
.  97  8  1 9  E  +  0  0 
.96595E+00 
. 95090E+00 
. 943  8  4  E  +  0  0 
.  982  9  5  E  +  00 
•98857E+00 
.98682E+00 
•97777E+0Q 
•95006E+00 
. 92342E+QG 

•  914  99E  +  00 

•  91435  E  +  00 

•  91 1 22  E  +  00 
.929506+00 
.  98654E  +  30 


PRESSURE  HISTORY  AT  TEST  STATION  NO.  007 


1.78  METERS  FR  OIV’T-NOZZLE  EXIT  * 

as*  aiJitsisirnsiiititirtsiiiiii 3  1X13 


TINE  OF  SHOCK  ARRIVAL,  . 

SHOCK  OVERPRESSURE,  . 

SIDE-ON  OVERPRESSURE  IMPULSE, 
(EQUIVALENT  WEAPON  YIELD, 
POSITIVE-PHASE  DURATION,  .... 
(EOUIV.  FREE-FIELD  VALUE, 


TA  *  .3556E-02  SEC 

P50  -  .1976E  +  03  KP A  (28.7  PSD 

ISO  »  . 1272E+01  KPA-S 

VSO  =  . 27A6E-02  KT) 

PPD  =  .1651E-01  SEC 

PPF  »  . 2301E— 01  SEC) 


SHOCK  DYNAMIC  PRESSURE,  . 

(EQUIV.  FREE-FIELD  VALUE,  ... 

DYNAMIC-PRESSURE  IMPULSE,  . 

(EQUIVALENT  WEAPON  YIELD,  ... 

DYNAMIC-PHASE  DURATION,  . 

(EQUIV.  FREE-FIELD  VALUE,  .. 


QS  =  .1076E+03  KP A  (15.6  PSI) 

QF  «  .1078E  +  03  KP  A ) 

IQ  *=  .A805E+00  KPA-S 

WQ  =  •  1816E— 02  KT) 

QPD  -  .1652E-01  SEC 

QPF  -  . 1728E+00  SEC) 


TIME(S)  PSO(KPA)  ISO(KPA-S)  PDYN(KPA) 

aa****aa**rs*asa*a*aBr*sstaseea****2sessst** 


IDYN(KPA-S)  VFLO( M/S )  RH0(KG/M**3) 

taeasxt&a 33a**saa*tt*at»seaa**s#i** s  & 


m OOOOE+OO 
•3280E-03 
. 6123E-03 
. 8749E—03 
•1151E-02 
•1429E-02 
•1711E-02 
« 1994E— 02 
. 2281 E-02 
•2570E-02 
•2362E-02 
•3153E-02 
•3455E-02 
.3754E-02 
•4061E-02 

•  436^-02 
• 4680E— 02 
•4996E-02 
•5318E-02 

•  *^649E— 02 
• 5978E— 02 

•  631SE-02 
• 6658E— 02 
•7000E-02 

•  -f341E-02 
•7678E— 02 
• 8023E—0Z 

•  83“78E—02 
• 8710E— 02 
•9069E-02 
•9441E-02 
• 9787E-02 
• 1015H— 01 
•1051E— 01 
•1087E-01 
• 112 2 E— 01 
• 1157E— 01 
•1192E-01 


•OOOOE+OO 
•OOOOE+OO 
• OOOOE+OO 
•OOOOS+OO 
.OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
-.8504E+00 
•1824E+03 
• 1798E+03 
•1768E+03 
• 1737E+03 
• 1670E+03 
•1633E+03 
•1619E+03 
•1572E+03 
• 1519E+03 
• 1480E+03 
• 1338E+03 
•1258E+03 
•1134E+03 
•1033E+03 
• 9699E+02 
•9681E+02 
• 9720E+02 
• 9363E+02 
• 8952E+02 
•8856E+02 
• 8690E+02 
• 8307E+02 
•7892E+02 
•7335E+02 
• 6655E+02 


•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•3581E-01 
•9168E-01 
• 1461E+00 
• 2011E+00 
• 2549E+00 
• 3080E+00 
• 3616E+00 
• 4145E+00 
• 4669E+00 
• 5180E+00 
•5672E+00 
• 6123E+00 
• 6526E+00 
• 6899E+00 
•7253E+00 
•7573E+00 
•7923E+00 
• 8279E+00 
• 8595E+00 

•  8919E+00 
• 9233E+00 
• 9538E+00 
•9827E+00 

•  1009E+01 
•1034E+01 


•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•2500E-02 
• 9324E+02 
• 9085E+02 
• 8817E+02 
• 8546E+02 
• 7960E+02 
«  7646E+02 
•7E29E+02 
•7135E+02 
• 6705E+02 
• 6468E+02 
• 6034E+02 
•5452E+0 2 
•4779E+02 
• 4127E+02 
. 3700E+02 
• 3665E+02 
• 3651E+02 
• 3384E+02 
• 3103E+02 
• 3023E+02 
• 2905E+02 
• 2667E+02 
• 2422E+02 
• 2118E+02 
• 1774E+02 


• OOOOE+OO 
• OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
. OOOOE+OO 
. OOOOE+OO 
• OOOOE+OO 
•OOOOE+OO 
• OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
• OOOOE+OO 
• OOOOE+OO 
•  1802  E— 01 
•4655E-01 
•7380E-01 
•1011E+00 
•1272E+00 
•1522E+00 
• 1773E+00 
• 2016E+00 
•2250E+00 
• 2474E+00 
. 2689E+00 
• 2885E+00 
•3057E+00 
• 3210E+00 
• 3348E+00 
• 3470E+00 
•3602E+00 
• 3733E+00 
• 3845E+00 
• 3956E+00 
• 4063E+00 
• 4163E+00 
• 4254E+00 
• 4333E+00 
• 4400E+00 


•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
•OOOOE+OO 
— •  2066  E  +  01 
• 2798E+03 
•2771E+03 
• 2740E+03 
• 2708E+03 
•2636E+03 
• 2596E+03 

•  2581E  +  03 

•  2529E  +  03 
•2469E+03 
• 2425E+03 
• 2319E+03 
• 2167E+03 
• 2019E+03 

•  1892E+03 

•  1809E+  03 
•1803E+03 
•1804E+03 

•  1754E+  03 
• 1697E+03 
• 1681E+03 

•  1657E  +  03 
. 1603E+03 
• 1544E+03 

•  1464E+03 
• 1364E+03 


•  1 178  E  +01 

•  1178  E  +  01 
•1178E+01 
• 1178E+01 

•  1178  E  +01 

•  1178  E+01 

•  1 178  E  +  01 

•  1178  E+01 

•  1 178  E+01 

•  1178  E  +  01 

•  1178  E  +  01 

•  1178  E  +01 
« 1171 E+01 

•  2381 E+01 

•  2  366E+01 
• 2348E+01 

•  2331 E+01 

•  2291 E  +01 

•  2  269  E  +01 
•2260E+01 
•2231E+01 

•  2  200  E+01 
• 2200E+01 

•  2243  E+01 

•  2323  E+01 
• 2346E+01 

•  2305  E+01 

•  2261 E+01 
•2255E+01 

•  2244  E  +01 

•  2 199E  +  01 

•  2 155  E+01 
•2139E+01 

•  2118  E  +  01 

•  2077  E  +01 
• 2032E+01 
• 1976E+01 

•  1908  E  +01 


.1227E-01 
•1262E-01 
•1296E-01 
• 1330E— 01 
• 1364E— 01 
•1398E-01 
• 1433E— 01 
•1468E-01 
• 1504E— 01 
• 1539E—01 
. 1575E— 01 
• 1611E— 01 
• 1649 E— 01 
. 1689E-01 
• 1730E— 01 
. 17^1E-01 
•1812E-01 

•  1853E-01' 
• 1896E— 01 
•1937E— 01 
.19^6E-01 
. 2017E-01 
. 2060E-01 
• 2102E-01 
• 2143E— 01 
• 2187E-01 
. 2230E-01 
• 2272E— 01 
•2314E-01 
•2355E-01 
•2398E-01 
• 2442E-01 
• 2487E-01 
•2532E-01 
• 2577E-01 
. 2624E-0I 
.26~H.E-01 
• 2718E— 01 
•2766E-01 
•2815E-OT 
. 2863E-01 
•2911E-01 
•2959E-01 
• 3008E— 01 
. 3058E-01 
• 3109E— 01 
. 3159E-01 
• 3209E— 01 
• 3260E— 01 
• 3312E-01 
•3366E-01 
. 3420E-01 

•  34">5E-01 
• 3528E— 01 
• 3581E-01 
• 3633E— 01 
•3686E-01 
• 3741E-01 
. 3795E— 01 


. 5979E+02 
• 5403E+02 
• 4897E+02 
• 4389E+02 
• 3933E+02 
• 3619E+02 
• 3654E+02 
• 3752E+02 
• 3611E+02 
•3472E+02 
• 3371E+02 
•3196E+02 
• 2941E+02 
.260^E+02 
• 2225E+02 
•1817E+02 
• 1422E+02 
• 1060E+02 
•7283E+01 
• 4476E+01 
• 1940E+01 
-.6224E+00 

—  •  3033E+01 
. 2229E+02 
• 2451E+02 
•2367E+02 

•  2264  E+02 
•2160E+02 
• 2062E+02 
• 1971E+02 
•1902E+02 
•1345E+02 
• 1678E+02 
•1382E+02 
• 9996E+01 
•5S12E+01 

•  3569E+01 

•  3104*i+01 
-•1354E+00 

—  • 3886E+01 

—  . 6805E+01 
-.9174E+01 
— •  1147E+02 
— • 1368E+02 
— • 1606E+02 
— • 1810E+02 
-*•  1963  E  +  02 
—• 2126E+02 
— • 2298  E+02 
— • 2354E+02 
— • 2500E+02 
— • 2742E+02 
— • 2947E+02 
— . 3156  E  +  02 
— •  3390E+02 

—  •  3627E+02 
— • 3838E+02 

4157E+02 
—•  4383E+02 


•1056E+01 
• 1075E+01 
•1093E+01 
• 11 09E+01 
•1123E+01 
• 11 36E+01 
•1148E+01 
• 1161E+01 
•1175E+01 
•1187E+01 
•1200E+01 
• 1211E+01 
• 1223E+01 
•1234E+01 
•1244E+01 
•1252E+01 
• 1259E+01 
•1264E+01 
• 1268E+01 
•1270E+01 
• 1272E+01 
• 1272E+01 
• 1271E+01 
•1273E+01 
« 1283E+01 
•1293E+01 
•1303E+  01 

•  1313E+01 
•1321E+01 
• 1330E+01 
• 1338E+01 
•1346E+01 

•  13  54E+01 
•1361E+01 
• 1367E+01 
•1370E+01 
« 1372E+01 
•1374E+01 
• 1375E+01 
• 1374E+01 
* 1371E+01 
« 1367E+01 
• 1362E+01 
•1356E+01 
•1349E+01 
•1340E+01 
•1331E+01 
•1320E+01 
•1309E+01 
•1297E+01 

•  12  84E+01 
• 1270E+01 
•1254E+01 
• 1238E+01 
•1221E+01 
•1202E+01 
•1182E+01 
• 1160E+01 
• 1137E+01 


• 1459E+02 
• 1211E+02 
•1010E+02 
• 8275E+01 
• 6796E+01 
•5858E+01 
• 5927E+01 
• 6177E+01 
•5767E+01 
• 5379E+01 
• 5106E+01 
• 4652E+01 
• 4026E+01 
•3271E+01 
. 2500E+01 
•1787E+01 
• 1206E+01 
•7710E+00 
•4551E+00 
• 2505E+00 
•1169E+00 
• 3310E— 01 
•9499E-03 
• 2320E+0 1 
• 2686E+01 
«  2494E+01 
•2289E+01 
• 2095E+01 
. 1921E+01 
• 1750E+01 
• 1538E+01 
• 1323E+01 
•1220E+01 
• 1261E+01 
• 1341E+01 
• 1415E+01 
• 1795E+01 
• 2529E+01 
• 2722E+01 
• 2764E+C1 
• 2931E+01 
•3165E+01 
• 3366E+01 
• 3560E+01 
• 3680E+01 
• 3872E+01 
•4154E+01 
• 4382E+01 
• 4558E+01 
•5119E+01 
• 5388E+01 
• 5289E+01 
• 5287E+01 
• 5197E+01 
• 4973E+01 
• 4705E+01 
• 4337E+01 
• 3938E+01 
• 3656E+01 


• 4457E+00 
• 4503E+00 
• 4541E+00 
• 4572E+00 

•  4598  E+00 
•4619E+00 
• 4639E+00 
• 4661E+00 
. 4682E+00 
• 4702E+00 
• 4721E+00 
. 4739E+00 
• 4755E+00 
•47^0E+00 
.4781 E+00 
, 4790E+00 
• 4796E+00 
«  4800E+00 
• 4803E+00 
• 4804E+00 
• 4805E+00 
• 4805E+00 
. 4805E+00 
• 4807E+00 
.  4818  E  +  00 
«  4830 E  +  00 
. 4840E+00 
. 4849E+00 
. 4858E+00 
. 4865E+00 
• 4872E+00 
• 4879E+00 

•  4884  E  +  00 
. 4890E+00 
• 4896E+00 
«  4902E  +  00 
• 4909E+00 
. 4920E+00 
. 4932E+00 
• 4946E+00 
• 4960E+00 
• 4974E+00 
. 4990E+00 
. 5007E+  00 
• 5025E+00 
.  5044E+  00 

•  5064  E  +  00 
• 5086E+00 
• 5108E+00 

•  5133  E  +  00 
. 5162 E+00 
• 5191E+00 
• 5220E+00 
. 5248E+00 

•  5275  E+00 
. 5300E+00 
• 5324E+00 
• 5347E+00 
• 5367E+00 


•  1261E  +  03 

•  1170E+  03 
• 1087E+03 
•1002E+03 
• 9219E+02 

•  8658  E  +  02 

•  8725E  +  02 

•  8904E  +  02 
• 8647E+02 

•  8389E  +  02 

•  8200E  +  02 
•7869E+02 
•7377E+02 
• 6719E+02 

•  5945E+02 

•  5093  E+02 
• 4240E+02 

•  3431E+  02 
• 2666E+02 

•  1998E+  02 

•  1377E+  02 

•  7394  E+01 
• 1264E+01 

•  5758E  +  02 
• 6161E+02 
.  5955E  +  02 
•5724E+02 
. 5495E+02 
. 5279E+02 
• 5055E+02 

•  4751 E  +  02 

•  4418E+  02 
• 4268E+02 

•  4384E+  02 
•4580E+02 
• 4777E+02 

•  54  29E  +  02 
• 6467E+02 
• 6802E+02 
• 6967E+02 
• 7277E+02 
.  7659E+02 
• 8002E+02 

•  8340E  +  02 
.  8601 E  +  02 

•  8939E  +  02 
• 9364E+02 
• 9734E+02 
• 1005E+03 
• 1072E+03 

•  1109E+  03 

•  1112E+03 

•  1124  E+03 

•  1123E+03 

•  1119E  +  03 

•  1105E+  03 

•  1078E  +  03 
•1045E+03 
•1024E+03 


•  1 835  E+01 
•1769E+01 

•  1708  E+01 
• 1650E+01 
• 1599E+01 
« 1 563  E+01 
• 1557E+01 

•  1 558  E+01 

•  1 543  E  +01 

•  1 529  E+01 
• 1519E+01 
.1503 E+01 

•  1 480  E  +01 

•  1449  E  +01 
• 1414E+01 
• 1377E+01 

•  1 342E+01 
• 1310E+01 

•  1 280  E+01 
• 1255E+01 

•  1 233E+01 
• 1211E+01 

•  1 190  E+01 
. 1400E+01 

•  1  415  E  +01 
• 1406E+01 
•1397E+01 

•  1 388  E  +  01 
• 1 379E+01 

•  1 370  E+01 
•1362E+01 

•  1 356  E  +01 

•  1 339  E  +  01 
• 1312E+01 
. 1278  E  +  01 
•1240E+01 

•  1218  E  +01 

•  1 209  E  +  01 
. 1177E+01 

•  1 139E  +  01 

•  1 107E+01 
• 1079E+01 

•  1 051 E  +01 

•  1 024E+01 
• 9949E+00 
• 9691E+00 

•  9475  E+00 

•  9248  E  +  00 
•9021E+00 
• 8910E+00 

•  8758  E  +  00 

•  8553  E+00 

•  8  373E  +  00 

•  8171 E+00 
• 7942E+00 
• 7713E+00 

•  7468  E  +00 
. 7208  E  +00 
• 6977E+00 


FRENCH  BLAST  SIMULATOR  CASE:  MOD.A2-26  -  PLOT  1 

XSTAi  =0.445 

Vary  =  0.004  m3  Offset,  Ay  =  0.010 

P4/P1  ~  107.8  Grid  Clustered  About  Xclust  =  0.140 

Lref  =  6.640  m 


DENSITY  vs.  DISTANCE 
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SHOCK  TUBE  W/STEPPED  DRIVER 


CASE  =  MODEL-A2 
L— ref  —  6.640  m 
L— drv  =  0.884  m 
V— drv  =  0.004  m3 
L-dvn  =  5.461  m 
L— rwe  =  0.00  m 


X-sta  =  1.778  m 
P— drv  =  107.8  atm 
P-amb=  101.4  kPa 
T-amb=  299.8  K 
T4/T1  =1.59 


P— SO  = 

197.6 

kPa 

( 

28.66 

Psi) 

t— a  = 

3.556 

ms 

PPD  = 

0.017 

s 

< 

0.023 

3) 

I—so  = 

1.272 

kPa-s 

( 

0.003 

kT) 

Q-s  = 

107.6 

kPa 

I-dyn  = 

0.481 

kPa-s 

( 

o 

o 

o 

w 

kT) 

PRESSURE-TIME  HISTORY 
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